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AVERTISSEMENT. 



J’ai réuni dans cet ouvrage les divers écrits publiés 
par Gauss sur la Méthode des moindres carrés. 

L’illustre géomètre, que les sciences viennent de 
perdre, attachait une grande importance à cette partie 
de ses travaux, et la meilleure manière de combiner 
les observations était, à scs yeux, un des problèmes les 
plus importants de la philosophie naturelle. 

Gauss n’ignorait pas les critiques dont sa théorie 
a été l’objet; mais son opinion bien arrêtée était que 
les géomètres adopteraient entièrement ses idées lors- 
que ses Mémoires, aujourd’hui fort rares, sc seraient 
répandus davantage. C’est pour cela, sans doute, 
qu'il a bien voulu m’écrire qu’il me verrait avec le 
plus grand plaisir en publier la traduction. Ces 
Mémoires forment un Traité complet de la combinai- 
son des observations , qui n’exige ni commentaires ni 
annotations. I-es questions de priorité sur lesquelles 
se sont engagées des discussions assez vives, y sont 
traitées brièvement, mais de la manière la plus nette 
et la plus loyale. J’ai donc dû me borner au rôle de 
traducteur : c’était le seul qui fût utile, et le seul 
d’ailleurs que Gauss m’eût autorisé à prendre. 

i. HEirrnAiSD. 
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MÉTHODE DES MOINDRES CARRÉS. 


THÉORIE 


LA COMBINAISON DES OBSERVATIONS 

CCI EXPOSE AUX «BUES ERREURS, 


Par Cn.-FRv.ntRic GAUSS. 



PREMIÈRE PARTIE, 

PRÉSENTÉE A LA SOCIÉTÉ. ROYALÉ DE COTTINGl'E , LE l5 FÉVRIER l8ai. 


\. 

Quelque soin que l’ou apporte aux observations qui con- 
cernent la mesure des grandeurs physiques , elles sont for- 
cément soumises à des erreurs plus ou moins considérables. 
Ces erreurs, dans le plus grand nombre des cas, ne sont 
pas simples, mais découlent à la fois de plusieurs sources 
distinctes qu’il est bon de distinguer en deux classes. 

Certaines causes d’erreurs dépendent , pour chaque obser- 
vation, de circonstances variables et indépendantes du ré- 
sultat que l’on obtient : les erreurs qui en proviennent 
sont nommées irrégulières ou fortuites , et de même que les 
circonstances qui les produisent, leur valeur n’est pas sus- 
ceptible d’être soumise au calcul. Telles sont les erreurs 
qui naissent de l’imperfection de nos organes et toutes celles 
qui sont dues à des causes extérieures irrégulftres , comme, 
par exemple, les trépidations de l’air qui rendent là vision 
moins nette; quelques-unes des erreurs dues à l’imperfec- 

c. i 
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lion inévitable des meilleurs instruments appartiennent à 
la même catégorie. Nous citerons, par exemple, la rugosité 
de la partie intérieure du niveau, le defaut de rigidité 
absolue, etc. 

Il existe, au contraire, d’autres causes qui, dans toutes 
les observations de même nature, produisent une erreur 
identique, ou dépendant de circonstances essentiellement 
liées au résultat de l’observation. Nous appellerons les er- 
reurs de cette catégorie, deserreurs constantes ou régulières. 

Il est du reste évident que cette distinction est jusqu’à un 
certain point relative et dépend du sens plus ou moins 
large que l’on veut attacher à l’idée d’observations de même 
nature. Par exemple, si l’on répète indéfiniment la mesure 
d’un même angle, les erreurs provenant d’une division 
imparfaite du limbe appartiendront à la classe des erreurs 
constantes. Si, au contraire, on mesure successivement 
plusieurs angles ditférents, les erreurs dues à l'imperfec- 
tion de la division seront regardées comme fortuites tant 
que l’on n’aura pas formé la table des erreurs relatives à 
chaque division. . 

2 . 

Nous excluons de nos recherches la considération des 
erreurs régulières. C’est à l’observateur qu’il appartient 
de rechercher avec soin les causes qui peuvent produire 
une erreur constante , pour les écarter s’il est possible, ou 
tout au moins apprécier leur effet, afin de le corriger sur 
chaque observation, qui donnera alors le même résultat 
que si la cause constante n’avait pas existé. Il en est tout 
autrement des errfeurs irrégulières: celles-là, par leur na- 
ture, se refusent à tout calcul, et il faut bien les tolérer dans 
les observatiims. On peut cependant, par une combinaison 
habile des résultats, réduire autant que possible leur in- 
fluence. C’est à cette question importante que sont consa- 
crées les recherches suivantes. 
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Les erreurs qui , tlans des observations de même espèce, 
proviennent d’une cause simple et déterminée se trouvent 
renfermées entre certaines limites que l’on pourrait sans 
aucun doute assigner, si la nature de cette cause était 
elle-même parfaitement connue. Dans la plupart des cas, 
toutes les erreurs comprises entre ces limites extrêmes 
doivent être regardées comme possibles. Une connais- 
sance approfondie de chaque cause apprendrait si toutes 
ces erreurs ont une facilité égale ou inégale, et, dans le 
second cas, quelle est la probabilité relative de chacune 
d’elles. La même remarque s’applique à l’erreur totale qui 
provient de la réunion de plusieurs erreurs simples. Cette 
erreur sera , elle aussi , renfermée entre deux limites dont 
l’une sera la somme des limites supérieures, l’autre celle 
des limites inférieures, correspondant aux erreurs simples. 
Toutes les erreurs comprises entre ces limites seront possi- 
bles, et chacune pourra résulter, d’une infinité de manières, 
de valeurs convenables attribuées aux erreurs partielles. On 
comprend néanmoins , en écartant les difficultés purement 
analytiques , qu’il y a possibilité d’apprécier la probabilité 
plus ou moins grandede chaque résultat, si l’on suppose con- 
nues les probabilités relatives à chacune des causes simples. 

Certaines causes pourtant produisent des erreurs qui ne 
peuvent pas varier suivant une loi continue, mais qui, au 
contraire, sont susceptibles d’un nombre fini de valeurs: 
nous pouvons citer, comme exemple, les erreurs qui pro- 
viennent de la division imparfaite des instruments (si tou- 
tefois on veut les classer parmi les erreurs fortuites), car 
le nombre des divisions, dans un instrument donné, est 
essentiellement fini. 11 est clair néanmoins que, si toutes les 
causes qui concourent à produire l’erreur totale ne sont pas 
supposées dans ce cas, leur somme formera une série sou- 
mise à la loi de continuité, ou, tout au moins, plusieurs 
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séries distinctes , s’il arrive qu'en plaçant par ordre de gran- 
deur toutes les valeurs possibles des erreurs discontinues, 
la différence entre deux termes consécutifs de la série soit 
moindre que la différence entre les limites extrêmes des 
erreurs soumises à la loi de continuité. Dans la pratique, un 
pareil cas ne se présentera presque jamais; il supposerait 
des défauts trop grossiers dans la construction de l’instru- 
ment. 

4 . 

Désignons par la notation q> (x) la facilité relative d’une 
erreur x : on doit entendre par là , à cause de la conti- 
nuité des erreurs, que <p (x) dx est la probabilité que l’er- 
reur soit comprise entre les limites x et x dx. Il n’est pas 
possible , en général, d’assigner la forme de la fonction œ, 
et l’on peut même affirmer que cette fonction ne sera jamais 
connue dans la pratique. On peut néanmoins établir plu- 
sieurs caractères généraux qu’elle doit nécessairement pré- 
senter : y (x) est évidemment une fonction discontinue; 
elle s’annule pour toutes les valeurs de x non comprises 
entre les erreurs extrêmes. Pour toute valeur comprise entre 
ces limites, la fonction est positive (en excluant le cas in- 
diqué à la fin du paragraphe précédent); dans la plupart 
des cas , les erreurs égales et de signes contraires seront 
également probables , et l’on aura 

?(■*)=?(— *)■ 

Enfin, comme les petites erreurs sont plus facilement com- 
mises que les grandes, <f (x) sera en général maximum pour 
x — o et diminuera sans cesse lorsque x croîtra. 

I/intégralc 

h 

rf (x) lix 

exprime la probabilité pour que l’erreur, eucore inconnue, 
tombe entre les limites a et b. On en conclut que la valeur 
de cette intégrale prise entre les limites extrêmes des erreurs 
possibles sera toujours égale à l’unité. Et comme <f (x) est 
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nulle pour les valeurs non comprises entre ces limites, on 
peut dire , dans tous les cas , que 


/: 


t f[x)dx — - i . 


D. 


Considérons l'intégrale 


/: 


xif(x)dx 


et représentons sa valeur par k. Si les causes d’erreur sont 
telles, qu’il n’y ait aucune raison pour que deiix erreurs 
égales et de signes contraires aient des facilités inégales, on 
aura 

?(*) = ?(—•*•)> 
et , par suite, • . 

k = o. 

Nous en conclurons que, si k ne s’évanouit pas et a, par 
exemple, une valeur positive, il existe nécessairement une 
cause d’erreur qui produit uniquement des erreurs positives 
ou qui , tout au moins , les produit plus facilement que les 
erreurs négatives. Cette quantité A, qui est la moyenne de 
toutes les erreurs possibles, ou encore la valeur moyenne 
de x , peut être désignée commodément sous le nom de partie 
constante de l’erreur. Du reste , on prouve facilement que 
la partie constante de l’erreur totale est la somme des par- 
ties constantes des erreurs simples qui la composent. 

Si la quantité k est supposée connue et qu’on la retranche 
du résultat de chaque observation , en désignant par x' l’er- 
reur de l’observation ainsi corrigée, et la probabilité cor- 
respondante par cf' ( x ' ), on aura 

x' = x — k, f' (x') — <f(x) 

et, par suite, 

/»oû . . /» 0 O; f * 20 

1 x'f'[x I xtf(x)dx — k I tf(x)dx — k-~k = o-, 

a/— -CO a J — 00 */ — 00 
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en sorte que les erreurs des observations corrigées n’ont pas 
de partie constante. Ce qui, du reste, semble évident à 
priori. 

6 . 


La valeur de l’intégrale 




e’est-à-dire la valeur moyenne de x , fait connaître l’exis- 
tence ou la non-existence d’une erreur constante, ainsi que 
la valeur de celte erreur; de môme l’intégrale 



( x ) -rfr , 


c’est-à-dire la valeur moyenne de a:*, parait très-propre à 
définir et à mesurer, d’une manière générale, l’incertitude 
d’un système d’observations ; de telle sorte qu’entre deux 
systèmes d’observations inégalement précises, on devra re- 
garder comme préférable celui qui donne à l’intégrale 



une moindre valeur. Si l’on objecte que cette convention 
est arbitraire et ne semble pas nécessaire, nous en conve- 
nons volontiers. La question qui nous occupe a, dans sa 
nature même, quelque chose de vague et ne peut être bien 
précisée que par un principe jusqu’à un certain point arbi- 
traire. La détermination d’une grandeur par l’observation 
peut se comparer, avec quelque justesse, à un jeu dans lequel 
il y aurait une perte^ à craindre et aucun gain à espérer : 
chaque erreur commise étant assimilée à une perte que l’on 
fait, la crainte relative à un pareil jeu doit s’exprimer par la 
perte probable, c’est-à-dire par la somme des produits des 
diverses perles possibles par leurs probabilités respectives. 
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Mais quelle perte doit-on assimilera une erreur détermi- 
née? C’est ce qui n’est pas clair en soi; cette détermination 
dépend en partie de notre volonté. Il est évident, d'abord, que 
la perte ne doit pas être regardée comme proportionnelle à 
l’erreur commise; car, dans cette hypothèse, une erreur 
positive représentant une perte, l’erreur négative devrait 
être regardée comme un gain : la grandeur de la perte doit, 
au contraire, s’évaluer par une fonction de l’erreur dont la 
valeur soit toujours positive. Parmi le nombre infini de 
fonctions qui remplissent cette condition , il semble natu- 
rel de choisir la plus simple, qui est, sans contredit, le 
carré de l’erreur, et, de cette manière, nous sommes con- 
duit au principe proposé plus haut. 

Laplace a considéré la question d'une manière analogue, 
mais en adoptant, pour mesure delà perte, l’erreur elle- 
même prise positivement. Cette hypothèse, si nous ne nous 
faisons pas illusion , n’est pas moins arbitraire que la nôtre : 
faut-il, en effet, regarder une erreur double comme plus 
ou moins regrettable qu’une erreur simple répétée deux 
fois, et faut-il, par suite, lui assigner une importance 
double ou plus que double? C’est une question qui n’est 
pas claire, et sur laquelle les arguments malhématiques 
n’ont aucune prise ; chacun doit la résoudre à son gré. 
Un ne peut nier pourtant que l’hypothèse de Laplace ne 
s’écarte de la loi de continuité et ne soit, par conséquent, 
moins propre à une élude analytique ; la nôtre, au contraire, 
se recommande par la généralité et la simplicité de ses con- 
séquencès. ... . 

7 ; - 


Posons, en conservant les notations précédentes, 

X « • 

y (jt) x‘ dx — ni 1 : 


nous appellerons m V erreur moyenne à craindre ou, sim- 
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plement, Y erreur moyenne des observations considérées. 
Nous ne limitons pas, du reste, cette dénomination au ré- 
sultat immédiat des observations , mais nous l’étendons, au 
contraire, à toute grandeur qui peut s’en déduire d’une 
manière quelconque. Il faut bien se garder de confondre 
celte erreur moyenne avec la moyenne arithmétique des 
erreurs , dont il est question dans l’art. 5. 

Si nous comparons plusieurs systèmes d’observations ou 
plusieurs grandeurs résultant d'observations auxquelles on 
n’accorde pas la même précision , nous regarderons leur 
poids relatif comme inversement proportionnel à /«’, et leur 
précision comme inversement proportionnelle à m. Afin de 
pouvoir représenter les poids par des nombres , on devra 
prendre, pour unité, le poids d’un certain système d’obser- 
vations arbitrairement choisi. 

i ' • . . 

8 . 


Si les erreurs des observations ont une partie constante, 
en la retranchant de chaque résultat obtenu , l’erreur 
moyenne diminue, le poids et la précision augmentent. En 
conservant les notations de l’art. 5, et désignant par m' 
l’erreur moyenne des observations corrigées , on aura en 
effet 

m ,7 = f x' 1 y'(x') djt? = Ç (x — Æ)’<f> (x) dx 

J — OO J — 00 

X OO /» X /»« 

x J <p(x)rfx — 2 A I x <p (x) dx + A* I f(x)dx 
-00 J — oo J — 00 


— m 1 — 2 k 7 -f- A* = m 1 


AK 


Si , au lieu de retrancher de chaque observation le nom- 
bre h , on retranchait un autre nombre /, le carré de 
l’erreur moyenne deviendrait 

m» — 2 Al + P = m' 1 + [l — k)\ 
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Soient X un coefficient déterminé et p. la valeur de l'in- 
tégrale 


X 


/ m 


(f(x)dx, 


— im 

fj. sera la probabilité que l’erreur d’une certaine observa- 
tion soit moindre que Xw en valeur absolue; i — u sera, 
au contraire, la probabilité que cette erreur surpasse Xm. 

* j ■ t 

Si, pour p = -> Xm a la valeur p, il y aura probabilités 

égales pour que l’erreur soit plus petite ou plus grande 
que p : p pourra donc être appelé V erreur probablç. La re- 
lation qui existe entre X et p dépend de la nature de la 
fonction <f, qui est inconnue dans la plupart des cas. Il est 
intéressant d’étudier cette relation dans quelques cas parti- 
culiers. 

I. Si les limites extrêmes des erreurs possibles sont 
-I- a et — a, et si, entre ces limites, toutes les erreurs 
sont également probables , la fonction <p (x) sera constante 

entre ces mêmes limites, et, par conséquent, égale à — ' • 
Par suite, on aura „ 

0 . i 4 ^ 

m = a \v ^y/r 

tant que X sera inférieur ou égal à y ^3 ; enfin , 

p = m = 0,8660254»», 

et la probabilité pour que l’erreur ne surpasse pas L’erreur 
moyenne est 

0 , 5 ^ 3503 . 

II. Si — a et -+-« sont encore les limites des erreurs 
possibles , si l'on suppose , de plus , que la probabilité de 
ces mêmes erreurs aille eadécroissanl à partir de l’erreur o 
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comme les termes d’une progression arithmétique, on aura 

. a — x 
T(*) = -^r- 

pour les valeurs de x comprises entre o et -f- a , et 

a 4 - x 


f(*) : 


a 2 


pour les valeurs comprises entre o et — a : de là on déduit 

m=a \fv r=^V3 s v * 

tant cjue / est compris entre o et y'6’ ; 

\ — \/6 — — 6 fi, 

tant que (i est compris entre o et i; et, enfin, 

p = m (\/6 — y 3 ) = 0,7 ■ 74389 w. 

Dans ce cas, la probabilité que j’erreur restera inférieure 
à l’erreur moyenne sera 

y/ï - g ~ 0,6498299 ■ 

III. Si nous supposons la fonction <p (x) proportionnelle 

JT* 

À e h' [ce qui, eu réalité, n’est vrai qu’approximative- 
menl (*)], elle devra être égale à 


T (*) = 


on en conclut 


h y jt 


= Vi 


(voir Disquisitiones generales circa seriem infinilant, 
art. 28, Mémoires de GoUingue , tome II). 

(*) II faut se reporter, pour comprendre cette remarque, à un chapitre 
du Theoria Motus Corporum cœlestium , dans lequel M. Gauss montre que 
cette loi de probabilités est la plus vraisemblable que Ton puisse adopter. 
A la fin du volume nous reproduisons ce chapitre, dans lequel l'illustre au- 
leur a fait connaître pour la première fois la méthode des moindres carres. 

» J. B. 
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Si l’on désigne par &z la valeur de l’intégrale 


e ‘ dz, 


U 


Le tableau suivant donne quelques valeurs de cette quan- 
tité : 


X 

V 1* 

0,6744897 

o ,5 

0 , 84 i 6 ai 3 

0,6 

1 ,0000000 

0 ,6826895 

1 ,o 364334 

°>7 

1 ,2815517 ' 

0,8 

1 ,6448537 

0,9 

2,5758293 

°>99 

3,2918301 

°>999 

3 ,8905940 

°.9999 

CO 

1 


10 . 

Quoique la relation qui lie X à p dépende de la nature 
de la fonction y, on peut cependant établir quelques résul- 
tats généraux , qui s’appliquent à tous les.cas dans lesquels 
cette fonction ne sera pas croissante avec la valeur absolue 
de la variable x ; alors on aura les théorèmes suivants : 

X ne dépassera pas u V3 toutes les fois que p sera infé- 


rieur a - : 


X ne dépassera pas 
a 

passera 


3 s /~ 


toutes les fois que p sur- 
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Lorsque fi=-r» les deux limites coïncident et X ne peut 

T J 

pas être supérieur à y/j" 

Pour démontrer ce théorème remarquable, représentons 
par j la valeur de l’intégrale 


£ 


(p(z)rfi; 


alors y sera la probabilité pour qu’une erreur soit com- 
prise entre — x et -+-#• 

Posons 

* =+(r). {y) = fUHri 'iY(y) = Y(r)‘ , J r > 

on aura * 

. +(»)=. = (-T ) 5 

on en conclut, en ayant égard aux hypothèses qui ont ,été‘ 
faites, que, depuis ^ = 0 jusqu’à y= 1, ip' (y) est tou- 
jours croissant, ou du moins 11’ est pas décroissant, et que, 
par suite., <p" ( y) est toujours positif, ou du moins n’est 
pas négatif. Or nous avons 

d .y S/' (y) = -p’ (y) dy y (y) dy , 

par suite, 

yŸ(r) — -p{y)=f r yY(y)dy, 

\ *o 

y 'Y {j) — 'P (y) a *l° nc une valeur constamment positi\e, 
ou du moins cette expression ne sera jamais négative. 


Il suit de là que I ■ 


+ (r) 


sera toujours positif et inoin- 


/*'(/) 

dre que l’unité. Soit f la valeur de cette différence pour 
y = pj à cause de ’| (u) = X/n, on a 


fszt — 
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d’où l’on conclut ' 

' f(p). 


( «3 ) 


Cela posé, considérons la fonction 

u 

que nous désignerons par F (^), et posons 

' d.Y{ r )=r{y)d r - J 

on aura évidemment 

F (f*) = = 

■ rw =(T%; =+ ' w v. . .. 

Or, puisque ifi' ( y ) croit continuellement (ou du moins ne 
décroît pas, car c’est ainsi qu’on doit toujours l’entendre) 
lorsque y croît, et que, d’un autre côté, F'(^) est con- 
stant, la différence 

y(»-rtr)- ^ w r TW l 

• ' \ 

sera positivepour toutes les valeurs de y pl us grandes que p, 
et négative pour les valeurs de y plus petites que p. On en 
.conclut que la différence <p (y) — F (jp) est toujours posi- 
tive, et, par suite, ^ (y) sera certainement plus grand que 
F (_^) en valeur absolue, tant que la fonction F(_y) sera 
positive, c’est-à-dire depuis y = pf jusqu’à y = i. La 
valeur de l’intégrale 


f\ 

J U./ 




sera donc inférieure à celle de l’intégrale 

f [+(/)]’<*/, 

■W 

et à fortiori moindre que 

• ». .4 

f [+(rff‘Ô r f 

. , . J O 
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c’est-à-dire moindre que ni*. Or la première de ces inté- 
grales a pour valeur 


donc 


— f if ) 1 

" 3 ^( 1 -/)* 5 


3f*’(i — /Y 
(«-(*/) 


f désignant, on ne l’a pas oublié, un nombre compris 
entré o et î . 

Si nous considérons f comme variable , la fraction 


3 

(»— f */) 1 ’ 


aura pour différentielle 



Vil -zl). 


a — 3f» 4- f */) df : 


cette fraction sera donc continuellement décroissante lors- 


2 

que croîtra de o à i, et que l’on aura, en outre, £<i ; 

sa valeur maximum correspondra à f = o et sera égale à 
3p% de sorte que, dans ce cas, le coefficient X sera cer- 
tainement inférieur, ou du moins ne sera pas supérieur 
à p Ce qu’il fallait démontrer. 

Lorsqu’au contraire p est plus grand que la valeur 
de la fonction sera maximum lorsque 


* — 3p -+- p/= o, 


c’est-à-dire lorsque 


f — 3 




et cette valeur maximum sera 

4 , 

- 9 (* — f*) ' 
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par conséquent, dans ce cas, le coefficient X n’est pas plus 
grand que 


3 — F 

Faisons, par exemple, 


> comme nous l’avions annoncé. 


pas surpasser V* c’est-à-dire que l'er- 


alors X ne peut 

reùr probable ne peut pas surpasser 0 , 86602:14 n* , à la- 
quelle elle devient égale dans le premier cas examiné 
(art. 9) : on conclut facilement de notre théorème que p 

n’est pas moindre que Xt^^, tant fl 116 ^ est moindre . 
que et qu’au contraire, il n’est pas inférieur à 


1 — lorsque X est plus grand que 


Vr 


L’intégrale 


r 


U. 


x 1 ij) (.r ) dx , 


se présentant dans plusieurs problèmes que nous aurons à 
traiter, il ne sera pas inutile de l’évaluer dans quelques cas 
particuliers. 

Posons 


£ 


I. Lorsque 


x' <f (x) dx =x n' 


?(*) = ■ 


pour les valeurs de x comprises entre • — a et -+- a , on a, - 


n , __ “ _ g , 
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11. Lorsque 


« (*) 


a 7 


(II e cas , art. 9), x étant encore compris entre — a et -t- a ; , 
on a 


i . 12 , 

"'=75° =-5 m 


III. Dans le troisième cas, lorsque 


?(*) = 


h' 


h v/7 


on trouvera , d’après les résultats obtenus dans le Mémoire 
cité plus haut, 

n' = 7 A* = 3 m'. 

4 

On peut d’ailleurs démontrer qu’en restant dans les liypo- 

, . «* * 
thèses admises au paragraphe precedent, le rapport — 

. . q 

n’est jamais inferieur a g- 

12 . 

Désignons par X, x', x", etc., les erreurs commises 
dans des observations de même espèce, et supposons que 
ces erreurs soient indépendantes les unes des autres. Soit, 
comme plus haut, <p(x) la probabilité relative de l’erreur 
x ; considérons une fonction rationnelle y, des variables 
.r, x', x", etc. 

L’intégrale multiple 

(i) j' ?(*)?(*') • • dxdx’dx ..., 

étendue à toutes les valeurs des variables X, x', x", etc., 
pour lesquelles la valeur de y tombe entre les limites don- 
nées o, n- représente la probabilité que la valeur de y soit 
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comprise entre o et y. Or cette intégrale est évidem- 
ment une fonction de y. Représentons sa différentielle par 
(y) dy , de sorte que l’intégrale considérée soit égale à j 


f 4> 

t/O 


et que, par conséquent, '('.(y;) représente la probabilité 
relative d’une valeur quelconque de y. x pouvant être re- 
gardé comme une fonction des variables y, x', x", etc., 
que nous désignerons par f( y, x', x 1 ', . . .), l’intégrale (1) 
prendra la forme 




? [/', » r '> )] — — — ? (*' ) T (*" ) ■ • (, J , , 


où J doit varier depuis y — o jusqu’à y— y, et les au- 
tres variables reçoivent toutes les valeurs pour lesquelles 
J'{y,x',x"...) est réelle. 

On aura donc 

4- (r)=Jt [/(/i*'»*" • • •)] ■■ t ( r '; ? K)- Mit** 

l’intégration , dans laquelle y doit être regardé comme une 
' constante, s’étendant à toutes les valeurs des variables 
x', x", etc., pour lesquelles /(y, x 1 , x",...) prend une 
' valeur réelle. 

13 . 

h 

L’intégration précédente exigerait, il est vrai, la con- 
naissance de la fonction <p, qui est inconnue dans la plupart 
des cas. Lors même que celte fonction serait connue, le 
. calcul surpasserait, le plus souvent, les forces de l’ana- 
lyse. Dès lors il sera impossible d’obtenir la probabilité de 
chacune des valeurs de y ; mais il en sera autrement si l’ou 
désire seulement la valeur moyenne dey, qui sera donnée 


t 
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par l’intégrale 


( >8) 


J y ÿ(r) d y, 

étendue à toutes les valeurs possibles de y. 

Si, par la nature de la fonction, ou à cause des limites 
imposées à x,x', x",elc.,y n’est pas susceptible de recevoir - 
toutes les valeurs , on devra supposer que t Ji ( y) s’annule , . 
pour toutes les valeurs que y 11e peut atteindre, et l’on 
pourra alors étendre l’intégration de — 00 à -t- 00 . 

Mais l’intégrale 

/ y-Hy)dy, 

prise entre des limites déterminées r, et est égale à 


?(•*') f{x").,.dydx'dx"..., 


prise depuis j - = n jusqu’à y = W et étendue à toutes les 
valeurs des variables x', x", etc., pour lesquelles^ est réelle. 
Cette intégrale est égale, par conséquent, à l’intégrale V, 


f 


yi( x ) t( x ') ?( ■*")<•• • dxdx' dx" . . . , 


dans laquelle y sera exprimé en fonction de x, x\ x", etc., 
la sommation s’étendant à toutes les valeurs des variables 
qui laissent y compris entre n et yi'. D’après cela , l’inté- 
grale 

/■“ 

y^{y) d r 


£ 


peut se mettre sous la forme 

y y (x) tf(x') <p ( x ”) . . . dxdx’ dx" . . , 


£ 


l’intégration s’étendant à toutes les valeurs réelles de 
x, x', x", c’est-à-dire depuis x = — ao à x = + 00 , 
x' = — 00 à x' = -+- ac , etc. 
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U. 

Si la fonction y se réduit à une somme de termes de la 
forme 


kx^xï x"'t . 


la valeur de l’intégrale 


/ 


y\[y)<iy. 


étendue à toutes les valeurs de y, c’est-à-dire la valeur 
moyenne de y, sera égale à une somme de termes de la forme 


c est-à-dire que la valeur 'moyenne de y est égale à une 
somme de termes déduits de ceux memes qui composent y, 
en y remplaçant x ,/3 , x">, - etc., par leurs valeurs - 
moyennes. La démonstration de ce théorème important * 
pourrait facilement se déduire d’autres considérations. 


15 .' 

Appliquons le théorème précédent au cas où l'on a 

X* ■+■ X 11 ■+■ x" J -+- . . . 

CT — — ’ 

ct désignant le nombre des termes du numérateur. 

On trouve tout de suite que la valeur moyenne de y est égale 
à rrf (la lettre m ayant toujours la signification de l’art. 5). 
La véritable valeur dey peut être inférieure ou supérieure 
à sa moyenne, de même que la vraie valeur de a: 8 peut, dans 
chaque cas , être inférieure ou supérieure à m*; mais la pro- 
babilité pour que la valeur fortuite de y ne diffère pas sen- 
siblement de m 8 , s’approchera sans cesse de la certitude à 
mesure que a deviendra plus grand. Pour le montrer plus 
clairement, comme il est impossible de chercher exactement 
cette probabilité, nous chercherons l’erreur moyenne à' 

' f 2 . 
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craindre en faisant^ = m}. D’après ce qui a etc dit (art. 6), 
rctte erreur sera la racine carrée de la moyenne de la fonc- 
tion 

f x" -+- x ' 2 -f- x" 




.Pour la trouver, il suffit d’observer que la valeur moyentic 
d’un terme tel que — est égale à (n ayant la mèmesigni- 
fication que dans l’art. H ), et que la valeur moyenne d’un 

2 X‘ X * 2 /M* 

terme tel que — — est égale à — — ; par conséquent, la 

valeur moyenne de cette fonction sera 
n' — m' 


De là nous concluons que si le nombre des erreurs irré- 
gulières est suffisamment grand , la valeur de m sera repré- 
sentée, avec une grande certitude, par la formule 


— y/ X ’ + - C " ~ 


et l’erreur moyenne h craindre dans la détermination du 
carré de m , sera égale à 


V- 


Comme cette dernière formule contient la quantité n, 
si l’on veut seulement, se faire une idée du degré de préci- 
sion de cette détermination, il suffira d’adopter pour la 
fonction <f une hypothèse particûlière. 

Si nous prenons , par exemple , la troisième hypothèse des 

art. 9 et H, cette erreur sera égale à m* \Æ Si on le 

préfère, on pourra obtenir une valeur approchée de n* au 
moyen des erreurs elles-mêmes, à l’aide de la formule 

x* x' 1 -f- x"' . 
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Un peut affirmer généralement qu'une précision deux fois 
plus grande dans cette détermination exigera un nombre 
d’erreurs quadruple, c’est-à-dire que le poids de la déter- 
mination est proportionnel au nombre a. 

On verrait de la même manière que si les erreurs des 
observations renferment une partie constante, on déduira de 
leur moyenne arithmétique une valeur delà partie constante, 
et cette valeur sera d'autant plus approchée que Je nombre 
des erreurs sera plus grand. Dans cette détermination, 


, . , * , /ni 1 — k 1 

1 erreur moyenne a craindre sera représentée par i / » 

h désignant la partie constante, et in l’erreur moyenne des 
observations non corrigées de leur erreur constante. Elle 

m 

sera représentée simplement par , si ni représente. l’er- 
reur moyenne des observations corrigées de la partie con- 
stante ( voyez art. 8). 


16 . 


Dans les art. 12 à lo nous avons supposé que les 
erreurs .r, x', x", etc., appartenaient au même genre d’ob- 
servations, de sorte que la probabilité de chacune de ces 
erreurs était représentée par la même fonction. Mais il est 
évident que les principes généraux exposés dans les ar- 
ticles 12 à 14, peuvent facilement s’appliquer au cas plus 
général où les probabilités des erreurs x , x\ .t", etc., sont 
représentées par des fonctions dillércntcs 


.*?£*), ?'(•*')» ?’(■*")> ■ • ■ » 

c’est-à-dire lorsque ces erreurs appartiennent à des obser- 
vations qui n’ont pas le même degré de précision. Suppo- 
sons que a? désigne l’erreur d’une observation dont l’erreur 
moyenne à craindre soit m\ x x", etc., celles d’autres 
observations dont les erreurs moyennes à craindre soient 
respectivement m', m'\ etc, alors la valeur moyenne de 
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la somme 


( 23 ) . 


sera 


x 3 •+• x ' 1 -+- x " 3 -+- . . . 
m* -+- m' 7 + m" 7 -4- . . . . 


Or, si l’on sait, par ailleurs, que les quantités le., 

sont respectivement proportionnelles aux nombres, i , p', 
u", etc., la valeur moyenne de l’expression 

x 3 -+-x' 3 -i- x" 3 4- ... 
i -+- fï' 3 -t- fi" 3 -h , . . . 

sera égale à m*. Mais si nous adoptons pour ni 3 la valeur que 
prendra celte expression, en y substituant les erreurs 
x,x', x ", etc., telles que le hasard les offrira, l’erreur 
moyenne qui aflcctc cette détermination sera, d’après l’ar- 
ticle précédent, 

-+- «'* - 4 - «* 4 — m * — /»'* — »<"• — ... 

, + ^ + p», + ... ' 

où n', n", etc., ont la même signification, par rapport à la 
seconde et à la troisième observation , que n par rapport à la 
première; et si l’on peut supposer les nombres etc., 

proportionnels à m , m' , ni", etc., celte erreur moyenne à 
craindre sera égale à. 


-• 


V^( n‘ — m*) fi ft' 1 -4- fi" 4 
i + fi' 3 -t- f»" 3 


niais celte manière de déterminer une valeur approchée 
de m n’est pas la plus avantageuse. 

Considérons l’expression plus générale 

x 3 -f- a' x ' 1 -t- x".r " 3 -f- . . . 


a /i ■ 


dont la valeur moyenne sera aussi m*, quels que soient les 
coefficients a', a", etc. L’erreur moyenne à craindre lors- 
qu’on substitue la valeur m° à une valeur de y, déterminée - 
d’après les erreurs fortuites .r, x x", etc., sera, d’après 
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les principes précédents, donnée par la formule 

v /(„* - a!‘ („'< — ) -H a"> — /«"• ) + ...' 

i -+- a' fi 1 -t- a" fi" 2 -+-... 

Pour que celte erreur soit la plus petite possible, il faudra 
poser s - 


n' — in ' 


ni 
- in 


-, A 







Ces valeurs ne pourront pas s’évaluer tant qu’on ne eon- 
naitra .pas les rapports y, — etc. Dans l’ignorance où l’on 

est de leur valeur exacte (*), le plus sur sera de les sup- 
poser égaux entre eux ( noyez art. H), et l’on aura alors 



c’est-à-dire que les coefficients a v , etc., doivent être sup- 

posés égaux aux poids relatifs des diverses observations, eu 
prenant pour unité le poids de celle à' laquelle correspond 
l’erreur x. Ceci posé, désignons, comme ci-dessus, par <7 le 
nombre des erreurs proposées; la valeur moyenne de l’ex- 
pression 

X* -f- a' x' x -+- a." x "' 1 -f - . . . 

G 

sera égale à m*, et lorsque nous prendrons, pour la vraie 
valeur de «t*, la valeur de cette expression déterminée 
au moyen des erreurs fortuites x, x', x", etc., l’erreur 
moyenne à craindre sera 

^ n 1 -t- a' , n ' 1 -t- a" ? /»" 4 -+-... — a ni' 


(*) On ne conçoit la possibilité de déterminer exactement p . , ( u', /a", etc., 
que dans le seul cas où , pai* la nature de la fonction ç>, les erreurs x , x', 
x", etc., proportionnelles à /x, /x', /x", etc., seraient également probables, 
c’est-à-dire le cas où • . 

p (*) = #»y (/***) =>V (*'*) • • • • 

, ( Noie de M. Qavu.) 
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et, enfin, s'il est permis de supposer les quantités >>', 
u ", etc., proportionnelles à m , ni', ni ", etc., celte ex- 
pression se réduira à * 




résultat identique à celui que nous avons trouvé dans le cas 
•où les observations sont toutes de même espèce. 

. • ' 17 , 


Lorsqu’une observation dont la précision n’est pas abso- 
lue, fait connaître "une certaine quantité liée analytique- 
ment à une grandeur inconnue, le résultat de celte obser- 
vation peut fournir pour l’inconnue une valeur erronée, 
mais dans la détermination de laquelle il n’y a, rien d’ar- 
bitraire (jui puisse donner lieu à un choix plus ou moins 
vraisemblable. • . 

Mais si plusieurs fonctions de la même inconnue sont 
données par des observations imparfaites, chaque observa- * 
lion fournira une valeur de l’inconnue, et l’on pourra éga- 
lement obtenir des valeurs, par la combinaison de plusieurs 
observations. Il y a évidemment une infinité de manières 
d’y parvenir; le résultat sera soumis, dans tous les cas, à 
une erreur possible. Selon la combinaison adoptée, l’erreur 
moyenne à craindre pourra être plus ou moins grande. 

La même chose aura lieu si plusieurs quantités obser- 
vées dépendent à la fois de plusieurs inconnues. Selon .que 
le nombre des observations sera égal au nombre des incon- 
nues, ou plus petit ou plus grand que ce nombre, le pro- 
blème sera déterminé, ou indéterminé, ou plus que déter- 
miné (du moins en général), et, dans ce troisième cas, les 
observations pourront être combinées d’une infinité de 
manières pour fournir les valeurs dés inconnues. Parmi ces' 
combinaisons, il faudra clioisirles plus avantageuses, ç’esl- 
à-dirc celles qui fournissent des valeurs pour lesquelles 
l’erreur moyenne à craindre est Ja moindrfc possible.. Ce 
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problème ejt certainement le plus important tjue présente 
l’application des mathématiques à la philosophie naturelle. 

Nous avons montré, dans la Théorie du Mouvement des 
Corps célestes , comnient on trouve les valeurs les plus pro- 
bables des inconnues lorsque l’on connaît la loi de proba- 
bilité des erreurs des observations , et comme, dans presque 
tous les ras, celte loi par sa nature reste hypothétique, 
nous avons appliqué celte théorie à l’hypothèse très-plau- 
sible, que la probabilité de l’erreur x soit proportionnelle 
à e~ h ' x '\ de là cette méthode que j’ai suivie, surtout dans 
les calculs astronomiques, et que maintenant la plupart 
des calculateurs emploient sous le nom de Méthode des 
moindres carrés. 

Dans la suite-; Laplace, considérant la question sous un 
, autre point de vue, montra que ce principe est préférable 

• à tous les autres, quelle que soit la loi de probabilité des 

erreurs, pourvu que le nombre des observations soit très- 
grand. Mais. lorsque ce nombre est restreint, la question 
est encore intacte; de sorte que, si l’on rejette notre loi 
hypothétique , la méthode des moindres carrés serait préfé- 
rable aux autres, par la seule raison qu’elle conduit à des 
calculs plus simples. ' . 

Nous espérons donc être agréable aux géomètres eu dé- 
montrant, dans ce Mémoire, que la méthode des moindres 
carrés fournit les combinaisons les plus avantageuses des 
observations, non-seulement approximativement, mais en- 
core d une manière absolue, et cela quelle que soit la loi de 
probabilité des erreurs et quel que soit le nombre des ob- 
servations, pourvu que l’on adopte pour l’erreur moyenne, 
non pas la définition de Laplace, mais celle que nous avons 

donnée dans les art. 6 et 6. 

• - . 

Il est nécessaire d’avertir ici que, dans les recherches 
suivantes , il ne sera question que des erreurs fortuites dimi- 
nuées de leur partie constante. C’est à l’observateur qu’il 
appartient d’éloigner soigneusement les causes d’erreurs 
constantes/ . > • 


•J 
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Nous réservons pour une autre occasion l’examen du cas 
où les observations sont affectées d'une erreur constante 
inconnue, et nous traiterons cette question dans un autre 
Mt '■moire. 

18 . 

PROBLÈME. 

Soit U une fonction donnée des inconnues V, Y', Y", etc. 
on demande l’erreur moyenne M à craindre dans la dé- 
termination de la 'valeur de U , lorsque, au lieu des véri- 
tables valeurs de Y, Y', Y 7 ', etc., on prend les valeurs 
déduites d’observations indépendantes les unes des autres; 
m , m', tu",’ etc., étant les erreurs moyennes qui corres- 
pondent à ces diverses observations. 

Solution. — Désignons par e, e', e", etc., les erreurs des 
valeurs observées X , V', Y", etd.j l’erreur qui en résultera, 
pour la valeur de là fonction [J, pourra s’exprimer par la 
fonction linéaire 

A c — t- A’ e s" e” . . . — E , 

. . ,, ... , 7 ,, . , r/U r/U r/U 

ou A, A , A , etc. , représentent les dérivées —, —,■> -jÿ,f etc. , 

lorsqu’on y remplace Y, Y 7 , Y", etc., par leurs vraies valeurs. 

Cette valeur de E est évidente si l’on suppose les obser- 
vations assez exactes pour que les carrés cl les produits des 
erreurs soient négligeables. 11 résulte de là que la valeur 
moyenne de E est nulle, puisque l’on suppose que les er- 
reurs des observations n’ont plus de partie constante. Or 
l’erreur moyenne M, à craindre dans la valeur de U, sera 
la racine carrée de la valeur moyenne de E*,' c’est-à-dire 
que M 2 sera la valeur moyenne de la somme 

Ve'+ X'V’-t- ..+ aXX' cc'+ 2 II" et" + 2 YVVe"-+-. . . ; . 

mais la valeur moyenne de X’e* est A s wt ! , celle de À/’e' 1 
est V i m ,î , etc., enfin les valeurs moyennes des produits 
2 Xk'cc' sont toutes nullcs; donc on aura 

M = {JV/ir-f- X' v m'’ 4- X" J /«"'-t-. . 
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11 est bon {l’ajouter plusieurs remarques à cette solution. 

I. Puisqu’on néglige les puissances des erreurs qui sont 
supérieures à la première , nous pourrons , dans notre for- 
mule, prendre pour A, V, )/', etc., les valeurs des coefficients 

différentiels » etc., déduites des valeurs observées V, V, 
a y 

V", etc. Toutes les fois que U est une fonction linéaire, 
eette substitution est rigoureusement exacte. 

II. Si, au lieu des erreurs moyennes, on préfère intro- 
duire les poids des observations, supposons que f>ip\ 
p'\ etc., soient les poids respectifs, l’unité étant arbi- 
traire, et P le poids de la valeur de U ; on aura 


V V» X"’ 



III. Soit T une autre fonction de V, V', V /; , etc. 5 posons 

• dT , dT _ „ '• 

dv~ *’ d\'~ *•' dy"~ r '" ' 

L’erreur commise sur T, en adoptant pour V, V, V 7 , etc., 
les résultats fournis par l’observation , sera 

x<? + 7.V + .. = E', 

et l’erreur moyenne à craindre dans cette détermination sera 

* -, _____ ___ __ 

-+- y." 1 ni' 2 x ,/I /n" , 4- . . .. 

11 est évident que les erreurs E et E' ne seront pas indépen- 
dantes l’une de l’autre, et que la valeur moyenne du pro- 
duit EE' ne sera pas nulle comme la valeur moyenne de ce'; 
elle sera égale à 

• xXm 3 -+- •/.' W m' 1 4- x" À" m" 1 + . . . . 

IV. Le problème comprend le cas où les valeurs des quan-, 
tités V, V', \ ?// , etc., ne sont pas données immédiatement 
par l’observation , mais sont déduites de combinaisons quel- 
conques d’observations directes. Pour que cette extension 
soit légitime , il faut que les détermi nations de ces quanti tés 
soient indépendantes , c’est-à-dire qu’elles soient fournies 
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parties observations différentes. Si cette condition d’indé- 
pendance n’était pas remplie, la formule qui donne la valeur 
de M ne serait plus exacte. Si , par exemple , une même ob- 
servation était employée, à la fois, dans la détermination ' 
de V et dans celle de V, les erreurs e et e' ne seraient plus 
indépendantes, et la valeur moyenne du produit ce' ne se- 
rait plus nulle. Si l’on'connaît, dans ce cas , la relation qui 
lie V et V aux résultats des observations simples dont ils 
dérivent, on pourra calculer la valeur moyenne du pro- 
duit ee', comme il est indiqué dans la remarque 111, et dès 
lors corriger la formule qui donnera M. 

19 . 

Soient V, \ !/ , V', etc., des fonctions des inconnues x, 
y , z , etc.; soient ta le nombre de ces fonctions, p le nom- 
bre des inconnues; supposons que des observations aient . - 
donné, immédiatement ou médiatement, L, L', L"', etc., 
pour valeurs des fonctions V, V, V', etc., de manière ce- 
pendant que ces déterminations soient absolument indé- 
pendantes les unes des autres. Si p est plus grand que ta, la 
recherche des inconnues est un problème indéterminé. Si 
p — ta, chacune des inconnues .r, y, z, etc., peut être re- 
gardée comme calculée en fonction de V, V, V', etc.; de, 
sorte que les valeurs des premières peuvent être déduites 
des valeurs observées de ces dernières, et l’article précédent 
nous permettra de calculer la précision relative de ces di- 
verses déterminations. Si p est plus petit que ta, chaque 
inconnue x, y, z , etc., pourra être exprimée d une infinité 
de manières , en fonction de \ , \ ', Y 7 ', etc., et, en général, 
ces valeurs seront différentes; elles devraient coïncider si 
les observations étaient, contrairement à nos hypothèses, 
d’une exactitude rigoureuse. 11 est clair, d’ailleurs, que les 
diverses combinaisons fournissent des résultats dont la pré- 
cision sera, en général, différente. 

D’ailleurs si, dans le deuxième et le troisième cas, les 
quantités V, Y ", etc., sont telles que ta — p i d’entre 


• % 
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elles, ou davantage, puissent être regardées comme des fonc- 
tions des autres, le problème est plus que déterminé relative- 
ment à ces dernières fonctions et indéterminé relativement 
aux inconnues x, y, z, etc. ; et l’on ne pourrait même pas 
déterminer ces dernières inconnues, quand bien même les 
fonctions V, etc. , seraient exactement connues : mais 

nous excluons ce cas de nos recherches. 

Si V, V, \", etc., ne sont pas des fonctions linéaires des 
inconnues, on pourra toujours leur attribuer cette forme, 
en remplaçant les inconnues primitives par leur différence 
avec leurs valeurs approchées, que l’on suppose connues ; 
les erreurs moyennes à craindre dans les déterminations 
V = L, V' = L', V" = L", . . . 
étant désignées respectivement par m , m\ m", etc., et les 
poids de ces déterminations, par p,p', p", etc , de telle sorte 
que * 

pm 7 = p' m '* — p" m' h 

Nous supposerons connus les rapports des erreurs 
moyennes ainsi que les poids, dont l’un sera pris arbitrai- 
rement. Si nous posons enfin 

(V-L )\[p ~v,. [Y — L') \fp' = v',. . 

les choses se passeront ensuite comme si des observations 
immédiates , également précises et dont l’erreur moyenne 
aurait pour valeur m \/p, avaient donné 

i>= o, </=o, p" = o,.. . 

20 . 

PROBLÈME. 

Désignons par v, v', v", etc., les fonctions linéaires * 
suivantes des indéterminées x, y, z, etc., 

u = ax -+- by ■+• cz -f - . . . -+- /, 
v'z=a'x •+■ b’y-hc'z+.. 4- 1' , 
v" r=z a" x - i- b" y -f- c" z-y - . . ;-t- l", ■ 
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Parmi tous les systèmes des coefficients x, x', x", etc., 
qui donnent identiquement 

r. v -+- *V -+- x"p" -+-...= x — 

k étant indépendant de x, y, z, etc., trouver celui pour 
lequel x* -+- *" -+- y." 1 + . . . est minimum. 

Solution. — Posons 

ai’ -+- a' d + a" v" • . — I , 

bv +. b'v'4- b"v"-\ 

cv -+- cV c"v" -(-... = Ç , 


£ , > 7 , £ seront des fonctions linéaires de x,y, z, et l’on 
aura 

5=x^8 ! + r2^ + 2 ^ «c + . . . -+- 2 ni, 
ri = x'^ab 4- b 1 + z y. + , ■ • + y £/, 

Ç =: x V nr + ^ + 2 ^ c ’ + — 4- ^ * 

OÙ 

o 3 = a 7 -+- a' 7 4- o" 3 -4- . . , 

et de même pour les autres S. 

Le nombre des quantités £,»),£, etc., est égal au nom- 
bre er des inconnues x, y, z, etc.; on pourra donc obtenir, 
par élimination, une équation de la forme suivante (*), 

x — A 4- (a») l -4- •+• (ay) £ + 

qui sera satisfaite identiquement lorsqu’on remplacera 
j-,yj,£, par leurs valeurs (3). Par conséquent, si l’on 




(*) On verra plus loin la raison qui nous a conduit à désigner les coeffi- 
cients de cette formule par la notation («a), ( a/8), etc. 

(Note de M. Gauss.) 
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a (aa) 4- b (a(B) H- c (ay) 4- . . , 
a' (aa) 4- b' (a|î) 4- c' (ay) 4-, . 
a" (aa) 4- b" (a.\ J) 4- c" (ay) 4- . . 


on aura identiquement 


(5) 


a v 4- a! v 4- a" »" 4- . . . = x — A. 


Cette équation montre que parmi les différents systèmes 
de coefficients x , x', x", etc., on doit compter le système 


On aura d’ailleurs, pour un système quelconque, 

• (x — à) «> 4* (x' — a') 0'4-‘(*"- a")e"4-. . .= A — k; 
et cette équation, étant identique, entraîne les suivantes : 
(x — a)«4- (*' — a' ) a' 4- (x" — a" ) a" 4- . . . = o , 


(»-«){ + (x' — a') b' 4- (*" — a") b" 4- 
’ ( x _ a) c 4-(x'— a')c'4- (x"— a")c" 4- 


..= o, 

. . — o. 


Ajoutons ces équations après les avoir multipliées, respec- 
tivement, par (a*), (a(3 ),'(ay), etc., nous aurons, en vertu 
du système (4) , • 

(x _ «) a 4- (x' - a') a' 4- (x" — a") a" 4- . .= o, 
c’est-à-dire 

x ! 4- x' 1 4- . . . = a 1 4- a' 5 -4- . . . 4- ( x — a) 1 4- ( x' — a')’ 4- . . . j 
. par conséquent , la somme 

x» 4- 4- x" 1 4- . . . 

aura une valeur minimum, lorsque l’on aura 

^ / f H " 

x = a, x — a , x = a , , . . . 

D’ailleurs celte valeur minimum s’obtiendra de la manière 
suivante. • 
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L’équation (5) montre que l’on a 

nu 4 - n'a' 4- n"~a" 4-. . .= l , 
êa -h b'a'+b" a N 4-...= o, 
ca 4- c a' 4- c" a" 4- . ..= U, 


Multiplions ccs équations, respectivement, par (aa) , 
(«(3), (ay), etc., et ajoutons; en ayant égard aux rela- 
tions (4), on trouvera 


a’ a - 4- a 1 4- . , 


(aa). 


21 . 

Lorsque les observations auront dortné des équations 
approximatives • • 

• <> = 0, <>':=0, v" — O, •• .'. , 

il faudra, pour déterminer l’inconnue X, choisir une com- 
binaison de la forme suivante, 

7 r 

xv 4- xV 4- x" v" 4- . . . = o , 

telle que l’inconnue x acquière un coefficient égal à i, et 
que les autres inconnues se trouvent éliminées. 

Le poids de cette détermination sera , d’après l’art. 18, 


D’après l’article précédent, on obtiendra la détermina- 
tion la plus convenable, eu prenant 


alors x aura la valeur A.. On obtiendrait évidemment la 
même valeur sans connaître les multiplicateurs «, a', 
a", etc., en effectuant l’élimination sur les équations 
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le poids de cette détermination sera 

, i 

(*“1 ’ 

et l’erreur moyenne à craindre 

m \[p (aa) = m' (aa) = m" \[p" (aa) = . . . . 

Une marche analogue conduirait aux valeurs les plus 
convenables des autres inconnues y, z , etc., qui seront 
celles que l’on obtiendrait en effectuant l’élimination sur 
les équations 

5 = 0 , n = o, C = o,.... 

Si nous désignons par fl la somme 

ou, ce qui revient au môme, 

p{V — L )’-+-/>' ( V' — L')’ + p" (V" — L" , 
on aura évidemment 


- dSX 
2Ç = ^’ 


da da 

■-T~1 2 Ç = — , 

dy dz 


par conséquent, les valeurs des inconnues, déduites de la 
combinaison la plus convenable, et que nous pouvons 
appeler les valeurs les plus plausibles , sout précisément 
celles qui donnent à fl une valeur minimum. Or V — L 
représente la différence entre la valeur observée et la valeur 
calculée ; donc les valeurs les plus plausibles des inconnues 
sont celles qui rendent minimum la somme des carrés des 
différences entre les valeurs calculées et observées des 


quantités V, V, V ,; , etc., ces carrés étant respectivement 
multipliés par le poids des observations. J’avais établi 
depuis longtemps ce principe par d’autres considérations 
( Theoria Motus Corporum cœlestium). 

Si l’on veut assigner la précision relative de chacune des 
déterminations, il faut déduire des équations (3), les valeurs 
o. . , 3 
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de x, y, z , etc., qui se présenteront sous la forme suivante : 

I X — A -+- (aa) Ç ■+■ (ap) » -+- (ay) Ç -H. . 

y — B + (P«)5 + (PP)« -+■ (Pv) Ç 

« = C -f- (7a) 1 -+- (yp)» (77) £-+-• . •• 

Les valeurs les plus plausibles des inconnues x , y, z, etc., 
seront A , B , C , etc. Les poids de ces déterminations seront 

1 1 1 

(“*]’ (PP)’ (~77) ’ ' " ’ 

et les erreurs moyennes à craindre 

pour x , ni \jp [au) = m' \Jp' [ax ) > . . . , 

pour y, m 4 p [$$) — m ' 4F (PP) » • • • » 

pour z, m \fp (77) = m[ \[p' (77), • . • , 

ce qui s’accorde avec les résultats obtenus antérieurement 
( Theoria Motus Corporum cœlestium). 


22 . 


Le cas où il n’y a qu’une seule inconnue est le plus fré- 
quent et le plus simple de tous. On a alors 

V — x, V' = x, V" = *,...; 

il sera utile d’en dire quelques mots. 

On aura 

" = 4p, a’ = 4F> «" = 4F' >•••» 

1 ~ — l 4p > /'=- v4?, r = — l" 4F ' . 

et, par conséquent, 

5 = (P ■+■ P' + P" + • • • ) x — [p'L -H P f L' 4 - p n L" + . . . ) ; 
d’où 

[au] = 


A = 


P+P +p'... 
p~L + p'V + p‘ L" . 
P 4- p' + p" + ... 
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Ainsi , si , par plusieurs observations qui n’ont pas la 
même précision et dont les poids respectifs sont />,»', 
p", etc., ou a trouvé, pour une même quantité, une pre- 
mière valeur L, une deuxième L', une troisième I/, etc., 
la valeur la plus plausible sera 

pL 4-// I/ + />"L" -+•■■■ 

/<-+-/ -H ’ 

et le poids de cette détermination sera 

p -4- p' -H p ,f -H • » • • 

Si toutes les observations sont également plausibles , la 
valeur la plus probable sera 

L -t- L' -t- L" -t- . . . 

, 

GJ 

« 

c’est-à-dire la moyenne arithmétique entre les valeurs ob- 
servées; en prenant pour unité le poids d’une observation 
isolée, le poids de la moyenne sera ta. 


SECONDE PARTIE, 

I-HÉSENTtE LC 1 FÉVRIER l8l3, A LA SOCIÉTÉ ROYALE DE COTTINOCE. 


23 . 

Il reste encore à exposer quelques recherches destinées à 
étendre et à éclairer la théorie précédente. 

Cherchons d’abord si l’élimination qui fournit les varia- 
bles x , y, z , etc. , en fonction de £ , r? , £, etc. , est toujours 
possible. Puisque le nombre des équations est égal à celui 
des inconnues, on sait que cette élimination sera possible 
si £,»),£, etc. sont indépendants les uns des autres; dans 
le cas contraire, elle serait impossible. 

3 . 
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Supposons, pour un instant, que £, n , Ç , etc. , ne soient 
pas indépendantes , mais qu’il existe entre ces quantités l’é- 
quation identique 

o = Fi;4-G»4~H!54-»..4-K; 

nous en conclurons 

F 2 fll + & 2 ab + H ^ ac -t- . . . = o, 

F ^ 4- G 2 ^' + H Jic + ... = °, 

F2 ac + G 2 ic + H 2 C ’ + --=°> • 

F^aZ + G^ bl + cl 4-.. = — K. 

Posons 

! a F 4 - 4 G + c H 4-. . .= 0, 
a'F4-4'G4-c'H4-...= S', 
a" F 4- 4" G 4- c" H 4- . • . = 0", 


il viendra 

«84-a'fi' + a"0" + ...^o, 

40 4- 4'0' 4- b" 6" 4- . . .= o, 
c8 + c'0'+ c" 0" 4- . . .= o, 

Z0 4- Z' 0' 4- Z" 0" 4- ... = — K. 

En multipliant les équations (i), respectivement par 0, &, 
0", etc., et ajoutant , il vient 

0’4-0' 3 4-0"’4-...= o, 
et cette équation entraîne les suivantes : 

0 = o, 0' = o, 0" = o, . . . . 

De là nous concluons , en premier lieu, K = o. En second 
lieu, les équations (i) montrent que les fonctions 
t'", etc., sont telles, que leurs valeurs ne changent pas lors- 
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que les variables x, y, z, etc., prennent des accroissements 
proportionnels à F, G, H, etc. Il en sera évidemment de 
môme des fonctions V, etc. : or cela ne peut avoir 

lieu que dans le cas où il serait impossible de détermi- 
ner x,y, z, etc. , à l’aide des valeurs de V, \", V ", etc. , 
lors même que celles-ci seraient exactement connues; mais 
alors le problème serait indéterminé par sa nature , et nous 
exclurons ce cas de nos recherches. 


24. 

Désignons par (3, (3', j3", etc., des multiplicateurs qui 
jouent le même rôle relativement à l’inconnue^, que les 
multiplicateurs «, etc., relativement à l’inconnue x, 
c’est-à-dire tels, que l’on ait 

a (P * ) + b (pp) + c (p 7 )-t-...= p, 
a '(p»)+*'(pp) + c'(p 7 ) + ...= p', 

u"(P«) + à''(PP) + C "(p 7 )-K..= P", 


on aura identiquement 

P« -+- pV -f- p'V'. . . — y — B. 

Soient y, y , y ", etc. , les multiplicateurs analogues relatifs 
à la variable z tels, que l’on ait 

a ( 7 a) + é ( 7 p)-+-c ( 77 ) + ...= 7 , 
a' (yx) -+- b’(yp) + c'(yy) + ...= ÿ. 
a"( 7 a)-|-à"( 7 p)+c"( 77 )+.., = 7 ", 

et, par suite, 

y » -h y' e'-t- y" v" s — G. 


De la même manière que l’on a trouvé (art. 20) 



* J 

i 

■ i 

•1 
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nous trouverons ici 


^pa = o, ^pà = i, 

M 

c* 

il 

0 

., 2p /= - b > 

2v a = °. £n b ~ 0 ’ 

2t c = '»• 

• • » 2^7 1 = ~ c ’ 

et ainsi de suite. 


t 


On aura aussi , comme dans l’art. 20, 


^a> = (aa), 2P’ = (PP)> ^^(w) 

Multiplions les valeurs a, a", etc. (art. 20), respecti- 
vement, par (5, (i 1 , (i", etc., et ajoutons; nous aurons 

*P + « P' -t- a*p w ■+■ . . = ( ) , 

c’est-à-dire 

2>=(«p)- 

En multipliant /3, fi', (5", etc., respectivement, par a, a', 
a", etc., et ajoutant, on trouvera 

«P -h&'p' -h a"P"+ . . . = (P*); 

donc 

(ap) = (pa). 

On trouverait, de la même manière, 

(“?) =■ M = 2 “'/• (Pt) = (tP) = 2 Pt » — 

2o. 

Désignons par 1 , V, V, etc., les valeurs que prennent les 
fonctions y, *T, y", etc., quand on y remplace x,y, z, etc., 
par leurs valeurs les plus plausibles, A, B, G, etc., c’est-à- 
dire posons 

(i A -f- h B -f- c C -t - ... 4* / ~ \ , 
a' A-f- b'B -t-c'C + . .+ /' = V, 
a"A + J"B + c"C + ... + ;" = r, 

Si nous faisons 

x>-t- v* + r»+. . ,=m, 


A ■ 
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île telle sorte que M soit la valeur de la fonction fl, qui 
répond aux valeurs les plus plausibles des variables; M 
sera (art. 20) la valeur minimum de fl. 

Par suite, 

n 1 + «' V + fl" r + . . . 
sera la valeur que prend \ , lorsque 

x — A , /= B, z — C,.... 

Celte valeur est nulle, d’après la manière môme dont A , B, 
C, etc., ont été obtenus. On aura donc 

^aX = o; 

on obtiendrait de même 

2 e * =°’ 

et 

2*1=°» 2P X==0 ’ 2 vV=0 ’ 

Enfin, en multipliant les valeurs de 1 , X*, X' v , etc., respec- 
tivement, par X, X', X", et ajoutant, il viendra 

n + e v+ t'v-t — = v 4 - v> + x"’ ■+■ . . 

c’est-à-dire 

2 a = M - 

/ 

26 . 

Remplaçons , dans l’équation 

p = ax + by -+- es -4- . . . -t- /, 

x , j, r, par les expressions (7) [art. 21], on trouvera, en 
employant des réductions faciles, 

v = a l+p u -t- V ÇH-...- 4 -X, 

„'= a'Ç 4- p'ïi + y'Ç 4- l-X', 

r"=a"Ç+ p%+ 7''5 + - .. + X*, 
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Multipliant respectivement, ou ces équations ou les équa- 
tions (i) de l’art. 20, par A,X', 1", etc., et ajoutant en- 
suite, on obtient l’identité 

lp + XV-t- 1 V+...sM. 

27. 

La fonction îî peut se présenter sous plusieurs formes 
qu’il est important d’indiquer. 

Elevons au carré les équations (i) [art. 20], et ajoutons- 
les membre à membre, nous trouverons 

O. = ■*’ a’ -Pj' 3 i’ -P z 3 c 3 -p. . . -p ai/ ah 
-t- 2 xz 2 ac + 2 / z ^ be + . . . -p 

-paz^cZ-p. . .-p 

c’est la première forme. 

Multiplions les mêmes équations, respectivement, par e, 
té, té', etc. , et ajoutons, on aura 

a = Çx-f- ny -t- Çz + . .-P Iv -P /V -p {' v" -p. . . ; 

remplaçons k, té, té 7 , etc. , par les valeurs indiquées dans 
l’article précédent , nous trouverons 

tt ™ Ç x — p 7i y -p Çs - P . . . — AÇ — Bn — C Ç — . . . —P M, 

OU 

• U=Ç(* — A) + r>(y — B>+Ç(z — 0-P...-P M: 

c’est la seconde forme. 

Enfin , remplaçons , dans cette seconde forme , x — A , 
y — B, z — C, etc., par les expressions (y) [art. 21], nous 
obtenons la troisième forme : 

a = («a) V -P (pp) -+- ( ri ) ?+■ • •-+- 2 (*P) U 

-p 2 («v)?Ç + 2 (p7) , iÇ -P...-PM. 

On peut donner une quatrième forme qui résulte évi- 
demment de la troisième et des formules des articles précé- 
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dents, 

a = (e - \y + (/ - V)* + {*"— V')’ M , 

c’est-à-dire 

ü = M+2(" _ X)’. 

Sous cette dernière forme on voit clairement que M est la 
valeur minimum de 12. 

28 . 

Soient e, e', e", etc., les erreurs commises dans les ob- 
servations qui ont donné 

V = L, V' = L', V" = L" , . . . . 

Les vraies valeurs des fonctions V, \ T ', V", etc., seront alors 
L — e, V-e’, L" — e",. . . , 
et les vraies valeurs de v, f', v" , etc. , seront respectivement 
— —e’yjp', — e" 

par conséquent, la véritable valeur de x sera 

A — «e — a! e' \[p' — a." e" yfp" — . . 

et l’erreur commise dans la détermination la plus conve- 
nable de l’inconnue x sera , en la désignant par Ear , 

E ar = ae y fp a.' e' \j p -H a." e" \[f? •+- . . . . 

De même, l’erreur commise dans la détermination la plus 
convenable de la valeur dej - sera 

E / =5 p e + P'e' 4P' + P'V' >/?'+■ • • • 

La valeur moyenne du carré (Ex)’ sera 

m'p (a’ -+- a” -f- a"’ + ...)= m'p ( aa). 

La valeur moyenne de (Ej)’ sera de même 
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comme nous l’avons déjà reconnu. On peut également 
assigner la valeur moyenne du produit Ex. Ey, qui sera 

/?<’/? (a P + a'p' 4- a" P" 4-. . .) = m 7 p (ap). 

On énonce ces résultats plus brièvement de la manière sui- 
vante : 

Les valeurs moyennes des carrés (Ex)*, (Ey)*,etc., sont 
respectivement égales aux produits de par les quo- 
tients différentiels partiels du second ordre 
d'-a. d'ci 

TTç 7 ’ 7i ï»’ 

et la valeur moyenne d’un produit tel que Ex . Ey est le 
produit de ^ m'p par ei > regardant fl comme fonc- 

tion de £, ïî, £, etc. 

29 . 

Soit t une fonction donnée et linéaire des quantités x, 
y, z, etc., par exemple, 

t =fx 4- gy 4- hz 4 - . . . 4- h ; 

la valeur de t déduite des valeurs les plus plausibles de 
x, y, z, etc., sera 

fk 4- g'B 4- A C 4- h : 

nous la désignerons par K. En désignant par Et l’erreur 
commise en l’adoptant, on aura 

Et = / Ex + jEy 4- àEz 4-. . 

la valeur moyenne de cette erreur sera évidemment nulle, 
c’est-à-dire que l’erreur ne contiendra pas de partie con- 
stante, mais la valeur moyenne de (Et)*, c’est-à-dire de la 
somme 

/>(Ej)’ 4- 2 f g E.r.E_j - 4- 2 fhV.x.Y.z 4-, . . 

■+■ £’( e /)’ -f- 2 g/i Ey.Ez 4-. . . 

4- à’(Ez)’ 4-. . 
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sera, d’après l’article précédent, égale au produit de m*p 
par la somme 

/■(«*) + 2 /s(“P) + 2 / a (“7) -+-••• 

2 ^'(f7) + ••• 

■+• /j ’(77) +• • 

c’est-à-dire au produit de m* p par la valeur de la fonction 
il — M, lorsqu’on y fait 

?=/> » = £> K = /i, 

Désignons parta cette valeur de fl — M 5 l’erreur moyenne 
à craindre, lorsque l’on prend < = K, sera wjypw, et le 

poids de cette détermination sera -• 

1 W 

Puisque l’on a identiquement 
n — M = (x — A) Ç 4- (/ — B) » 4- (a — C)Ç 4*. . .» 
ta sera égal à la valeur de l’expression 

(x— B) g- 4 - (2 — C)/j-H... 

[qui représente (t — K)], dans laquelle on remplacera 
x, y, z, etc., par les valeurs correspondantes à £ =f, 
ri=g, Ç=h, etc. 

Enfin , observant que t , exprimé en fonction des quan- 
tités £, n, etc., aura K pour partie constante, si l’on sup- 
pose 

t — Fl; 4- G » 4- HÇ 4 - ... 4 - K., 

on aura 

SI±/F + ffG — t— /( II 4- . .. 

30 . 

Nous avons vu que la fonction fl acquiert son minimum 
absolu M, lorsque l’on y fait 

x — A , y — R , s 1 ,..., 

• * r ' k * 

ou, ce qui revient au même, 


Ç = o, w = o, Ç= o,. . . . 
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Si l'on attribue à l’une des inconnues une autre valeur ; que 
l’on fasse, par exemple, 

i = A + 4, 

les autres inconnues restant variables, fi pourra acquérir 
une valeur minimum relative , qui s’obtiendra à l’aide des 
équations 

rfft HCl 

X = A + A> — = o, — = 
et, par suite, 


or, puisque 


>1 = 0, Ç = o,. . . ; 


x= A 4- (««) ï 4- (<*P) >1 -h(a.y)Ç + . . ., 
on en conclut 

ç = — • 

(aa) 

On trouvera de même 


,=C + S2l»,. 

La valeur minimum relative de fi sera 


(aa)Ç’ + M = M+ —, :• 

Nous en conclurons, réciproquement, que si fi ne doit pas 
surpasser M -f- p.% la valeur de x est nécessairement com- 
prise entre les limites A — p\/[ôa), et A-f-py^(aa). H 
est important de remarquer que p yqao) devient égal à 
l’erreur moyenne à craindre dans la valeur la plus plau- 
sible de x, si l’on pose 

fl =: m \fp ; 

c’est-à-dire si p est l’erreur moyenne d’observations telles, 
que leur poids soit l’unité. 

Plus généralement, cherchons la plus petite valeur de la 
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fonction fi qui puisse correspondre à une valeur donnée 
de £, t désignant, comme dans l’article précédent, l’expres- 
sion linéaire 

f* + gy -H ht + . . . 4- *, 

dont la valeur la plus plausible est K ; désignons par K-t-x 
la valeur donnée de t. D’après la théorie des maximum et 
minimum, la solution du problème sera donnée par les 
équations 

d fl ^ dt 

dx dx ’ 


d fl dt 

dy ~ dy 
dQ. „ dt 


dz 


dz ’ 


ou 

5 = 9/, 

» = 9 g, 
Ç = 9 h, 


Q désignant un multiplicateur encore indéterminé. 

Si, comme dans l’article précédent, nous posons identi- 
quement, 

t=F5-t-Gn-l-H(;+... + K ) 

nous aurons 

K-t-*=0(/F + gG + /<Tl + ...) + K; 

d’où 



W 


ü) ayant la même signification que dans l’article précédent. 

Puisque fi — M est une fonction homogène du second 
degré, par rapport aux variables Ç, >7 , £, etc., sa valeur 
pour 

5 = 9/, »! = 9g, ç = e/i..., 

sera évidemment 


9’m , 
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et, par conséquent, la valeur minimum do fl, lorsque 
f = K + », 

sera 

y} 

M -f- 0 7 « = M-i 

W 

Réciproquement, si fl doit rester inférieure à une valeur 
donnée M-l-p.’, la valeur de Usera nécessairement com- 
prise entre les limites R — [X y'o), K-f-py^o), et p\/w sera 
l’erreur moyenne à craindre dans la valeur la plus plau- 
sible de f, si f* représente l’erreur moyenne d’observations 
dont le poids serait l’unité. 

31 . 

Lorsque le nombre des inconnues a:, y, z, etc., est un 
peu grand , là détermination des valeurs numériques de 
A , B, C, etc., par l’élimination ordinaire, est assez pénible. 
C’est pourquoi nous avons indiqué, dans la Théorie du 
Mouvement des Corps célestes , et développé plus tard, 
dans le Mémoire sur les éléments de l’orbite de Pallas 
( Commentaires de Gottingue, t. I), une méthode qui sim- 
plifie ce travail autant que possible. 

La fonction fl doit être ramenée à la forme suivante : 



où les diviseurs A» 0 , O", CD'", etc., sont des quantités 

déterminées; u ° , u" , etc., sont des fonctions linéaires 

de x,y, z, etc., telles que la seconde u' ne contient pas x, 
la troisième u " 11 e contient ni xni y, la quatrième ne con- 
tient ni x, ni y , ni z, et ainsi de suite, de sorte que la der- 
nière u(® — 0 ne contient que la dernière des inconnues 
x, y, z, etc.; enfin les coefficients de x, y, z, etc., dans 
«°, u", etc., sont respectivement égaux à *L°, us/, ©",etc. 

Alors on pose 

u° = o , u' = o, u" = o , u m — o, . . . , 
et l’on aura très-facilement les valeurs de x, y , z, etc., 
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en résolvant ces équations et commençant par la dernière. 
Je ne crois pas nécessaire de répéter de nouveau l'algorithme 
qui conduit à cette transformation de la fonction il (*). 

Mais l’élimination qu’il faut effectuer pour trouver les 
poids de ces déterminations exige des calculs bien plus 
longs encore. Nous avons montré, dans la Théorie du Mou- 
vement des Corps célestes, que le poids de la détermination 
de la dernière inconnue, qui entre seule dans id a ~ 
est égal au dernier ternie de la série des diviseurs X», 
il!/, e", etc. Cette recherche est facile; aussi plusieurs cal- 
culateurs, voulant éviter une élimination pénible, ont eu 
l’idée, faute d’autre méthode, de répéter la transformation 
indiquée en considérant successivement chaque inconnue 
comme la dernière. J’espère donc que les géomètres me 
sauront gré d’indiquer, pour calculer les poids des déter- 
minations, une méthode nouvelle qui ne me semble plus 
rien laisser à désirer sur ce point. 


Posons 


32 . 


x° x -+- tPo 0 ^ -+- e° z -+- 
Hl/ y -4- G' z -+• 
e" z + 

on aura identiquement 

^ dtl — \dx -+- n dy -f- Ç dz -+- . . . 

du • u' du’ u" du" 

f/0 _J_ - . —f- - - 

ju* ns.' e" 



+ C', 
4--C, 


/ , .ifM' G" , \ 

= “ [ dx + X- djr + x’ dz+ “) 

-t- u" [dz -f- . . . ) 


(*) On trouvera ces calculs ilnns une Note & la tin du volume. J B. 


Digitized by Google 


( 48 ) 


d’où nous déduirons 


(a) 


5 = «% 

Tft.' 

7i = — a° 4- « , 

a\>° 


g* 

i = — «• 


G' 

ifl/ 


Les valeurs de u°, u", etc., déduites de ces équations, 

se présenteront sous la forme suivante : 

( «’ = C, 

1 u' = A'|4- », 

«" = A" 5 4- B"« + Ç, 


M, 


(3) 


De la différentielle complète de l’équation 

fî = Ç(x — A) -4- n(y — B) 4- Ç(z — C)-4-. 

retranchons l’équation 

-dCï = Idx 4- rtdy 4- X,dz + . . . , 

2 

il viendra 

~ dn = ( x — K)dl 4 -{y — B) dn 4 - (z — C) Ç 

Cette expression doit être identique avec celle que l’on 
obtient à l’aide des équations (3), c’est-à-dire 

j r. rf ï + -^(A'‘*5 + *) + S»( A ' <, « + B'*+rfÇ)+.. i î 

«'A3 Ut) ^ 

on aura donc 


x — 


— H-A' 


( 4 ) 


y — 


Z = 


-A" 


G" 


—7 4- B — + 


G" 

a~ 

G 


-A, 
t- B, 
T, +• ■ *4- C, 


En substituant dans ce 3 expressions, les valeurs de u° , 
h', u", etc., tirées des équations (3), on aura effectué l’éli- 


Digitized by Google 


( 49 ) 

minaiioi). Pour déterminer les poidi 


(5) 


i A' 1 A"’ 

i(««)=Xî + S 7 + ër 

\ T B" 1 

(PW = ^T' + ô7' 


( 77 ) = 


1 

W' 


nous aurons 

AT* 

■ ■ ■ * > 
B’"' 



La simplicité de ces formules ne laisse rien à désirer. On 
trouverait des formules également simples pour exprimer 
les' autres coefficients (œj jÿ), (a/), ((2y), ete. ; mais, comme 
leur usage est moins fréquent, nous nous dispenserons de 
les exposer. 

. . 33. ’ 

L’importance du sujet nous a engagé à tout préparer 
pour le calcul et à former les expressions explicites des 
coefficients A', A", A!", B w , B w , etc. 

Ce calcul peut être abordé de deux manières : la pre- 
mière consiste à reporter, dans les équations (a), les valeurs 
de u 0 , u\ etc., déduites du système (3), qui doivent 
rendre ces équations identiques; et la seconde à exprimer, 
au contraire, que le système ( 2 ) devient identique lorsqu’on 
y substitue les valeurs de £, w, £, déduites du système (3). 

La première méthode conduit aux formules suivantes : 


A = °’ 
e° & 

_ + + = 

;ë + «? A+ ê» A +A =°' 


Ces formules feront connaître A 1 , A"> A w , etc. 

g. 4 
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On aura ensuite, 

O / 

îü? ~ 


B" = 1 




qui donneront B 7 , B w , etc. ; puis 

S+C*==o, 


qui feront connaître G", etc.; et ainsi de suite. 

La seconde méthode donne le système suivant : 

X® ,V 4- Di.* — 9 , 

d’où l’ou lire A'; 

X" A" 4 - Hb® B' 4 - e» = O , 

Db' B" 4- & =o, 

d on l’on tire B 1 et A 7 ; 

X® M" 4- ii!>° B'" 4- S® G" 4 - (D° = o , 
ub' B'" 4- & G" 4- = O , 

&>G" 4- Æ" = o, 

d’où l’on déduira G", B'", A'"; et ainsi de suite. 

Les deux systèmes de formules offrent des avantages à peu 
près égaux, lorsque l’on veut les poids des déterminations ' 
de toutes les inconnues x, y,, s, etc.; mais lorsqu’on ne 
cherche qu’une seule des quantités (**), ((3(3), (77), etc., • 
le premier système est bien préférable. 

D’ailleurs, la combinaison des équations (1) et (4) con- 
duit àux mêmes formules, et fournit, en outre, un 
second moyen d’obtenir les valeurs les plus plausibles 
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A, H, C , etc. , qui 

sont 




ê 

> 

fl 

1 

1 

- A'—, 

ï H 

— À" ^ — 

<m 

y ii x 



vit»' 

o" 

(ü" r •• 1 

• > 

B .=*_ ■ 


y" 

— B"~ — 

, put 

B" — 7 - — 



itt>' 

e" 


• y 

C = 


C 

O» 

yi " ■ 
cm x 

C ffi" * 

• • y 


L’autre calcul est identique avec le calcul ordinaire dans 
lequel on suppose 

u“ zx n , h' = o , u" — o , . 


34. 

Les résultats obtenus dans l’art. 32 ne sont que des cas 
particuliers d un théorème plus général qui peut s’énoncer 
de la manière suivante : 

Théorème — Si t représente la fonction linéaire sui- 
vante des inconnues x, y, z, etc., 

t =fx -t- g}- + hz-h . . . + k, 

dont V expression en Jonction des variables u°. u', u", etc.. 
-soit 

t -f- k r u! -f- h" u ’ K. , 

K sera la valeur la plus plausible de t , et le poids de cette 
détermination sera 

1 i 

A s, ^ub'X ,, 4^s ? 'A / ' ! +7T. ‘ 

Démonstration. — La première partie du théorème est 
évidente, puisque la valeur la plus plausible de t doit cor- 
respondre aux valeurs 

11“ = O, II' = O, li'=o 

Pour démontrer la seconde partie, remarquons que l’on a 

- r/fi = l;</x r>dy -f- Çrfz-f-, . . , 

• dt =fdx+ gdjr + /tdz + . . . ; 

• : ' 4 - 
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»*l par conséquent, lorsque 

?=/, ■» — g, l — h\..., 

on a 

(ICI — idt , 

quelles que soient les différentielles dx, dy, dz, etc. Il suit 
de là qu’en supposant toujours , 

% —fi r > — g i ? = h i • • • » 


on aura 


= À Vu® H- h' du' -h y du" -y .... 


Ur on voit facilement que si les différentielles dx , dy, 
dz , etc., sont indépendantes les unes des autres, il en sera 
(le même de du 0 , du', du", etc..-, nous aurons, par conséquent, 
pour 

% =/) r ‘—gl Ç 

u» = X® À®, u' = tli>'A', u" = S "A", 

Par conséquent, la valeur de fl correspondant aux mêmes 
hypothèses, sera 

X® A®’ -f- -A' A" 4- £"A"' 4- ... 4- M : 

ce qui , d’après l’art. 29, démontre l’exactitude de notre 
théorème. 

Si d’ailleurs on désire effectuer la transformation de la 
fonction t , sans avoir recours aux formules (4) (art. 32), on 
a immédiatement les relations 

/=x«A», 
g = ifi>®A®4-s)h , A', 

A = e° A® 4 - e'A'+ e"*", 


qui permettront de déterminer A®, A”', A", etc., et nous au- 
rons enfin 

K = — Ç A» — £*'— £'A" — . . . . 
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Nous traiterons particulièrement le problème suivant , 
taut à cause de sou utilité pratique, que de la simplicité de 
la solution : 

Trouver les changements que les valeurs les plus plau- 
sibles des inconnues subissent par l’adjonction d'une nou- 
velle équation, et assigner les poids de ces nouvelles 
déterminations. 

Conservons les notations précédentes. Les équations pri- 
mitives, réduites à avoir pour poids l’uni lé, seront 
i’ = o , d — o, ; 

on aura 

a = .< 1 + i>" 4- , 

-, K, etc., seront les dérivées partielles 
dSl .do dSl 

I o.dx‘ idy 2dz* ’ 

et enfin on aura, par l’élimination, 

x = A 4- («a) Ç 4- (a{5) » 4- (cry) £ 4- ■ • . , 

/ = B + (ap) Ç 4- (pp) >t 4 (p7) Ç + . • • > 

2 = C 4- (»y)Ç + {p7) r> 4- (77) Ç+. T . . , 

Supposons maintenant que l’on ait une nouvelle équa- 
tion approximative, 



dont nous supposerons le poids égal à l'unité. Cherchons 
les changements que subiront les valeurs les plus plausibles 
A , B, C, etc., et celles des coefficients (**)>. (Z 3 / 3 )» etc. 
Posons 

si -+- <**’ == Si*, 

1 rfü* _ 1 da* _ ; 1 da* _ . # 

ï dx ~ 5 ’ ’ 5 HT ~ ç ” ’ ” 

et soit 

le résultat de J élimination. 



( 3 4 ) . . 

Soit enfin, 

v*=fx + gy + hz- k, 
qui deviendra, en ayant égard aux équations (i), 
v* = F£ 4- G* + HS H- . . . 4- K , 

et posons 

Tf+Gg + HA 

' % • * 

K sera évidemment la valeur la plus plausible de la fonc- 
tion t'*, telle qu’elle résulte des équations primitives, sans 
avoir egard à la valeur o fournie par la nouvelle observa- 
tion, et - sera le poids de cette détermination. 

Or nous avons 

5* — ? »* = » ■+•' gv* , Ç* =.Ç 4- /iv*,- . . , 

et, par suite, 

F f * + G 7j* -h H + . . , + K = v* ( i + F/+ G # 4- Tl /< + , . . ) ; 

d’où l’on déduit 

- Fr+G /H- H E* -4- ... 4- K . 

I “f* b) 

, On a, en outre, 

* = A -+- (»«) (“?)>!*•+- (ay) Ç*4- ... 

— *’*[/(“*)] + (?(«#) 4- Æ(st?) 4-.. .] 

— AH-(aa)Ç*H-(ap)>i*H-...— Fc*, 

. =A + (aa)?*H- («p}„*4-... ( Fç* + G„*4- Il Ç* H-...R). 

I-f- w 

Nous déduirons de la*. 

A* = A--ÜL, 

I -f- w 

qui sera la valeur la plus plausible de x, déduite de toutes 
les observations. 
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nr conséquent. 


<«) — 


K J 


! W 

sera le poids de cette détermination. 

On trouvera de la meme manière, pour valeur la plus 
plausible de déduite de toutes les observations, , 

i -f- w *, 

le poids de celte détermination sera 


i 

W) 


G’ 


I -+- w 

et ainsi do suite. 

Le problème est donc résolu. 

\ joutons quelques remarques. 

1. F.n substituant les nouvelles valeurs A*, K*, etc., 

la fonction v* obtiendra la valeur la plus plausible 

*•*-». . * 

K K - (F/+G^+ 11// + ...)• = — K — , 

I +« I 61 

et, puisque l’on a, identiquement, 

F 0 * . H . 


I -f- W 


i H - w 


I H- o> 


K 

I -f- M 


le poids de celle détermination sera (art. 2G) 


1 "f* W 


~ ~ + l. 


F f -h G g — t— Il // -t- 

■•y " 

Ces résultats pourraient se déduire immédiatement des rè- 
gles exposées à la fin de l’art. 21. I, 'ensemble des équa- 
tions primitives avait, en effet, fourni la détermi nation 

t*=K, dont le poids était -* une observation nouvelle 
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donne une autre détermination v* == o, indépendante de la 
première, ayant pour poids l’unité; leur combinaison don- 
nera la détermination 

K > . ' -, 


qui aura pour poids 


1 -4- t*» 


~f~ I • 


II. On conclut de ce qui précède que, pour 
x = A*, y = B*, 3 — C , . , 

on devra avoir 

5* = o, **=o, C* = o, .... 

et, par suite, 

x _ /K gK- c = ~ àK 

1 -f- w l ï -+- <*» 

Puisque d’ailleurs 

ft = !• (x — A) -hv(y — B) H- Ç(z — C) M, 

' n* = n-t- p*% ■ 

on devra avoir, pour ces mêmes valeurs, 

« = ( F/ -h G g + H h + . . . ) + M = M + * 


et 


('+«)* 


» „ , K- 

< Si — M -f 7—; — ^ -H - — : — r, — M + 




i -i- w 


111. En comparant ces résultats avec ceux de l’ar- 
ticle 30, nous voyons ici que la fonction fl a la plus petite 
valeur qu’elle puisse obtenir lorsqu'on s’impose la condition 


K 
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Nous donnerons seulement ici la solution du problème 
suivant, qui a beaucoup d’analogie avec le précédent; mais 
nous nous abstiendrons d’indiquer la démonstration , à la- 
quelle le lecteur suppléera facilement en s’aidant de ce qui 
précède. 

Trouver les changements des valeurs les plus plausi- 
blés des inconnues et les poids des déterminations nou- 
velles, lorsque l'on change le poids de l'une des observa- 
tions primitives. 

Supposons qu’après avoir achevé le calcul on vienne à 
remarquer qu’on a attribué à une observation, à la pre- 
mière par exemple qui a donné V = L , un poids trop fort 
ou trop faible, et' qu'il serait plus exact de lui attribuer 
le poids p*, au lieu du poids p : il ne sera pas nécessaire 
de recommencer le calcul , mais il sera plus commode de 
former les corrections à l’aide des formules suivantes. 

Les valeurs les plus plausibles des inconnues deviendront 

(p*—p)«X _ > 


x == A — 


.r = K 


* = .c — 


■ ( P* — p) {a a -+- l> 0 -4 
{P* ~ P) PL 


•) 


(p* —r)( a * -+- *P + n -+-•••) 

(p' — p) 7* 

p) [ax + èp ■+• cy -H. - .) ’ 


■[P* 


les poids de ces déterminations s’obtiendront en divisant 
l’unité, respectivement , par 

(P*— /»)“’ 


(«0 
(PP) — 
irt) — 


p -+- (p* — p) («a ■+■ ép -f- cy '•+•. . .)’ 
(P'-P)P _____ 


p -h (p* — p) [ax + afi + cy-\-, , 

_ -(p* — p)r 


P •+• (P* ~ P) (« a 4- b ? + cy 4- . . . ) 

Cette solution convient au cas où, après avoir achevé le 
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calcul , il faudrait rejeter tout à fait l une des observations, 
puisque cela revient à faire p*=o; de même, p* = oc 
conviendra au cas où l’équation V =L, qui dans le cal- 
cul avait été regardée comme approchée, serait rigoureu- 
sement exacte. 

Si, après le calcul terminé, plusieurs équations nouvelles 
venaient s’adjoindre aux proposées, ou si les poids attribués 
à plusieurs d’entre elles étaient erronés, le calcul des cor- 
rections deviendrait trop compliqué, et il serait préférable 
de tout recommencer. ' 

' 37. 

iNous avons donné, dans les art. 15 et 16, une méthode 
pour déterminer, approximativement, la précision d’un 
système d’observations ('*'); mais cette méthode suppose 
connues exactement les erreurs réelles que l’on a effective- 
ment rencontrées dans une suite nombreuse d’observa- 
tions; or cette condition n’est remplie que bien rarement, 
pour ne pas dire jamais. 

Si les quantités dont l’observation fournit les valeurs 
approchées dépendent, suivant une loi donnée, d’une ou de 
plusieurs inconnues, on pourra trouver, par la méthode 
des moindres carrés, les valeurs les plus plausibles de ces 
inconnues; si, dàs lors, on calcule les valeurs correspon- 
dantes des grandeurs observées, ces dernières pourront être 
regardées comme différant peu des véritables : de sorte que 
leurs différences avec les valeurs observées, représenteront 
les erreurs commises avec une certitude d’autant plus 
grande, (pic les observations seront plus nombreuses. 
Telle est la marche suivie dans la pratique par les calcu- 

(*) Les recherches sur le même sujet insérée* par nous ( Zeitschrift fur 
Astronomie und verwanâte Wisscrischqften , vol. I , page i85) sont fondées 
sur l’hypothèse relative à la probabilité des erreurs à laquelle nous avions 
été conduit dans 4a Théorie du Mouvement des Corps célestes. 

{Mo te de M. Gauss.) 

Oi» trouvera pe Mémoire à la lin du volume. J. H. 
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latcurs, qui oui essayé, dans des eas compliqués, d évaluer 
à posteriori la précision des observations. Quoique suffi- 
sante dans bien des cas, cette méthode est, théoriquement, 
inexacte et pourrait quelquefois conduire à de graves er- 
reurs; c’est pourquoi il est très-important de traiter la 
question avec plus de soin. 

Conservons les notations de l’art. 19. La méthode dont 
il s’agit consiste à regarder A, B, C, etc., comme les véri- 
tables valeurs des inconnues x, y, z, etc., et X, X', X", etc., 
comme celles des fonctions e, v ' , v", etc. Si toutes les obser- 
vations ont une égale précision et que leur poids commun 


P = P ~ P =■ 


soit pris pour unité, ces mêmes quantités, changées désigne, 
représentent, dans celte supposition, les erreurs des obser- 
vations, et, par conséquent, d’après l’art. 15, 



v/ 


x , +v a ~+ V ,i -t-. . 

CI 



sera l’erreur moyenne des observations. Si les observations 
n'ont pas la même précision , — X , — X', — X", etc. , re- 
présenteront les erreurs des observations, respectivement 
multipliées par les racines carrées. des poids, et les règles 
de r art. 16 conduiront à la même formule. 



qui exprime déjà l’erreur moyenne de ces observations 
lorsque leur poids est égal à l imité. 

Mais le calcul exact exigerait évidemment que l’on rem- 
plaçât X, X', X", etc. , par les valeurs de e, t>', v", etc. , dé- 
duites des véritables valeurs des inconnues X,y, z , etc., 
et la quantité M par la valeur correspondante de il. 
Quoique l’on ne puisse pas assigner cette dernière valeur, 
nous sommes certain pourtant qu’elle est plus grande 
que M qui est son minimum : elle n’atteindrait cette li- 
mite que dans le cas, infiniment peu probable, où les va- 
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leurs véritables des inconnues se confondraient avec les plus 
plausibles. Nous pouvons donc affirmer, en général, que l'er- 
reur moyenne calculée par la pratique ordinaire est plus pe- 
tite que l’erreur moyenne exacte, et que l’on attribue, par 
conséquent, aux observations une trop grande précision. 

Voyons ce que donne une théorie rigoureuse. 

38. 

Avant tout, il faut chercher comment la quaulité M dé- 
pend des véritables erreurs des observations. Désignons ces 
erreurs, comme dans l'art. 28, pare, e! , e", etc., et po- 
sons, pour plus de simplicité, 

e\fp =zt, tf. \jp' — s’, e " \jp" = s” , . . , 


et 


Soient 


i sjp = m' V^7 = m" vV' = . . . = fi. 
A-x>, B-y, C — z% , 


les vraies valeurs des inconnues x, y , z, etc., pour les- 
quelles £ , k, £, etc. , soient, respectivement, — £°, ■ — n*.' 
— Ç°, etc. Les valeurs correspondantes de v, e', c", etc., se- 
ront évidemment 


B , — 6 , . . . 


de sorte qu’on aura 


t® . — fl t -f- fl — 1~ fl ■ s -+- 

r,’ =b E + b' «' H- //' t" H- 
V — c s -t- c'e' + r'V' ■+ 


•K® 5 (6 Ce £ — f— Ct £ — f— 

>•' = ^ E -t- <SV -+- 

= v* -+- 7V -t- 7" t" •+■ 


tfin 


il” = e : + e'’ 4-e'"+. 


sera la valeur de la fonction fl, correspondant aux vraies 
valeurs des variables x, y, z , etc. 
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Puisque l’on a identiquement 
fl.= M -+-(* — A) §4- (jr -B)„+(* — C)Ç+.. . , 
ou aura aussi 

M = — _>•*»• — z‘ £• — 

De là, il résulte évidemment que M est une fonction homo- 
gène du deuxième degré des erreurs e, e', e" , etc. ; cette fonc- 
tion pour diverses valeurs des erreurs pourra devenir plus 
ou moins grande. Dans l'ignorance où nous sommes de ces 
valeurs, il est bon d’examiner attentivement la fonction IM, 
et de calculer d’abord sa valeur moyenne d'après les prin- 
cipes du calcul des probabilités. Nous obtiendrons celte 
valeur moyenne en remplaçant les carrés e 1 , e", etc., par 
m m'*, etc., et en omettant les termesen ee', ee", etc., dont 
la valeur moyenne est zéro-, ou, ce qui revient au même, 
en remplaçant chaque carré s*, e", e"’, . , . , par p’, et en né- 
gligeant u', es",. . . . D’après cela, leterine fl° fournira et p’; 
le terme — x " £ 0 donnera 

— [a a a' a' + a" a" -)- . . . } p’ = — p’ ; 

chacune des autres parties donnera également — p . 7 , de sorte 
que la valeur moyenne totale sera 

[ér — p) p’, 

a désignant le nombre des observations, et p le nombre des 
inconnues. La vraie valeur de M pourra, suivant les cas 
que le hasard présentera, être plus grande ou plus petite 
que cette valeur moyenne, mais la différence sera d’autant 
moindre que le nombre des observations' sera plus grand; 
de sorte que . * ■ 

/" M 

• • Y m — p 

pourra être regardécomme une valeur approchée de p : par 
conséquent la valeur de p. , fournie par la méthode erronée 
dont nous avons parlé dans l’article précédent, devra être 
augmentée dans le rapport de — p à \ss. 



( )- 

39 . 

Afin de faire voir plus clairement jusqu’à quel point il est 
permis de regarder la valeur de M , fournie par les obser- 
vations, comme égale à la valeur exacte, il faut chercher 
quelle est l’erreur moyenne à craindre lorsque l’on fait 

M • • 


Celte erreur moyenne est la racine carrée de la valeur 
moyenne de la quantité 

' — .rT— jV — z'V— .. . — (o — p)j* 8 V 


y*if - 


que nous écrirons ainsi : 

/a •— — y* tf — z ■ 


CT p 






— [a*-**?— jrV — ... — (cr — — ft* ; 

et comme la valeur moyenne du second ternie est évidem- 
ment nulle, la question se réduit à chercher la valeur 
moyenne de la fonction 

¥ £= (n*— — y*rf r- z'X?—:. .y. 

Désignons cette valeur moyenne par N, l’erreur moyenne 
cherchée sera 


V; 


N 


(a—ç,y 


Si l’on développe la fonction Ÿ, on voit qu’elle est une 
fonction homogène des erreurs e, e\ e"-. etc., ou, ce qui re- 
vient au même, des quantités s, s', z", etc.; ou trouvera 
donc la valeur moyenne : 

i°. En remplaçant les quatrièmes puissances e\ e' 4 , 
e w , etc., par leurs valeurs moyennes; 
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■i". En remplaçant les produits e'c'*, e’ e"% etc-., par 
leurs valeurs moyennes, c’est-à-dire par in* in ”, m*m etc.; 
.3°. En négligeant les produits tels quee’e', c*e'e", etc. 
Nous supposerons (art. 1 (>) les valeurs moyennes de e‘, 
e'*, e"*, etc., proportionnelles à cn‘, m' 4 , m' n , etc, de sorte 

que les rapports des unes aux autres soient — i v* désignant 

la valeur moyenne des quatrièmes puissances des erreurs 
pour les observations dont le poids serait l’uni lé. 

Les règles précédentes pourraient se traduire de cette 
autre manière : 

Remplacer chaque quatrième puissance s‘, s H , s ,A , etc.j 
par v*; chaque produit e’s 1 , s 3 s ,/3 , etc., pat- u 4 , et négliger 
tous les ternies, tels que i‘ s i' ou î’s'e", - 

Ces principes étant compris, on verra facilement que : 

I. La valeur movenne de fl 03 est ' 

et»' 4-'(o’ — a) g*. * 

* * * . * * r 

II. La valeur moyenne du produit est 

a an* 4- ( a' a! -4- a" a." 4- . . .) g* = « a ( v* — g') 4- g'. 


Car 

a * 4- a' x 4- a" a" 4- . . . — i . 
De même, la valeur moyenne de e' ! .z° est 
flV(v‘ — g‘) 4- g‘ ; 
la valeur moyenne de e ,,! est 

a"x ( v* — g 1 ) 4- g‘; 


et ainsi de suite. 

Donc la valeur moyenne du produit 


c» -g J» 4- «" 3 4-, , .) x* ? ou a" JC*? 


sera 


1 — g 4- og , 


Les produits y 0 r° ou £l'' z !} £°, etc,, auront la même 
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valeur moyenne; donc le produit 

n*(x« •+• y , ïi t .+ z’K" -K . .) 

aura pour valeur moyenne 

pv 1 -+- p (tT l) Ji‘. 

III. Afin d’abréger les développements qui vont suivre, 
nous adopterons la notation suivante. Nous attribuerons à la 
caractéristique 2 un sens plus étendu que nous ne l’avons 
fait jusqu’ici, en lui faisant désigner la somme des termes 
semblables, mais non identiques, qui proviennent de toutes 
les permutations des observations. Nous aurons, d’après 
cette notation , • 

,r« =2«, 

ar 01 a 1 «’ 2 ^aa'ss'. 


Calculant par parties la valeur moyenne du terme .r 0 ’ £ oï , 
nous aurons d’abord, pour valeur moyenne du produit 

avrv . . 

a’aV + a 1 (fl' 1 -+- fl" 1 + . . . ) a* 

• = a 1 a’ (v‘ — pé ) -+- 

4 ’* * , _ ( ' . « 

De même, la valeur moyenne du produit a'* e'* £ 0 ’, est 

/t' i «'’fv* — pt‘) H- a' 1 (*‘2 fl ’ 1 » 

et ainsi de suite. 

Par conséquent, la valeur moyenne du produit 
. sera . ■ ; . , 

( v* — p* ) 2 : *’ ■+• f 1 ’ 2 

Or la valeur moyenne de est 

aaa'oa'pt*. 
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La valeur moyenne de *a"se"£°* est 
?. ax au fit 1 ; 

et ainsi de suite. D’où Ton conclut facilement que la valeur 
moyenne du produit 




est 
2 f*’ 


, 2 a<, “ v=ft '[(2' Ta ) ~ S" 1 *’] =fx< (•— 

Ceci posé, nous aurons pour valeur moyenne du pro- 
duit x 01 g 0 ’, 

(v* — 3 fi’ ) at’ 2 -+- pt* 

IV. On trouvera d’une manière analogue, pour valeur 
moyenne du produit x°j° £ 0 r; 0 , 

v‘ ajî + i' j3' 4- (*'2 a b x> P’ + (*’ ^ a fi b' a'. 

Or on a 

^flaè'p'=2' 7a \j£ i b 

y^aba' p' =2 « b -2“ i 3 ~ 2 " * “ ? • 

2"“ — ’’ 2 fl ?== °» ^bx-o-, 

cette valeur moyenne sera donc 

(»' — 3 f*’)2 fl4a ? + (*' ^'+2 aA ’2 a ^)‘ 

V. On trouverait, par un calcul semblable, que la valeur 
moyenne de x°z 0 £ 0 £° est 

(v> — 3p 4 ) 2«c«y + .f*'( ,+ 2# 


c. 
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et ainsi de suite. En additionnant, on obtient la valeur 
moyenne du produit 

* , $ 4 (x*ï* + .rV + 

cette valeur est 

(■>' — 3fi 4 )^ [o“( aa H- bp cy +. . .)] 

-+- (p-Hi)(* , + (*' •2 a P 4 "2 fltf , 2“ 7+ "') 

= (*' — 3p')2 («“ + H- cy +-•••)] + (p + 2 ) f 4 *- 

VI. On trouverait de la même manière 
\ (v‘ — 3 f* 4 ) 2 [6 p («a 4- bp -t- cy -+- . . .)] 4- (p 4- 2 ) p‘ 

pour valeur moyenne du produit 

7 V (x 4 ’?* 4-/V 4- zT4-.. •), 

et 

(v* — 3fi 4 ) y [cy (a et 4- bp 4- cy 4-. ■ .)] 4- (p 4- 2 ) ft 4 
pour valeur moyenne du produit 

4-.yV 4- *•?•+...); 

et ainsi de suite. 

Nous aurons donc, par l’addition, la valeur moyenne du 
carré 

(*•?• -t-yV 4- z»Ç» -H. 

elle sera , 

(v 4 — 3 fi 4 } 2 [(<*“ 4- bp 4- cy -t- . • . )* ] 4- (p* 4- 2 p)fi 4 . 

VU. Nous concluons enfin de tous ces préliminaires, 

N = (ej — 2 p) y 4 4- («r* — o — 2 np 4- 4 p 4- p’) ft 4 
4_ (./ 3fi 4 )2[( oa + bp 4- cy 4-. . ) J ] 

» ’ . . = (n — p) (y‘ — p 4 ) 4- (n — p)’ fi 4 

— (y 4 — 3 fi 4 ) j^p— 2 (a* 4- bp + cy 4-. J- 
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Donc l’erreur moyenne à craindre, lorsqu'on prendra 



i(>. 


La quantité 

2 [_(<*<* H- b$ + cy +. . .)*], 

qui entre dans l’expression précédente, ne peut pas géné- 
ralement se réduire à une forme plus simple. Cependant 
on peut assigner deux limites entre lesquelles sa valeur doit 
nécessairement être comprise : 

i°. On déduit facilement des relations précédentes, 

(fla 4- b p -t- cy — . . .)* 4- (fl a' 4- !>$' - f- cy' -f-., .)’ 

4- (fla" 4- b 8" ■+■ cy" -f-, . .y ... . = fla -t- bp 4- cy 4- . . . ; 

d’où nous concluons que 

aaq-èpq-cyq-... 

est une quantité positive plus petite que l’unité, ou du 
moins qu’elle n’est pas plus grande. Il en sera de même de 
la quantité .♦ ’ 

do! q- b'p + c'y' 
qui est égale à la somme 

(a'a q- ê' (5 4- c'y q- . . . )* q- ( a' a' q- b’ P' -f c’y' -i-, . .)’ 
q- (fl' a" q- i'p»q. e yq-...)>+.‘; ; 

de mente 

q!' a" q- b ' |ï ' q- c" y"q- . • • 

5. 
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sera plus petit que l’unité ; et ainsi de suite. Donc 

+ ip 4- C7 H )•] 

est nécessairement plus petit que ra. 
a°. On a 

^ («« + ip -t-cy = p, 


car 


= 2 c7=i,- " : 

d'où l’on déduit facilement que 

^[(aa + +c 7 +•. . .)’] 


• • C* 

est plus grand, ou du moins n’est pas plus petit, que -• 
Par conséquent le terme 

v* — 3 [*■ 

( 5J — p) 

est nécessairement compris entre les limites 


y, [p — 2 +-c 7 -+-. . J 


»* — V et v< — V t, 

TJ p CT “* p TJ 


ou bien, entre les limites plus étendues, 


y» — 3p' 

st — p 


et 


— P 


Donc le carré de l’erreur moyenne à craindre pour la va- 
leur 

' M 

f* = - — r 


est compris entre les limites 

2 v* — 4 fl 1 

SJ — p 


et 


2f l' 
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de sorte qu’on pourra atteindre un degré de précision aussi 
grand que l’on voudra , pourvu que le nombre des observa- 
tions soit suffisamment grand. 

Il est très-remarquable que dans l’hypothèse de l’art. 9 
(III) , sur laquelle nous nous étions autrefois appuyé pour 
établir la théorie des moindres carrés, le second terme 
du carré de l’erreur moyenne disparaît complètement (car 
on a v» — 3p* = o); et comme, pour trouver la valeur 
approchée p, de l’erreur moyenne des observations, il faut, 
dans tous les cas, traiter la somme 

V -I- X" + X" ! -I- . . . = M, 

comme si elle était égale à la somme des carrés des a — p 
erreurs fortuites, il en- résulte que, dans cette hypothèse, 
la précision de cette détermination devient égale à celle que 
nous avons trouvée, art,, lo, pour la détermination déduite 
de a — p erreurs vraies. 
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SUPPLÉMENT 

• . A LA 

THÉORIE DE LA COMBINAISON DES OBSERVATIONS 

QUI EXPOSE AUX MOINDRES ERREURS. 

PRÉSENTÉ LE l6 SEPTEMBRE 1 8 >G A LA SOCIETE ROTALE DE COTTI NOCE. *. 

I. 

Dans le Mémoire précédent, on a supposé que les quan- 
tités à déterminer, à l’aide d’observations imparfaites, 
fussent dépendantes de certains éléments inconnus, en fonc- 
tion desquels on sût les exprimer : le problème consistait 
alors à déduire des observations, le plus exactement pos- 
sible, la valeur de ces éléments. 

Dans la plupart des cas, c’est effectivement ainsi que la 
question se présenté; mais quelquefois il en est un peu 
autrement; et l’on pourrait même se demander, au premier 
abord, si le problème se ramène au précédent. Il n’est pas 
rare,. en effet, que lés quantités auxquelles se rapportent 
les observations ne soient pas explicitement exprimées en 
fonction de certains éléments, et qu’elles ne paraissent ré- 
ductibles à une telle forme que par des opérations difficiles 
ou ambiguës. Il arrive souvent, d’autre part, que la nature 
du problème fournit certaines conditions auxquelles les 
valeurs observées doivent rigoureusement satisfaire. 

Cependant, en y regardant de plus près, on aperçoit que 
ce cas ne diffère pas essentiellement du précédent, et qu’il 
peut s’y ramener. Si l’on nomme, en effet, ta le nombre des 
grandeurs observées, et n celui des équations de condition , 
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en choisissant, parmi les premières, ct — a quantités, rien 
ne nous empêche de les considérer comme nos seules in- 
connues, les autres étant regardées comme des fonctions de 
celles-là, que les équations de condition définissent. Par cet 
artifice, nous rentrerons dans le cas du Mémoire précé- 
dent. 

Néanmoins, quoique cette marche conduise souvent au 
résultat d’une manière assez commode, on ne peut nier 
qu’elle ne soit moins naturelle, et il est, par cotisé-* 
quent, désirable de traiter le problème sous l’autre forme, 
qui admet, du reste, une solution très- élégante. Il y a 
plus : cette solution nouvelle, conduisant à des calculs plus 
rapides que la précédente toutes les fois que a est moindre 

que ^ ct , ou, ce qui revient au même, toutes les fois que le 

nombre des éléments que nous avons désigné par p dans le 


Mémoire précédent, est plus grand que il faudra, dans 


ce cas, préférer la solution nouvelle, quand bien même il 
serait facile, par la nature du problème, de, faire dispa- 
raître sans ambiguïté les équations de condition. 


2 . 

Désignons par e, y', v" , etc., les quantités, en nombre n, 
dont les valeurs sont fournies par l’observation. Supposons 
qu’une iiiconnue dépende de celles-là et soit exprimée par 
une fonction connue u,de y, y', v", etc. Soient /, l", etc., 

ce que deviennent les quotients différentiels 

du du du 

dv ’ dv' ’ du " ’ 

lorsqu’on y substitue à v , y', y ", etc., leurs valeurs véri- 
tables. Si l’on substituait à y, y', y", etc., dans la fonc- 
tion n, leurs valeurs véritables, on obtiendrait aussi la 
véritable valeur de «; mais si les observations sont alleo- 
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tées d’erreurs e, e', e", ele., il en résultera pour u une 
erreur totale représentée par 

le te' + 

pourvu que l’on puisse, comme nous le supposerons tou- 
jours (lorsque u n’est pas linéaire), négliger les carrés et 
les produits des erreurs e, e', e", etc. 

Quoique la grandeur des erreurs e, e', e", etc., soit in- 
certaine, l’incertitude attachée à la valeur trouvée pour u 
peut généralement se mesurer par l 'erreur moyenne à 
craindre dans la détermination adoptée. D’après les prin- 
cipes développés dans le premier Mémoire, cette erreur 
moyenne est 

V^/H* -f- rv/f' 3 -h~l" 1 m"‘ -+-777 , 

m, m', m", etc., étant les erreurs moyennes des diverses 
observations. Si toutes les observations sont affectées du 
même degré d’incertitude, cette expression devient 

m \/r- -t- r -+- . . 

Il est clair d’ailleurs qu’au degré d’approximation auquel 
nous nous arrêtons, on peut remplacer /, l\ etc., par 
les valeurs que prennent les coefficients différentiels 

dit tlu du 
dv ’ dv' ’ dv" ’ 

lorsqu’on y remplace e, v' , v" , etc., par leurs valeurs ob- 
servées. 

, 1 . 

/ ’ . . , 

Lorsque les quantités n, e', u", etc., sont indépendantes, 
l’inconnue ne peut se déterminer que d’une seule manière, 
et l’incertitude attachée au résultat ne peut être ni évitée 
ni diminuée. Les observations fournissent une valeur de 
l’inconnue qui n’a , dans ce cas, rien d’arbitraire. Il en est 
tout autrement lorsque les quantités ^, e', e", etc., sont 



liées par des relations nécessaires, que nous supposerons 
exprimées par a équations, 


X = o , Y — o , Z =r o , . . , , 

X, Y, Z, etc., désignant des fonctions données des indétermi- 
nées v , v ' , v " , etc.; car, à la fonction n, on peut, dans ce cas, 
substituer toute autre fonction U telle, que la différence 
U — u s’évanouisse identiquement en vertu des équations 

X = o, Y = o, Z=o,.... 


Si les observations étaient rigoureusement exactes, cette 
substitution ne changerait en rien le résultat; mais, en rai- 
son des erreurs inévitables, à chaque forme adoptée pour u 
correspondra un résultat différent, et f erreur commise, au 
lieu d’ôtre 

le + /V + K . . , 


deviendra 


en désignant par L, L', L", etc., les quotients différen- 
tiels 

rl U (/U d\j 
tlv, ’ dv' ’ dv" * 

Quoiqu'il soit impossible d’assigner la valeur des diverses 
erreurs , nous pouvons cependant comparer les erreurs 
moyennes à craindre dans les diverses combinaisons. La 
combinaison la plus avantageuse sera celle qui donnera, à 
l’erreur moyenne, la valeur minimum. Cette erreur étant 

^L’iw’-t- . . , 

«■ • X 

nous devons chercher à rendre la somme 


Wm‘ -+• V'm“ H- + . • . , 

aussi petite que possible. 
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4 . 


Les fondions U, en nombre infini, par lesquelles on peut 
remplacer «,‘ne différeront les unes des autres, dans nos 
recherches, que par les valeurs qu’ elles fourniront pour 
L,L', L", etc.: nous devons donc, avant tout, chercher 
les relations qui existent entre les systèmes de valeurs que 
peuvent prendre ces coefficients. Désignons par 




les valeurs que prennent les coefficients 


dX 

dX 

dX 


dv' 

dv '* 

dv" ' ' 

» 

d Y 

d Y 

dY 


‘ dv ’ 

dv’' 

dv" ’ ' 

* * » 

dZ 

dZ 

dZ 


de’ 

dd ’ 

dv" ’ - • 

» 1 


si l’on y substitue pour e, v', v>" , etc., leurs Valeurs véri- 
tables. Il est clair que si Ton donne à e, v' , v" , etc., des 
accroissements di 1 , dv 1 , dv" , etc., qui ne changent pas 
X, Y, Z, etc., et leur laissent, par conséquent, la valeur 
zéro, ces accroissements, qui satisferont aux équations 

o = adv + n'dv' -f a" d v" 4- . . . , 
o = bdv -f- b'dv’ -+- b" dv" 

O = cdv c'dv' c"dv" - . . , 


ne changeront rien à la valeur de U — h, et I on aura , par 
conséquent, 

o = (/ — L) dv + (/' — L') dv' -p (/" — L" j dv" + 

Un en conclut facilement que L , I/, L", etc., doivent avoir 
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L = l 4- ax -y b y + cz + , 

L' = ï 4- a x 4- b' j- 4- c'a 4 - . . , 

L"== /"+ a" x- 1 - b" y -+- c "z 4 - . . . , 

x, j, z, etc., désignant des multiplicateursdétcrminés. Réci- 
proquement, il est clair que, pour toutes les valeurs de .r, 

y, z , etc., on pourra former une fonction U, pour laquelle 
les valeurs L, L', L", etc., seront précisément celles que 
fournissent ces équations , et cette fonction pourra , d’après 
ce qui précède, être substituée à «. La forme la plus simple 
qu’on puisse lui donner est 

U — n 4- x X 4- y Y 4~ z Z 4“ ... 4- w', 

u' désignant une fonction de e, v", etc., qui s’annule 
identiquement lorsque X , Y, Z, etc., sont nuis, et dont la 
valeur, dans le cas actuel, sera maximum ou minimum (car 
scs dérivées devront toutes s’annuler). Mais l’introduction 
de cette fonction n’apporte aucune différence dans les ré- 
sultats. 

5 . 

Il est maintenant facile d’assigner à x,y, z, etc., des 
• valeurs telles, que la somme 

L- m 1 4- L ' 1 m' 1 4 - L" m"’ + . . . . . . 

soit un'miniraum. •*» -.* 

Il est clair que pour atteindre ce but, la connaissance des 
erreurs moyennes absolues m, m', m ", etc., n’est pas néces- 
saire; il suffit de connaître leurs rapports. Introduisons, en 
effet, au lieu de ces quantités, les poids des observations . p, 
//, p , etc., c’est-à-dire des nombres réciproquement pro- 
portionnels à «i*, /7i ,s , m"*, etc. Les quantités x, jr t s , etc. , 
devront être déterminées de telle sorte que le polynôme 
(ax 4- by 4- rz 4- . . . 4 - 1 )’ i (a'x+ b'y+ c'a. ..4- /')’ 

P P 

acquière une valeur minimum. Supposons que x°, y", 
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x°, etc., soient les valeurs déterminées de x, y, z , etc., 
auxquelles répond ce minimum, et adoptons les notations 
suivantes : 


à 1 a' 7 a" 1 

~~ -4- • — 7 - 4~ ■ ■ ,, + . 
P P P 

a b a' b’ a" b" 

..— (aa), 
..= (ab), 
• • = (««)» 

■ 

f • // ' * 

P P P 

ne a'c' a"c" 

1 ~ H “TT "4“ . 

P P P 


b ’ b' 7 b " 7 



— -4- — •+- — 4- • 
P P P 

; 


c 7 c' 7 c" 7 



- H ; 4- — + • 

♦ - * 

P P P 


al a'V a"l" 



~ 4 ; — | - h. 


P P P 

‘ ‘ 

b t b'V b"l" 


, - 

1 H — + ■ 

P P P 

■■ — (bl), 


cl c'I' c'T 



~ H y -f- — — -+- . 

P P P 


“ 


La condition de minimum exige évidemment que l’on 
ait 

o =(aa)x‘ -4- [ab)y 7 -4- (rtc) z # -f-, . . -t- {al), 
o = ( ba) x" -+j [bb) y 0 ■+ ( bc) z* -4- . . . -f- ( bt), ' 
o = (ca) x° -4- (cb)y° - 4 - (rc) z" - 4 - ... - 4 - (cl), 

Après que ces équations auront fourni les valeurs de x 0 ,^ 0 , 
z°, etc., on posera 

a x" - 4 - b y° - 4 - r z" - 4 - . . . - 4 - / = L, 
a' x° -h b' y" -h c’ z" + . . . + l’ = V , 
a" x' - 4 - b" y* - 4 - c"z* 4 - . . . - 4 - Y— L", 

et la fonction la plus propre à déterminer notre inconnue, 
à laquelle correspond l’erreur moyenne minimum, sera 
celle dont les coefficients différentiels, pour les valeurs con- 
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sidérées des variables, seront égaux «à L, L', I.", etc. Le poids 
de la détermination ainsi obtenue, sera 



c’est-à-dire que ^ sera précisément la valeur que prend le 

polynôme considéré plus haut, pour les valeurs des variables 
x,y, z , etc., qui satisfont aux équations (i). 


6 . 

i 

Dans l’article précédent, nous avons montré à déterminer 
la fonction U qui fournit la détermination la plus conve- 
nable de l’inconnue u. Examinons actuellement quelle 
est la valeur q ui en résulte Désignons-la par K: on l’ob- 
tiendra en substituant, dans U, les valeurs observées des 
quantités u, m', u", etc. Soit A la valeur que prend « lors- 
qu’on y fait les mômes substitutions, soit enfin y. la véri- 
table valeur de cette inconnue, telle qu’on l’obtiendrait 
par la substitution des valeurs véritables de e, </, v " , etc., 
soit dans u , soit dans U. On aura 

h = * -t- le -+- té l" e" . . . , 

R = * + Le+ L V L" e” . . .. , 

et , par conséquent , 


K = A + (L — l)e + (L' — /') e'-h (L" — 

Substituant, dans cette équation, à la place de L — l, 
L' — L" — etc., leurs valeurs fournies par (a), et 
posant 


( 4 ) 


ae +■ a' e' -|- a" e" -f-. .. = A, 

be + b' e' + b" e" 

ce c' e' -f- c"e" -f-. . ,= S, 


nous aur&ns 



* » 


(78) 

Il n’est pas possible de calculer JL, us,, B, etc., par le moyen 
des formules ( 4 ), caries erreurs e,e',e", etc., qui y figurent, •’ 
ont des valeurs inconnues , mais on voit facilement que ces 
quantités «L, tP.,, S, etc., ne sont autre chose que les valeurs 
de X, Y, Z, etc., qui correspondent aux valeurs observées 
de v, e', e", etc., et alors le système des équations (i), (3), 

(5), forme la solution complète de notre problème. Il est 
clair, en effet, que l’on peut appliquer au calcul de a, a', 
a", . .. , b, b', b", etc., la remarque faite à la fin de l’ar- 
ticle 2, à l’occasion des quantités /, etc., c’cst-â- 

dire remplacer les valeurs véritables de e, e', e", etc,, 
par les valeurs observées. 

7. 

On peut substituer à la formule (3), qui représente le 
poids de la détermination la plus probable, plusieurs ex- 
pressions qu’il est utile d’indiquer; remarquons, d’abord, 
qu’en ajoutant les équations ( 2 ) après les avoir multipliées 
a a' a" 

par-, —^1 etc., on aura 

P P P , ' . 

. , , , . «L a' L' a" L" • 

(■ aa ) x"-t- 1- ( ac z“ 1 — H - K... 

' . p p p 

Le premier membre est nul ; en désignant donc, d’après 
la notation adoptée, le second membre par (aL), on aura 

(«L) = o, 

et de même 

(4L) = o, (cL)=o,.... . 


L 1/ L" 

Multiplions ensuite les équations ( 2 ) par -, —,■> etc.; 

nous aurons, en les ajoutant, 

/L l' V l" L" L’ L' 1 L"’ 

h— - H — + — + — 

P P P 


P P P 

et , par suite , nous avons cette seconde expression du poids, 
P = 1 


IL F L' l"L" 

1 1 TT ' 

P P P 
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Si, enfin, nous ajoutons les mémos équations (a) après 
l V l" 

les avoir multipliées par - , — , — -, etc. , nous obtiendrons 
P P P 

une troisième expression du poids, 

P [ 

(a/)x" 4 - ( bl ) y* 4 - (c/)z°- h .. . 4- (//)’ 
où l’on a posé, conformément à la notation précédente, 



On passera facilement de là à la quatrième expression 
du poids , 


i ■=(//) — (aa) x •’ — ( bb)y 01 — (rc) z” — ... 

— a (ab) x’ y — 2(ac)x°z' 1 — .2 (bc)y°zf — . . . . 


8 . 

La solution générale que nous venons d’exposer s’ap- 
plique principalement au cas où l’on n’a qu’une seule in- 
connue à déterminer. Lorsqu’on cherche, au contraire,- les 
valeurs les plus plausibles de plusieurs inconnues, dépen- 
dant des mêmes observations, ou lorsqu’on ignore quelles 
sont les inconnues qu’il faut, de préférence, déduire des- 
observations , il convient de procéder d’une manière diffe- 
rente , dont nous allons actuellement nous occuper. 

Considérons x, y, z , etc., comme des indéterminées, et 
posons 

(aa)x 4 - ( ab) y 4- (ac) z 4 - . . . = £ , 

(ab ) x -I- ( bb ) y 4- { bc jz -t- . . = n , 

(ac) x 4 -(bc) y H- (cc)z 4-. . .= £, 

Supposons que l’on en déduise, par l’élimination, 

(aa)S 4 - (ajî) »! 4- ( <*'/) ? 4- . . .= x, 

(MS + W)’i + (Pï)i; + ... = r, 

(•/“) 5 +{ïP)» + (77)?+" •=*;. 





( 80 ) 

et remarquons, avant tout, que les coefficients placés symé- 
triquement sont nécessairement égaux, c’est-à-dire que 

(P“) = (“?)» 

(7P) = (P7)> 


comme cela résulte de la théorie de l'élimination et comme 
d’ailleurs nous le démontrons plus loin 5 nous aurons 

x°— — (ctx)(al) — (*$)(bl) — (cty) [cl) , 

W »• =-(7«)(«0 -<7P)(«) “(77) W) — 


et ën posant 

! (o£0t ) X — f— ( oc(3 ) itL “H ( cty ) G 4- • • . 

(j3st) X 4- ((3(3) ift> 4 - (|3y) G -4- . . . 

(ya) X -f- (y(3) ifb 4- (yy) G 4 - . . . 

nous obtiendrons 

K = 4 — A(a/) — B(à/) — C(c/)— 
et si nous posons de plus , 

1 2 A 4 * b B — f— c C 4 * « * • — - p s , 
et' A 4- b B 4” c' C 4- . * ■ p ' z ' , 
a" A 4- i”B4-c"C 4-. ..=p"t", 

il viendra 

(n) K = A— lt — fi' — l“t" — 

9 . 

La comparaison des équations (7) et (9) nous apprend 
que les quantités auxiliaires A , B, C, etc., sont les valeurs 
que prennent les indéterminées x,jr t z , etc., lorsque l’on 
suppose 

£=**•> * = £ = G,...; 
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( 8 ‘ ) 

0 » en conclut évidemment, 

[aa) A -f- ( ab ) B («c) C -f-. . . = A, 

(ba) A 4 - (bb) B -+- (bc) C -+-• . . = Uh, 

(e«) A -f- ( cb ) B-f-(ce)C-f-...= G, 

En ajoutant les équations (io) après les avoir multipliées 

a a' a'' . 

par -» - 7 » - 77 » etc., on obtient 

r P P P 

A» — — a z — t— u z — f- c i s ” — r— . . . , 
et de môme 

US, = bt 4- b'e'-h b" s" -h... , 

G = ce -4- c'e' + c"z" . . , 


A, étant, comme on l’a dit, la valeur que prend X lorsqu’on 
substitue à v, v', v", etc., leurs valeurs observées; on aper- 
çoit facilement que si l’on applique à chacune de ces quan- 
tités les corrections — z, — s', — s", etc., la valeur deX de- 
vient égale à zéro, et que de môme Y, Z, etc., s’évanouis- 
sentpar cette hypothèse. L’équation (it) prouve aussi que K 
est la valeur que prend u par suite des mêmes substitutions. 

En nommant compensation des observations , l’applica- 
tion des corrections — e, — e', — e", etc. , aux grandeurs 
directement observées, nous sommes conduit à une con- 
séquence très-importante : 

Les observations compensées comme nous l'avons indi- * 
qué, satisfont exactement à toutes les équations de condi- 
tion, et font prendre à toute fonction des quantités obser- 
vées la valeur qui résulterait de la combinaison la plus 
convenable des observations non modifiées ; et puisque les 
équations de condition sont trop peu nombreuses pour 
qu’on puisse en déduire les valeurs exactes des erreurs, 
nous aurons au moins trouvé, par ce qui précède, les er- 
reurs les plus plausibles. C’est sous ce nom que les quan- 
tités e, t", etc., .seront désormais désignées. 

c. , 6 
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10. 


‘Le nombre des observations étant plus grand que celui 
des équations de condition , outre le système des correc- 
tions les plus plausibles, on peut en trouver un nombre 
infini qui rendent les équations de condition exactement 
satisfaites. 

Il importe d’examiner les relations qui lient entre eux 
ces divers systèmes. Soit — E, — E', — E", etc., un pareil 
système de corrections , autre que le système le plus plau- 
sible; nous aurons 

«E -+- a'E' 4- a v E"-h . . . = ,i>, 
bE 4- />' E' 4- i"E"-f- . . . = U!, , 
cE 4- e' E' 4- c"E"4- . . . = C , 


Multipliant ces équations par A , B, C, etc., et ajoutant, il 
vient, en ayant égard aux équations ( 10 ), 

ptE + /A'E' -4- p"t" E" = A.L-4-Bitl, H- Ce-t-.... 

Mais les équations (i3) combiuéesde la même manière don- 
nent 

(i4 )/>£* -f -/A'* +/)"t" , + ... = AA + Biti, CS 4-. . . ; 

de la combinaison de ces résultats on déduit facilement 
p F. 1 -4- p' E ,J 4- p" E" ! -4- . . .= p i- -4- p' s” 4- //'«"’ 4 -. .. 

-4- />(E - «)> 4 -p' (E'- *' }* *1- p"{ E" — «")’4-..., 

et, par suite, Ja somme 

p E’ 4- E' 1 4- p E" 2 4- . . . 

est, nécessairement, plus grande que 
pt 1 4- 4- />"*"’ 4- 

ce que l’on peut énoncer de la manière suivante : 

Théorème. — Les carrés des corrections qui peuvent 
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concilier les observations avec les équations dé condition 
étant respectivement multipliés par les poids des observa- 
tions correspondantes , donnent une somme minimum 
quand on adopte les corrections les plus plausibles . 

On reconnaît précisément le principe des moindres carrés 
dont, au reste, les équations (12) et (10) peuvent facilement 
se déduire. La somme minimum, que nous désignerons 
désormais par S, est égale, d’après l’équation ( r 4 ) , à 


11 . 

La détermination des erreurs les plus plausibles étant in- 
dépendante de /, Z", etc., fournit la préparation la plus 

commode, quel que soit l’usage que l’on veuille faire ulté- 
rieurement des observations. En outre, 011 voit sans peine 
que, pour atteindre ce but, il n’est pas nécessaire d'effec- 
tuer 1 ’éliminalion indéfinie, c’est-à-dire de calculer («a), 
(a/ 3 ), etc., il suffit de déduire des équations (12), par une 
élimination définie, les quantités auxiliaires A, B, C, etc., 
que nous nommerons, dans ce qui va suivre, les corrélatifs 
des équations de condition 

* X = o, Y=o, Z = o, ... j 

ces quanti tés seront ensuite substituées dans l’équation (tu). 

Celte méthode ne laisse rien à désirer lorsqu’on demande 
seulement les valeurs les plus plausibles des quantités four- 
nies par l’observation. Mais il en est autrement lorsqu’on 
désire, eu outre, le poids de chacune des valeurs trou- 
vées. Quelle que soit alors celle des quatre formules précé- 
dentes que l’on veuille employer, il est indispensable de 
connaître L, L', L", etc., ou t ce qui revient au même, 
x°, r°, z°, etc.; par cette raison, il sera utile d’étudier, de 
plus près, l’élimination qui fournit ces quantités et d’obtenir 
une méthode plus commode pour la détermination des poids. 

* 6 . 
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12 . 

Les relations qui lient les quantités dont nous nous 
cupons sont notablement simplifiées par la considération 
de la fonction indéfinie du second degré 
(aa)x’ 4 - 2 (ab ) xy 4 - i ( ac ) xz . . 

4- (bb)jr 2 4- 2 ( bc)yz -t- . . . 4- ( ce) z 7 + . . . , 

que nous désignerons par T. 

Cette fonction est évidemment égale à 
, g, (ax + by 4- cz 4-. ■ ■ )* . (a'x 4- b' y 4- c'z4- . • •)’ 

De plus , on a évidemment 

(i 6 ) T lx ■+■ ny-\- Çz -+-■ . 

et si, enfin, on exprime x , y, z , etc., au moyen des équa- 
tions (y), en fonction de £, «, f, etc., on aura 

T = («a) Ç 1 4- 2 («p) £« 4- 2 (ay) çÇ 4-. • . 4- (PP) 

4- 2 ( P 7 ) 4- . . . 4- ( 77 ) S 1 4- 

La théorie développée plus haut fournit deux systèmes 
de valeurs déterminées pour les'quàntités x,y, z, etc., 
£ , > 3 , £, etc. Le premier est 

x = x\ y = y*, z = z‘, . . . , 
l=—(al), n=—(bl) i l = — (cl), 

à ce système correspond la valeur 

T =(«)-î v - 

ainsi qu’on le voit en comparant à l’équation ( 16 ), la troi- 
sième forme du poids P, ou par la considération directe de 
la forme ( 4 ). 

Le second système de valeurs est 

x = A, J = B, z = C, 

£ = X , •« = !)!), £ = j 

la valeur correspondante de T est 

T = S , * 
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comme cela est évident par les formules (10) et (i 5 ), et 
ericore par les formules (14) et (i 6 ). 


• 13 . 

Nous devons, avant tout, faire subir à la fonction T 
une transformation semblable à celle qui a été indiquée 
(Theoria Motus , art. 182), et, avec plus de développe- 
ments, dans les Recherches sur Pal/as. 

Posons , à cet effet , 


<«•» =<“)-£$’ 


{6c, i) = (6c) — 


{«b) (ac) 
(aa) 


{bd, l)= (6rf)_(^lM ) , 

, , , , ( ac y (6c, i) : 


(cd, 2) = (cd) — 


(ac) (ad) (6c, 1) (bd, 1) 


("«) 


(bb, l). 


(m, 3) = m - w _ L - (S4, 

V ' (««) ( bb,t ) (cc, 2) 


et, ensuite (*) , 

* (66, \)y + (6c, 1) z 4- (bd, . .= 

(cc, 2) z - 1 - (cd, 2) ce Ç", 

(<M, 3 ) 


on aura 

P , . Ç’’ , 

T ~ (^) + [66~T) + (cc, 2) + (rfrf', 3 ) + • 1 • ’ 

et*/, <p w , etc., se déduiront de J-, r; , ç, etc., par les 

(* ) Dans les calculs précédents il suffisait de trois lettres de chaque séria 
pour faire apercevoir la loi des formules; il a paru nécessaire d’en faire ici 
figurer une quatrième, pour rendre l'algorithme plus manifeste. 

( Note de M. Gai ss,) 
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équations suivantes : 


(aa) ■’ 

Y" -y — (fi) t _ ( ég > ') 

(«a)* (W, .) 

( flrf ) g (W. ') , 

7 — t=Ô*— -CWTT)" ' 


(frf, a) 
(ff> 2 ) 


et l’on en tirera facilement toutes les formules utiles à notre 
but. Ainsi , pour déterminer les corrélatifs A , B, C, etc., 
nous poserons 

(ab) 


(, 8 ) 


Uï/ = tJt, ■ , , 

( aa ) 

e"=s- (ac) 


A, 


x - [hc ’ 

(aa) (bb, 1) ’ 


(aa) 


( bd , 1 ) 


(frf, 2) 


,ii ’ C", 

(W, 1) (cr, 2) 


et enfin A, B, C, D, etc., s’obtiendront par les formules 
suivantes , en commençant par la dernière : 

(aa) A-f- (rt£)B-)- (ac) C-f- (ad) D-f-. .. = J. , 

(bb, i)B-f-(ic, i)C-f- (bd, 1 ) D = 1 !!/, 

(*9) { ' (cc, 2 )C + (crf, a)D-f- 

(dd, 3)0+;.,=®"', 


Pour exprimer S, nous aurons la formule nouvelle 

JU ! ©"’ ®"’ 


(20) S = 


( aa)' ‘ (bb, 1) ^ (ce, 


(dd, 3 ) 


enfin , le poids P, qu’il faut attribuer à la détermination la 
plus plausible de la quantité u, sera donné par la formule. 

, ... _ i^Y _ (*',«)’ _ (frf, 2 y __ (di,3 y _ 

' P ' (aa) (bb, 1) (ce, 2) (dd,$) 
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dans celle formule, on fait 


(5.2) 


(W,. ) = («)- 

{cl, 7 .)=(cl) — 
(dl, 3 )=(dl)~ 


[“!>) ( al ) 

{a a) 

M («0 

(aa) 

{ad) (ni) 
(aa) 


(bc, i) (bl, i) 

(bb, I) ’ 

(bd, i) (bl, 1) _ (cd, 2) (cl, 2 ) 
(bb, 1) (cc, 2) 


Les formules (17), (21), dont la simplicité ne laisse rien 
à désirer, fournissent la solution complète de notre pro- 
blème. 


14 . 


Après avoir résolu le problème que nous avions en vue, 
nous allons aborder quelques questions secondaires qui 
éclaireront davantage l’ensemble de cette théorie. 

Nous chercherons, en premier lieu, s’il peut arriver que 
l’élimination qui fournit x, y, z , etc., on fonction de £ , 
V , £, etc., devienne, dans certains cas, impossible. Cela 
aurait évidemment lieu si les fonctions £, r,, f , etc., n’é- 
taient pas indépendantes les unes des autres. Supposons, 
pour un instant, qu’il en soit ainsi, et que l’une d’elles 
puisse s’exprimer en fonction des autres, de telle sorte que 
l’on ait la relation identique 

aÇ -+- fj>, -4- yï -+- . . « — O , 

a, (3, y, etc., désignant des noifibres déterminés. 

On aura alors 


et si 


x(aa) fi (ab) •+- 7 (ac) = o, 

a(<7é) + p(iè) -t-7(èc)-f-... = o, 
a ( tf c ) 4- (} ( 6<r ) + 1 ( cc )-+-... = o, 


nous posons 

a a + fi ê + yc-H. 

aa' •+■ fi b' 4- yc'-t-, . . = p' , 

aa" -\- fi b" -h yc"~{r. . . p" , 
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on en déduira 

o0 + «'0' + a"©"+.. = o, 

h& +- b' &' + b" . = o, 

c0+ e' 0' -+- c"%” -+- . . . = O , 

et, par suite, 

pW 4- p'<d ’’ -+• /j"©" 1 4- . . = o; 

/», p\ p" r etc., étant, parleur nature, tous positifs, cette 
équation exige 

© = o, 0' = o , ©" = o , . .. . 

Si nous considérons les différentielles complètes dX . , 
d Y, r/Z , etc., répondant aux valeurs de e, e', etc., 
immédiatement fournies par les observations, ces diilërcn- 
li elles 

adv 4- a' dv' 4- a" do" 4- . . . , 
bd? 4- b' dv 4~ b dv 11 4~« ■ ■ , 
cdv 4- c dd 4- e” dv" 4- • . . , 


d’après les résultats précédents, seront liées, les unes aux 
autres, de telle sorte qu’en les multipliant respectivement 
par «, P,y, etc., la somme des produits sera identiquement 
nulle, en sorte que, parmi les équations. 

X = o, Y = o, Z = o,.. ., 

il eu est une au moins <îuc r on peut regarder comme 
inutile , car elle sera satisfaite dès que les autres le se- 
ront. 

En examinant la question de plus près, on voit que 
cette conclusion n’est applicable qu’à des valeurs des 
variables infiniment peu ditlërentes de celles que fournit 
l’observation. Il j a , en effet, deux cas à distinguer : le 
premier est celui où l’une des équations 

X = o, Y=o, Z = o , . . 

est renfermée dans les autres d’une manière générale et 


Digitizedby Google 



( ây ) 

absolue, et peut, par conséquent, être supprimée; le se- 
cond est celui on , pour les valeurs particulières de e, 
e', e", etc., auxquelles se rapportent les observations, l’une 
des fonctions X, Y, Z, etc., X par exemple, acquiert une 
valeur maximum ou minimum, ou, plus généralement, une 
valeur dont la différentielle s’annule lorsque les autres 
équations restent satisfaites. 

Mais comme nous ne considérons pour nos variables que 
des variations dont les carrés soient négligeables, ce second 
cas (qui dans les applications ne se présentera que bien rare- 
ment) , pourra être assimi lé au premier, et l’une des équations 
de condition pourra être supprimée comme surabondante. 

Si les équations restantes sont indépendantes dans le sens 
que nous venons d’indiquer, on peut, d’après ce qui pré- 
cède, être certain que l’élimination est possible. Nous nous 
réservons, du reste, de revenir sur cette matière qui mé- 
rite d’être examinée comme subtilité théorique plutôt que 
comme question d’une utilité pratique. 

15 . ' . 

Dans le premier Mémoire, art. 37 et suivants, nous 
avons montré le moyen de fixer, à posteriori , d’une ma- 
nière très-approcliée, le poids d’une détermination. Si les 
valeurs approchées de a quantités sont fournies par des 
observations également précises, et qu’on les compare avec 
les Valeurs qui résultent pour elles des hypothèses les plus 
plausibles qu’on puisse faire sur les p éléments dont elles 
dépendent, on a vu qu’il fallait ajouter les carrés des dillé- 
rences obtenues, diviser la somme par cr — p, et que le 
quotient pouvait être regardé comme une valeur approchée 
du carré de l’erreur moyenne inhérente à ce genre d’ob- 
servations. 

Si les observations sont inégalement précises, la seule 
modification que l’on doive apporter aux préceptes précé- 
dents, consiste en ce que l’on doit multiplier les carrés des 
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différences par les poids respectifs des observations corres- 
pondantes, et l’erreur moyenne obtenue de cette manière 
se rapporte aux observations dont le poids est pris pour 
unité. 

Dans le cas actpel , la somme des carrés dont nous par- 
lons se confond évidemment avec la sommeS, et la diffé- 
rence tz — p avec le nombre a des équations de condition. 
Par suite, pour l’erreur moyenne des observations dont le 


poids est i, nous aurons l’expression y/-? et la détermi- 
nation sera d’autant plus digne de confiance que a sera plus 
considérable. 

Mais il est bon d’établir ce résultat, indépendamment des 
raisonnements du premier Mémoire; pour y parvenir, nous 
introduirons quelques notations nouvelles. Supposons 
qu’aux valeurs 


£ = «, n== b, Ç = c , . . . , 

répondent 

x — a, y-=p, Z= 

de sorte que l’ôn ait 

a = a (aa) -+- 6 (ajî) -t- c ( «7 ) -K . • , 
$ — a (a8) -f- b ( (ijï) -H c ( Py ) -+- • . i , 
7 = a («y ) -I- b ( Py ) + c { 77 ) 4- . . ; 

et, en outre, qu'aux valeurs 

\ " (i\ 7i — b , Ç = c j . . . , 

répondent 

x = a', X—P'i z = y ' j -*>; 
enfin, qu’aux valeurs 



répondent. 

* = «", y = P", Z 

et ainsi de suite. 
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La combinaison des équations (4) et ( 9 ) fournit 

A = ae 4- a! e' 4- a" e" 4 - . . . , 

B= Çc + pe' + p"c" + ..., 

C — ye + y'e' 4 - y" e" 4 -. . . , 

et, comme on a 

S - — A» A 4 ~ tïî> B 4" Q G “4" • • • , 

on aura 

S =(ne 4 - a' e' 4 - a" e" 4- . . . ) ( a e 4- a ' (f 4- a" e" 4- ■ . • ) 

4 - (i<? 4 - ÎV 4 - 4 V + ...) (pe + pV4-f«"+...) 

4- (ce 4- eV 4- e" e" 4 - . . . ) ( 71 ? 4- 7 V 4- y" c" 4 - . . .) 4 -.... 

16 . 

La série des observations qui fournissent les quantités 
n, v' , v", etc., affectées des erreurs fortuites e, e', e", etc., 
peut être considérée comme une épreuve qui ne fait pas 
connaître, il est vrai, la grandeur de chaque erreur, 
mais qui, par le moyen des règles exposées plus haut, per- 
met de déterminer la quantité S, fonction connue de toutes 
les erreurs. Dans une telle épreuve, les erreurs peuvent être 
les unes plus grandes, les autres plus petites; mais plus 
sera grand le nombre des erreurs employées, plus il y aura 
une grande probabilité que S diffère peu de sa valeur 
moyenne : la difficulté revient donc à trouver la moyenne 
de S. 

Par les principes exposés dans le premier Mémoire et 
qu’il est inutile de reproduire ici , on trouve pour cette va- 
leur moyenne 

(«a 4- b -t- cy 4- .. .)m’ 4- («'a' 4- è'ff 4- c' y' 4 - . .) ru' 7 4- • • ■ . 

En nommant p l’erreur moyenne qui coi respond aux obser- 
vations dont le poids est 1 , de telle sorte que l’on ait 

f * 3 = pm* p'm '* = p"m" t — . . . , 
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l’expression précédente peut s’écrire comme il suit : 

fa « n'a' a"«" \ , /AS A' 8' 6" 6" \ , 

_ + _l r f... U>+ M £- +••• U* 

X/» P P \p P P 


■+ 


£7 , f_7_ , 

/> P' P" 

Mais on a trouvé 

a cl a' a! a" cl” 

1 t S h. 

P 


+ 


F- 


.= (an) (aa) 4- (ab) (afi) -f- (ac) («•/) 


P P P 

or, le second membre est l’unité, comme on le reconnaît 
facilement par la comparaison des équations (6) et (7). 

On trouvera de même 

AS b'W , b" fi" 

— H 7" + — 7T-+ ■ •= I, 

P P P 


c' f' 


C"l" 


CJ_ 

P P’ ’ P" 
et ainsi de suite. 

D’après cela , la valeur moyenne de S devient ou’, et si 
l’on juge permis de regarder la valeur fortuite de S comme 
égale à la valeur moyenne, on en conclut 

/S 

<r 


^v/1 


17. 

On peut apprécier la confiance que mérite cette détermi- 
nation en calculant l’erreur moyenne à craindre, soit 
pour sa valeur propre, soit pour celle de son carré. La se- 
conde sera la racine carrée de la valeur moyenne de l’ex- 
pression 



dont le développement s’obtiendra par des raisonnements 
semblables à ceux qui ont été exposés dans le premier Mé- 
moire (art. 39 et suivants). Nous les supprimons pour abré- 
ger, en nous contentant d’indiquer le résultat. 
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-L’erreur moyenne à craindrç dans la détermination du 
carré p* s’exprime par 



v* étant la valeur moyenne des quatrièmes puissances des 
erreurs dont le poids est l’unité, et N la somme 

( e x -j- b ( 8 -H cy -K - • . )* -f- (a' x' -t- b' J5’ -f- c' y* -f- , . ) * 

(a" x" -t- b" fi" -t- c" 7 " + ...)’+••• - 

Cette somme ne peut pas, en général, se simplifier; mais, 
par un procédé analogue à celui dont on a fait usage au pa- 
ragraphe 40 du premier Mémoire, on peut montrer que sa 

valeur est comprise entre a et Dans l’hypothèse sur la- 
quelle nous avions primitivement établi la méthode des 
moindres carrés , le terme qui contient cette somme dis- 
paraît à cause de 

V* = 3 (*', 

et la précision que l’on doit attribuer à la détermination 



• , 1 » 

est par conséquent la meme que si l’on avait opéré sur <7 
observations entachées d’erreurs exactement connues, con- 
formément aux préceptes des art. 15 et 16 du premier 
Mémoire. 

18. 

Pour la compensation des observations, il y a, comme 
nous l’avons dit., deux opérations à exécuter: première- 
ment, il faut déterminer les corrélatifs des équations de 
condition, c’est-à-dire A, B, C, etc., qui satisfont aux équa- 
tions ( 12 ) ; secondement, substituer ces quantités dans l’é- 
quation (io). La compensation ainsi obtenue peut être dite 
parfaite et complète , par opposition à la compensation 
imparfaite ou incomplète. Nous désignerons de cette der- 
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nièrc manière celles qui résultent des mêmes équations (io), 
dans lesquelles on substituera des valeurs de A, B, C qui ne 
satisferont pas aux équations (is), c’est-à-dire qui satisfe- 
ront à quelques-unes seulement, ou à aucune. .Nous ne nous 
occuperons pas ici d'un tel système de corrections, et nous 
ne leur accorderons même pas le nom de compensation. 

Lorsque les équations ( 10 ) sont satisfaites, les sys- 
tèmes ( 12 ) et (i3) deviennent équivalents, et la différence 
dont nous parlons peut alors s’énoncer comme il suit : Les 
observations complètement compensées satisfont aux équa- 
tions de condition 

X — o, Y = o, Z — o, ... ; 

les observations incomplètement compensées ne satisfont 
qu’à une portion de ces équations, et peut-être à aucune; 
la compensation à la suite de laquelle toutes les équations 
sont satisfaites, est nécessairement complète. 

19. 

Il résulte de la définition même des compensations, que 
la réunion de deux systèmes de compensations peut en four- 
nir un troisième, et l’on voit qu’il importe peu que les 
règles données pour obtenir une compensation parfaite 
soient appliquées aux observations primitives ou atix ob- 
servations déjà imparfaitement compensées. 

Soient — 0, — 0', — 0", etc., un système de compen- 
sations incomplètes, résultant des formules 

Bp +B'i H-Pc + .., . 

. ^ e'p' — A -t- B "b' -h Ce' 

1 ' j B" p" — A» fl" 4- B \b" -t- Ce" + . . , 

Les observations ainsi changées ne satisfaisant pas à toutes 
les équations de condition, soient A>*, 1 )!)*, C 11 , etc., les va- 
leurs que prennent X, Y, Z, etc., quand on y substitue les 
valeurs ainsi obtenues pour f, e', v", etc. On devra cher- 
cher les valeurs A*, B*, C*, etc., satisfaisant aux équa- 
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A* = A* (aa) 4 - B* (ab) 4 - C* (ac) 

' , \ ** = A* (ab) -4- 15* (b b) -f- C* (4c) 

1 1 | ©* = A*(«c) 4-B*(Ac)4-C*(cc) +..., 

et cela fait, la compensation complète des observations 
ainsi modifiées se fera par les nouveaux changements — x, 
— — r.", etc.; x, x', x", etc., se déduisant des for- 

mules 

/ y. p = A* a -h B* b 4 - C*c 4-' '. . . , 

, | y p' “ A* a' 4- R* b' 4- C.*c’ 4 - . • . , 

( 1 *y=A t /4-B*J*4-tV4-, ., 

Cherchons comment ces corrections s’accordent avec la com- 
pensation complète des observ ations primitives. Il est clair 
d’abord que l’on a 

A * .= a — «6 — a!Q' — «"0" — . . . , 

Di* = l(î> — £0 _ b' 0' — 4" 0" , 

s* = e — c© — c'0' — c"0" — . . . , 


En substituant dans ces équations, pour 0, 0', 0 ff , etc., 
leurs valeurs fournies par le système (I), pour A.*, D!>*, 
©*, etc., celles que donne le système (II), il vient 

A = (A 0 -+- A*) (aa) -h (11° -t- H*) (ab) 

vit, = (A 0 - 4 - A*) ( ab) 4 - (B° 4- B*) (44) + . . . , 

e = (A 0 4- a*) (ac) 4- (B* 4- B*) (4c) 4 -. . . ; 

et il suit de là que les corrélatifs des équations de condi- 
tion ( 12 ) sont 

A = A» 4- A*, B = B» 4- B*, C = C° 4- C*, . . . , 
et alors les équations ( 10 ), (I) et (III) montrent que l'on a 
i = 64-«, t' — 4- s” = 0" 4 -r.",...; 

et, par suite, la compensation parfaite a la même valeur 
pour chaque inconnue, soit qu’on la calcule directement, 
soit qu’on l’obtienne médiatement en partant d’une com- 
pensation incomplète. 
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20 . 


Lorsque les équations de condition sont très-nombreuses, 
la détermination des quantités corrélatives A, R, C, ete., 
peut exiger des calculs tellement longs, que le calculateur v 
en soit rebuté-, il pourra être avantageux d’obtenir, dâns 
ce cas, une compensation complète, à l’aide d’une série 
d’approximations reposant sur le théorème de l’article pré- 
cédent. On partagera, pour cela, les équations de condition 
en deux ou plusieurs groupes, et l’on cherchera d’abord une 
compensation qui rende satisfaites les équations du pre- 
mier groupe. On traitera ensuite les valeurs modifiées par . 
ce premier calcul , et on les corrigera de nouveau en ayant 
égard seulement aux équations du second groupe. Ce se- 
cond calcul donnera des résultats qui , en général, ne satis- 
feront plus aux équations du premier groupe, et il faudra, 
si l’on n’a formé que deux groupes, revenir alors au premier 
ety satisfaire à l’aide de nouvelles corrections. Les observa- 
tions seront ensuite soumises à une quatrième compensa- 
tion , dans laquelle on n’aura égard qu’aux conditions du 
second groupe; et en opérant ainsi alternativement sur l’un 
et l’autre groupe d’équations, on formera des corrections 
qui seront nécessairement de plus en plus petites. Si le 
choix des groupes a été fait habilement , on arrivera bien 
vite à des valeurs que les corrections ultérieures ne chan- 
geront plus. 

Quand on forme plus de deux groupes, on. doit procéder 
de la même manière, les divers groupes étant employés suc- 
cessivement jusqu’au dernier, après quoi on revient au pre- 
mier pour les reprendre dans le même ordre. Il nous suffit 
d’avoir indiqué ce procédé, dont la réussite dépendra beau- 
coup de l’habileté du calculateur. 

21 . 

Il nous reste à donner la démonstration du lemme admis 
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dans l'art. 8. Adoptons, pour plus de clarté, des notations 
plus propres à mettre la démonstration en lumière. 

Soient x°, x', x", etc., des indéterminées; supposons 
que les équations 

«•«j* 4- n" x 1 4- n”x" 4- ■ • • = X% 
n'°x’ -I- n" x’ + n ” x” -t- . . . = X', 

«"x" 4- * !l x' + n™x" 4- . . .= X", 


aient donné, par élimination, les suivantes : 

N" X* -+- N 41 X' 4 - N” X" -t- .—x*, 
N'*X*4- N"X'4- N ,, X"4-... = x', 
S : *X*+ N” X' 4- N” X" 4 x”, 


En substituant, dans les deux premières équations du se- 
cond système, les valeurs de X°, X', X", etc., fournies par 
le premier, nous obtiendrons deux équations identiques : 

x • = N“° («•• x° 4- «"'x' + n al x" 4- . . .) 

4 - N” 1 (n l# x • + n'*x" 4 - . . .) 

4 - N" (n» x" 4- n n x' 4- n” x" 4- . . . ) 

4- N ,s («“x* 4- n 51 x' 4- n”x" 4- . .) 

4- 

x ' — N 1 * (n‘° x* 4- «“'x' 4- n n x" 4-. . .) 

4 - N" («“x" 4- «" x' + n"x" 4-.- .) 

4- N 11 (« M x* 4- n n x' 4- «"x" 4-. . .). 

H~ • 

Ces équations étant identiques , on peut y substituer telles 
quantités que l’on voudra à la place de x”, x', x", etc. Fai- 
sons dans la première 

x* = N 1 *, x' = N", x" = N 1 ’, . . , 

et dans la seconde 

x* = N M , x' = N 01 , x" == N t! , . . . . 

Eu retranchant ensuite lesdeux identités membre à membre, 
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il viendra 

N" — N* 1 = (N M N" — N“ N“) (*»' — n'”) 

-t- (N** N 1 ’ — N" N”) (/*»’ — n ») 

-4- (N" N' 3 — N' 0 N°>) («” — n") 

■+■ 

-t- (N 01 N 13 — N" N 03 ) (fl 13 — n 31 ) 

■+■ (N 01 N 13 — N" N 33 ) (/J 13 — n 3 ') 

-t- : . . . . 

-+■ (N" 3 N 13 — N” N' 3 ) (n 33 — n 33 ) 

4 - 

Ce que l’on peut écrire plus brièvement de la manière sui- 
vante : 

n "> — n » 1 =2 ( N ’“ N,/3 — N ' a N,i3 ) (» a/3 — ni3a ) ; 

or, (3 désignant deux indices pris au hasard ; on conclut de 
là que les égalités • 

/J 31 = «'•,. * 03 =« M , 

et, généralement , 

entraînent nécessairement 

N 1 * = N 01 ; 

et, comme l’ordre des indéterminées est arbitraire, il est 
évident que, dans la supposition admise, on aura généra- 
lement 

N*/ 3 — 

22 . 

La méthode exposée dans ce Mémoire devant surtout s’ap- 
pliquer utilement aux calculs de haute géodésie, le lecteur 
nous saura gré d’y joindre quelques exemples puisés dans 
cette partie de la science. 

Les équations de condition qui existent entre les angles 
d’un système de triangles, peuvent, en général, se par- 
tager en trois catégories. 

1. La somme des angles horizontaux, formés autour d’un 
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même sommet et embrassant la totalité de l'horizon , doit 
être égale à quatre droits. 

II. La somme des angles de chaque triangle peut toujours 
être regardée comme connue ; car, lors même que le triangle 
est situé sur une surface courbe, l’excès de la somme de ses 
angles sur deux droits peut être calculé avec une telle 
approximation, qu’il est permis de considérer le résultat 
comme absolument exact. 

III. Enfin, on obtient un troisième genre de relations en 
examinant les rapports des côtés dans les triangles qui for- 
ment un réseau fermé. Si, en effet, les triangles sont tel- 
lement placés, que le second triangle ait un côté a commun 
avec le premier, et un côté b commun avec le troisième; 
si le quatrième triaugle a deux côtés c et d respectivement 
communs avec le troisième et le cinquième, et ainsi de 
suite, jusqu’au dernier triangle, qui ait avec le précédent 
un côté commun h , et avec le premier de tous un coté com- 
mun l , les quotients 


a b c d l 



pourront se calculer au moyen des angles qui leur sont 
opposés dans le triangle dont les deux côtés comparés font 
partie , et comme le produit de ces fractions est évidemment 
l’unité, on aura une relation entre les sinus des divers 
angles mesurés (diminués du tiers de l’excès sphérique ou 
sphéroïdique lorsqu’on opérera sur une surface courbe). 
Du reste, dans les réseaux un peu compliqués, il arrive 
souvent que les équations de la seconde et de la troisième 
catégorie rentrent en partie les unes dans les autres, et que, 
par suite, leur nombre doit être réduit. Au contraire, il 
pourra arriver, mais seulement dans des cas assez rares, 
que l’on adjoigne quelques équations nouvelles à celles de 
la seconde catégorie; c’est ce qui aura lieu lorsque le ré- 
seau contiendra des polygones non divisées en triangles; ou 
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pourra alors introduire des équations relatives aux figures 
qui ont plus de trois côtés. Dans une autre occasion , nous 
reviendrons avec plus de détails sur ces diverses circon- 
stances, dont l’examen nous éloignerait en ce moment de 
'notre but. iVous ne pouvons cependant nous dispenser de 
faire ici une remarque indispensable à ceux qui voudraient 
faire l’application rigoureuse de notre théorie : nous suppo- 
sons toujours que les quantités désignées par o, e', v", etc., 
ont été observées immédiatement, ou déduites d’observa- 
tions telles, que leurs déterminations soient indépendantes 
les unes des autres, ou, au moins, puissent être regardées 
comme telles. Dans la pratique la plus ordinaire, on observe 
les angles que l’on peut regarder comme étant les éléments 
v, v\ v" y etc., eux-mêmes. Mais on ne doit pas oublier que 
si le système contient, en outre, des triangles dont les angles 
n’aient pas été directement observés et aient été déduits de 
ceux que l’on connaissait, par des additions ou soustrac- 
tions, ces angles ne devront pas être mis au nombre des 
grandeurs déterminées par l’observation, et l’on devra les 
faire entrer dans le calcul comme des fonctions des éléments 
qui ont servi à les former. Il en sera autrement si l’on 
adopte la méthode d’observations de M. Struve ( Astrono - 
mische Nachrichten, II, page 43 1 ), qui consiste à détermi- 
ner toutes les directions autour d’un même sommet , en les 
rapportant toutes à une seule et même direction arbitraire. 
Les angles mesurés ainsi seront pris alors pour o, y", etc., 

et les angles des triangles se présenteront tous comme des 
différences. Les équations de la première catégorie devront, 
dans ce cas, être supprimées comme superflues, car elles 
seront identiquement satisfaites. Le procédé que j’ai suivi 
moi-même dans les triangulations exécutées pendant ces 
dernières années, diffère des deux méthodes précédentes; 
on peut cependant l’assimiler, quant au résultat, avec le 
procédé de M. Struve, en ce sens que, dans chaque station , 
on doit regarder e, m', u", etc., comme les angles formés 
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par les directions qui en partent, avec une même ligne 
arbitrairement choisie. 

Nous donnerons deux exemples : le premier se rapporte 
au premier mode d’opération , et le second est relatif à 
des observations faites d’après la seconde méthode. 

23 . 

Le premier exemple nous sera fourni par l’ouvrage de 
M. Kraycnhof : Précis historique des opérations trigono- 
métriques faites en Hollande. Nous chercherons à com- 
penser la partie des observations relatives au terrain com- 
pris entre Harlingen , Sneek , Oldeholtpade , Ballum , 
Leeuvvarden, Dockum , Drachten , Oostenvolde et Gronin- 
gen. Entre ces points , on a formé neuf triangles numérotés,- 
dans l’ouvrage cité , i ai , 122, 123,124, 1 25 , 127, 128, 
i 3 i, i 3 a. Les angles observés sont les suivants : 

Trianclf. 121. 

0. Harlingen 5 o" 58 ' i 5",238 

1. Leeuwarden 82.47.1 5 , 35 1 

2. Ballum... 4 ^ 1 4 - 2.7 , 202 

TnisNCLE 122. 

3 . Harlingen 5 i . 5.39,717 

4 - Sneek 70 - 48 . 33,445 

5 . Leeuwarden 58 . 5.48,707 

TkIAHCLF. 123 . 

6. Sneek 49 3 o- 4 o,o 5 i 

7. Drachten 4 2 - 52. 5 g, 382 

8. Leeuwarden 87.36.21,057 

Thianclf. 124 . 

9 Sneek 45 - 36 . 7 , 49 ® 

10. Oldeholtpade 67.52. 0,048 

11 . Drachten 66.3 i.56,5i3 
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Triangle 

1 25 . 

12. 

Drachten 

53 ° 55 ' 24", 7 45 

i 3 . 

Oldeholtpade 

4 7 -48.52,58o 

» 4 - 

Oosterwolde 

78. 15.42,347 


Tri.\nc,i.f. 

I27. 

i 5 . 

Leeuwarden 

59.24. o ,645 

l6. 

Dockum 

76.34. 9,021 

> 7 - 

Ballum 

44 - i. 5 i,o 4 o 


Triangle 

128. 

18. 

Leeuwarden 

72. 6 32,043 

• 9 - 

Drachten 

46.53.27 , i 63 

20. 

Dockum 

61. 0. 4 ,494 


Triangle 

1 3 r . 

21 . 

Dockum 

5 7 1.55,292 

22. 

Drachten 

83 . 33 . i 4 , 5 i 5 

23 . 

Groningen 

39.24.52,397 


Triancle 

l 32 . 

24 * 

Oosterwolde 

81.54.17,447 

25 . 

Groningen. . 

3 i .52.46,094 

26. 

Drachten 

66.12.57, 246 


La considération de ces triangles montre que les vingt-sept 
angles, directement fournis par l’observation, ont entre eux 
treize relations nécessaires, savoir : deux de la première 
espèce, neuf de la seconde, et deux de la troisième. Mais 
il n’est pas utile d’écrire ici toutes ces équations sous 
leur forme finie, car pour le calcul nous avons besoin seu- 
lement des quantités, désignées dans la théorie générale, 
par X, a, a', a '',. . . , aPo, b, b ', b",. . . , c’est pourquoi 
nous écrirons immédiatement les équations (i3), qui met- 
tent ces quantités en évidence. Au lieu de e, ê', e", etc., 
nous écrirons simplement ici (o), (i), ( 2 ), etc. De cette ma: 
nière, aux deux équations du premier genre répondent les 
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suivantes : 

(.)+( 5 )+( 8 ) + (. 5 ) +(.8) = - a",, 97, 

( 7 ) + ( 11 ) + ( ,2 )+ ( l 9 ) + ( 22 ) + ( ? - 6 ) = — o", 436 . 

Nous trouvons ensuite, pour les excès sphéroïdiques des 
neuf triangles : i", 749 ; i'\i 47 ’, i", 698 ; o",873 ; 

i",i 6 y, i",i °4 ; a",i6i ; i", 4 ° 3 . Nous aurons alors lequa- 
tion de condition du second genre : 

(if) + p<‘) ■+■ (>(’) — i8o°,o' i",749 — °» 

et ainsi des autres, et nous avons les neuf équations sui- 
vantes : 

(o) +(.) + (*) = -3' >9 58, 

(3 ) + ( 4 ) + ( 5 ) = +0,722, 

(6) + ( 7 >+ (8)=-o, 7 53, 

(9) ( ,o ) + (n) = + 2,355,. 

(ta) + (i3) + (>4) = — 1 >aoi , 

(i 5 ) + (16) + (17) = — 0,461 , 

(•8) + (19) + (20) = +2,596, 

(21) + (22) -t- (a 3 ) = + o,o 43 , 

(24) + ( 25 ) -4- (26) = — 0,616. 

Les équations de condition du troisième genre s’expriment 
plus facilement par le moyen des logarithmes : la première 
est ' 

log sin (•>' — o", 583 ) — log sin (v' — o", 583 ) — log sin (v 5 — o",382) 
- 4 - log sin ( v ’ — o",382) — log sin (<4 — o", 4 i 4 ) 
-t- log sin (<>’ — o ",4 1 4 ) — log sin (v 18 — o", 38 g) ’ 
-4- log sin (fi”— o", 38 g) — logsin (v”— o'', 368 ) 
- 4 - log si n (v M — o" , 368 ) = o . 

Il semble inutile de développer l’autre sous forme finie. A 
ces deux cqualions répondent les suivantes, dans lesquelles 
les coefficients se rapportent à la septième décimale des lo- 
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garithmes vulgaires (*) : 

17,068(0) — 20,174(2)— 16,993 ( 3 ) -h 7,328(4)— 17,976(6) 
+ 22,672(7)— 5 , 028(16) +21, 780(17)— 19,710(19) 
+ 11,671(2.0) = — 371, 

1 7 ,976(6) — 0,880(8) — 20,617(9)+ 8 , 564 (io) — rg, 082(13) 
+ 4,375 (1 4) + 6,798(18)— 1 1 ,671(20) + 1 3 , 657(2 r) 
— 25 ,620 ( 23 ) — 2 ,9g5 (24) + 33 , 85 | ( 25 ) = + 370. 

Aucune raison ne nous portant à attribuer des poids iné- 
gaux aux diverses observations, nous supposerons 

p* = p' = p 1 = , . . = 1 . 

En désignant les corrélatifs des équations de condition dans 
l’ordre même où ces équations ont été écrites, par 
A,B,C,D,E,F,G,H,I,K,L,M,N, 
nous les déterminerons par les équations suivantes : 


~ 436 = 6 KL" M 

— 3 ,g 58 = A + 3 C — 3 ,io 6 ftI, 

+ 0,722 = A + 3 D — 9,665 M, 

— 0,753 = A + B + 3 E + 4 , 6 g 6 M + 17,0961^, 

+ 2,355 = B + 3 F — i2,o53 N, 

— 1,201 = B + 3 G — > 4 » 7°7 N, 

— 0,461 = A + 3 H + 16,752 M, 

j- 2 , 5 g 6 = A + Iî + 3 1 — 8,039 M — 4,874 N, 

+ o,o 43 =B+ 3 K — 11, g 63 N, 

— 0,616= B + 3 L+ 3 o, 85 gN, 

— 371= 2,962 B — 3 , 106 C — 9,665 D + 4 ,696 E 

+ 16,752 H — 8 ,o 3 g I + 2go2 ,27 M 

— 459,33 N, ■■ 

t~ 370 = 5,917 A + 17 ,096 E — 12 ,o 53 F — 14,707 g 

— 4,8741 — 1 1 ,g 63 K + 3 o, 85 g L 

— 459,33 M + 3385 . 96N 


(*) Ces coefficients ont été tous multipliés, après la différentiation, par 

. n et o l8o.Go.Go 

10, et divises par 200304, o— » pour convertir les erreurs en 


secondes. 


J. B. 
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Nous en déduisons par l'élimination : 

A = — 0,598 
B = — o ,255 

c = — 1 , 234 

U = + o, 086 
E = — 0,447 
F = 4 - 1 ,35i 
G = 40,271 

Les erreurs les plus plausibles sont enfin données par les 
formules 

(o) = C 4- 17,068 M 
• (i)=A + C 

(2) = C — 20, 174 M 

( 3 ) = D — 16,993 M, 


et nous obtenons les valeurs numériques suivantes , aux- 
quelles nous adjoignons, pour qu’on fasse la comparaison, 
les corrections adoptées par M. Krayenhof : 

' DE Kr. DE Kr. 

(0) = — 3", 108 —2", 090 (i4) = -t-o",795 4-2", 4°° 

(1) = — i ,832 4-0, 116 ( i 5 ) = 4- o , o 6 i 4 - 1,273 

(2) =4- 0,981 —1,982 (i6) = 4- 1,211 4-5,945 

(3) =4- 1,952 4- 1,722 I (17)=— 1,732 —7,674 

(4) =—0,719 4-2,848 (i8)=4 i,265 4-1,876 

( 5 ) = — o,5i2 — 3,848 (49) = 4- 2,959 4- 6 , 25 i 

(6) =4- 3,648 — 0,137 (20) = — 1,628 — 5,53o 

(7) = — 3,221 4- 1,000 {21 ) — 4" a,2ii 4- 3,486 

(8) = — l,l8o — I , 6l 4 ( 22 ) = 4- 0,322 — 3,454 

(9) =— 1,116 Q (23)=— 2,489 o 

(10) = 4- 2,376 4-5,928 (24) = — 1,709 4- o,4oo 

(11) =4- 1,096 — 3,570 (25) = 4- 2,701 4-2,o54 

(12) = -+- 0,016 4-2,4l4 (26) = — 1,606 —3,077 

(13) = — 2,oi3 — 6,014 

La somme des carrés de nos corrections estQ7,8845 ; l’erreur 
“ moyenne, telle que l’indiquent les 27 angles observés, est 
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Il = 4- o, 65 g 
I=4-i ,o 5 o 
K = 4- 0,577 
L = — 1 , 35 1 
M = — o , 109792 
N = 4- o, 1 19681. 
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par conséquent, 

v/® = ,,, 44 , 

La somme des carrés des corrections de M. Krayenhof, 
est de 34i ,4^oi . 

24 . 

Les triangles dont les sommets, dans la triangulation du 
Hanovre, ont été placés à Falkenberg, Breithorn, Hausel- 
berg, Wulfsode et Wilsede, nous fourniront un second 
exemple. 

On a observé les directions suivantes : 

A la station de Falkekbebg. 

0. Wilsede 1 87° 47 ' 3 o" 3 n 

1. Wulfsode 225.9.39.676 

2. Ilauselberg 266. 1 3 . 56 . 2.39 

3 . Breithorn .. 274 . i 4 - 43-634 

A la station de Beeituoen, 

4 - Falkenberg g4-33.4o.755 

5 . Hauselberg 122. 5 1.23 o 54 

6. Wilsede i 5 o.i 8.35 100 

A la station de Hauselbebg. 

7. Falkenberg 86.29. 6.872 

8 . Wilsede 154.37. 9 624 

g. Wulfsode.. 189. 2.56.376 

10. Breithorn 302.47-37.732 

A la station de Wulfsode. 

11. Hauselberg 9. 5 . 36 . 5 g 3 

12. Falkenberg. 4^-27 . 33.556 

1 3 . Wilsede 1 18. 44 - * 3 . 1 5 g 

A la station de Wilsede. 

14. Falkenberg 7 . 5 1 . 1.027 

1 5 . Wulfsode....... 298.29.49.519 

16. Breithorn 33 o. 3 . 7*3g2 

17. Hauselberg 334 «2.5. 26. 746 
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Ces observations permettent de former sept triangles. 


Tri angle I er . 


Falkenberg 

8 ° o' 4 7 ">395 

Breithorn . . 


Hauselberg. 



Triangle H. 

Falkenberg . 

86.27.13,323 

Breithorn . . 

55.44.54,345 

Wilsede. . . . 



Triangle III. 

Falkenberg. 


Hauselberg 


Wulfsode.. . 



Triangle IV. 

Falkenberg. 

. ... ' 78.26 25,928 

Hauselberg. 


Wilsede. . . . 



Triangle V. 

Falkenberg. 


Wulfsode . . , 


Wilsede . . . . 



Triangle VI. 

Breithorn . . . 


Hauselberg . . 

OC 

0 

w 

cc 

5 

cc 

Wilsede 



Triancle VII. 

Hauselberg . . 


Wulfsode . . 


Wilsede 




Fions avons ici sept écjualious de condition du second 
genre (il n’y a pas lieu évidemment d’en former du premier 
genre) -, pour les former, nous devons chercher, avant tout, 


( l<> 8 ) 

1rs excès sphéroïdiques des sept triangles, et pour cela il est 
indispensable de connaître la longueur d’un côté. Celui 
qui réunit Wilsede à Wulfsodeest 228^7 m ,94- On en con- 
clut, pour les excès sphéroïdiques des divers triangles : 
I...o",202; II.-. .2 ,/ ,44 2 j III - - • i' / ,2D7 ; IV. . . i",qi 9 ; 
V...i",y57; VI...o",32i; VII . . . i", 2 9 5. 

Si l’on désigne par F* 01 , e (,) , e! 1 *, «' <3) , etc., les angles qui 
déterminent les directions indiquées plus haut, et mar- 
quées des mêmes indices, les angles du premier triangle 
seront 

- ii(J) — ,»(•) , p( s ! — pi*), 36o -4- «(’) — 

et la première équatiou de condition est par conséquent 

— et 1 ) -4- p( s ) — «>(') •+• -t- — p('*) -t- 1 79 0 5q' 59", 798 = o. 

Les six triangles restants fourniront six équations analo- 
gues, mais un peu d’attention montrera que ces équations 
ne sont pas indépendantes; la seconde est en efl'et identique 
avec la somme de la première, de la quatrième et de la 
sixième; la somme de la troisième et de la cinquième est 
identique avec celle de la quatrième et de la septième : c’est 
pourquoi nous négligerons la seconde et la cinquième. Au 
lieu des équations restantes sous forme finie, nous écrirons 
ici les équations correspondantes ( 1 3 ) , en subsi tuant aux 
notations e, e', t " , etc., (o), (1), (2), etc. : 

-«", 368 = -(2)+ ( 3) — (4) -4- ( 5 } -4- (7 ) — (.o), 

-t- 1,773 = -(i)-I- (2) — ( 7 )-»- (9 ) — (>•) + (<*). 

■ -4- 1,042 = — (0)4- (2) — (7) -4- (8)4- (,4) -(! 7 ), 

— 0,8.3 = — (5)4-(6) — (8)4- (10) — (.6) 4- {«7) » 

— o.75o = — (8) 4- (9) — (11)4- (« 3 ) — (i 5 ) 4- (17).' 

Un peut obtenir, au moyen des triangles du système, 
huit équations du troisième genre, et pour cela il est per- 
mis de combiner trois des quatre triangles I, II, IV, VI, ou 
des triangles III, IV, V, VII ; cependant un peu d’attention 
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montre qu’il suffit d’en considérer deux appartenant res- 
pectivement aux deux systèmes de triangles et que celles-là 
comprendront toutes les autres. 

Nous aurons ainsi , pour sixième et septième équation de 
condition , 

Iogsin(i>C s ) — v( ] ) — o",o 6 ']) — logsin(p( i ) — — o",o67) 

+ logsin (e( H > — — o",64o) — logsin (pl'l — 1>(‘) — o",64o) ' 

+ logsin (eW — v(*> — o", 107) — logsin (v < v > — •>('•) — o",ioj) ~o, 

logsin («(“î — 1>(‘> — o",4i9) — logsin (pf ■*> — (.(>') — o",4'9) 

-i- logsin (r< M ) — (,(") — o",64o) — logsin (pW — p<") — o",64o) 

-f logsin (p( IJ ) — p("> — o", 43 2 ) — logsin (et 1 ’)— «PO — o", 43 2 )=o, 

auxquelles répondent les équations 

+ 2 5 = + 4 , 3 i(o) — 1 53, 88 (2) -t- 149.57 ( 3 ) 

-+■ 39,11 ( 4 )— 79 > 6 4 ( 5 )-+- 4 °. 53 ( 6 ) 

-t- 3 i,go (i 4 ) -4- 2 . 75 , 3 g(i 6 ) — 307, 5.9(17', 

- 3 = -|- 4 , 3 i(o) -24,16 (.) + 19, 85 (a) ' 

+ 36 ,i 1 (11) — 28,59(12) — 7,52 (i 3 ) 

• +31,90(14)+ 29,06(15)— 60,96(17). 

Si nous attribuons la même certitude aux diverses direc- 
tions, en supposant /> (0) = p (,) = p (,) ... = 1 ,les corrélatifs 
des sept équations de condition étant désignés par A , B, 
C, D, E, F, G, leur détermination dépendra des équations 
suivantes : 

— i ,368 =r + 6 A — 2B — 2 C — 2 D + 184,72F — 19,85 G, 
+ 1,773 = — 2A+6B + 2C+ 2 E — 1 53, 88 F — 20,696, 
-t- 1 ,042 =— 2 A + 2B + 6C — 2 D — 2E +■ 181,00F 

+ io 8 , 4 o G, 

— o, 8 i 3 = — 2 A — 2C + 6D+2E — 462,61 F — 60,96 G, 

— 0,750 =+ 2 B — 2 C 4-2D+6E — 307,29F — 1 33 , 656 ' 
+ 25 = 184,72 A — i 53 , 88 B + 181, 00C — 462,51 D 

— 307,29 F + 224868 F + 16694,1 G, 

— 3 = — i 9 , 85 A — 20,69 B + 108,40 C — 60,960 

— 1 33 , 65 E -f- 16694,1 F + 8752,396. 
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Nous en déduisons, par élimination, 

A = — 0,225, 

B = -t- o,344> 

C = — 0,088, 

D = — 0,171, 
fc, — o,323, 

F = -4- 0,000215915, 

G = — o,oo54'j46 20 > 

et les erreurs les plus probables sont données par les for- 
mules 

(o) = — C-t-4»3iF -+-4»3 *g, 

( 1 ) = — B — 24 ,i 6G, 

(2)=— A-4-B + C — 1 53,88 F -+- 19,850, 


d’où l’on déduit les valeurs numériques suivantes : 


(o) = -+■ o",o 65 
( 1 ) = — 0,212 

(2) = -H o, 33 g 

(3) = — o,i 9 3 

( 4 ) “ -F o , 233 

(5) = — 0,071 

(6) ==— 0,162 

(7) =— 0,481 

(8) = + 0,406 


(9) = + 0 ", 02 I 

(10) = + o,o 54 
(n) = — 0,219 

( 12 ) = -4- o,5oi 

(13) = — 0,282 

( 1 4) = — o,256 

(1 5 ) = -F 0,164 

(16) = -+- 0,230 

(17) = — o,i3g. 


I a somme des carrés de ces erreurs est égale à 1 ,2288 ; l’er- 
reur moyenne résultant des 18 directions observées est, par 
conséquent , 


/ 1 ,2288 „ / 

W — = O ,4190. 


25 . 


Afin de donner un exemple de la dernière partie de notre 
théorie, cherchons dans quelle précision les observations 
compensées déterminent le côté Falkenberg-Breilhorn , 


* 


Digitized by Google 


( 1 1 1 ) 

au moyen du côté Wilsede- Wulfsode. La fonction m, par 
laquelle il est exprimé , est, dans ce cas, 

„ m , sin (o( l3 > — *('*) — o", 652 ). sin(w(") — p{'*> — o", 8 i 4 ) 

77 »94 X sin ^ , j — — o",652). sin (p<* > — 1>(*> — o", 8 i 4 )’ 

sa valeur déduite des observations corrigées est 

26786’", 68. 

La différentiation de cette équation fournit, en exprimant 
de (0) , </e (1) ,.elc., en secondes, 

du = [o m , 16991 — rfcO] + o ,n ,o 8836 

— [o m , 03899 ( ) — </p( 13 )] -t- o m , 16731 ( dvi “ ) — </p(’*>) ; 


on déduit de là : 


(al)=- 

o,o 8836 , 

( ht) = -+- 

0 ,.13092, 

(cl) == — 

0,00260, 

(dt)= + 

0,07895 , 

(rl)= + 

0 ,o 38 g 9 , 

(/')=- 

4 o, i 3 i 5 , 

uo = -+• 

10,9957 , 

( // ) = 4- 

0, i 3238 , 


Les méthodes indiquées plus haut donnent, en prenant le 
mètre pour unité de longueur, 

^=0,08329 ou P =12,006. 

On en conclut que l’erreur moyenne à craindre dans la va- 
leur du côté Falkenberg-Breithorn est o m ,2886 ni (m dé- 
signant l’erreur moyenne à craindre dans les directions 
observées, cette erreur étant exprimée en secondes), et, 
par conséquent, si nous adoptons la valeur de m annoncée 
plus haut, cette erreur moyenne à craindre est o m , 1 20g. 

Au reste, l’inspection du système de triangles montre 
immédiatement qu’on pouvait complètement laisser de côté 
la station Hauselberg, sans rompre le réseau qui réunit les 
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quatre autres. Mais il ne serait pas permis pour cela de 
supprimer les opérations qui se rapportent à ce point, car 
elles contribuent certainement à augmenter la précision de 
l’ensemble. Pour montrer plus clairement quel accroisse- 
ment de précision en résulte, nous terminerons en faisant 
de nouveau le calcul, après avoir exclu tous les résultats qui 
se rapportent au poiut Hauselberg. Des dix-buit directions 
mentionnées plus haut, huit cessent alors de servir, et les 
erreurs les plus plausibles sur celles qui restent, sont 


(o) = + o",327 
( i ) = — o , 206 

(3) = — o, 121 

(4) = H-o, 121 
( 6 ) = — o , 1 2 1 


(12) = + o",2o6 

(1 3 ) = — o, 206 

(» 4 ) = + », 327 

( 1 5 ) = + o , 206 
( 1 6) = -t- o , 1 2 1 


La valeur du côté Falkenberg-Breithorn devient alors 
26766"', 63, résultat peu différent de celui qui a été obtenu 
plus haut. .Mais le calcul du poids donne 

^=o,i3o82 ou P ■= 7 , 644 > 


et l’erreur moyenne à craindre est , en mètres , 
/*- 

0,36169m — o m ,i5i5. 


On voit que par l’adjonction des opérations qui se rap- 
portent à Hauselberg, le poids de la détermination du 
côté Falkenberg-Breithorn est augmenté dans le rapport 
de 7,644 à 12,006, c’est-à-dire dans le rapport de l’unité 
à 1,571. 
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NOTES 


NOTE I. 

Exposition de la méthode des moindres carrés. 

. (Extrait du Theoria Motus Corporum eatlesiium.) 

1 . . . 

Abordons maintenant une recherche beaucoup 

plus générale et des plus fécondes dans toute application 
du calcul aux phénomènes naturels. Soient V, \', etc., 
U fonctions des v inconnues p , q, r, s, etc. , et supposons 
que des observations directes aient donné, pour ces fonc- 
tions , les valeurs 

V= M, V' — M' j V" = M", . . i . 

En général , le calcul de ces inconnues constituera un pro- 
blème indéterminé, déterminé ou plus que déterminé, sui- 
vant que l’on aura 

f* O) !>■=*, ou (*). 

Nous ne nous occuperons ici que du dernier cas, dans lequel 
évidemment il ne serait possible d’obtenir une représenta- 
tion exacte de toutes les observations, que si ces observa- 


(*) Si, dans ce troisième cas, p-\- i — y des fonctions V, V', V", etc., 
pouvaient être regardées comme des fonctions de toutes les autres , le pro- 
blème deviendrait plus que déterminé relativement à ces fonctions , mais 
indéterminé relativement a p, q, r, s, etc. On ne pourrait pas en déduire 
les valeurs de ces dernières , même si les valeurs des fonctions V, V', V", etc., 
étaient d’une exactitude absolue : mais nous excluons ce cas particulier de 
nos recherches. 

G 


8 
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lions n’étaient affectées d’aucune erreur. Mais comme cela 
n’a jamais lieu dans la nature, on devra regarder comme 
possible tout système de valeurs des inconnues p, <7, r, 
s, etc., desquelles résultent, pour les fonctions V — M, 

X' — M', X" — M", des valeurs qui ne surpassent pas les 
limites des erreurs que l’on peut commettre dans les obser- 
vations, mais on ne doit pas regarder tous ces . systèmes 
possibles comme jouissant du même degré de probabilité. 

Supposons d’abord, dans toutes les observations , un étal 
de choses tel, qu’il n’y ait pas lieu de regarder l’une d’elles 
comme plus exacte qu’une autre, c’est-à-dire, que l’on 
doive regarder des erreurs égales dans chacune d’elles 
commeégalcmenl probables. La probabilité qu’une erreur A 
soit commise dans l’une des observations sera une fonction 
de A, que nous nommerons »(A). Quoique cette fonction ne 
puisse être assignée d’une manière précise, on peut du 
moins affirmer qu elle doit devenir maximum pour A = o, 
avoir dans la plupart des cas la même valeur pour des va- 
leurs de A égales et de signes contraires, et, enfin, s’éva- 
nouir quand on donne à A une valeur égale ou supérieure 
à l’erreur maximum ; çp(A) doit donc, à proprement parler, 
être rapportée à la classe des fonctions discontinues, et, si 
nous nous permettons, pour la facilité du calcul , d'y sub- 
stituer une fonction analytique, il faudra que cette dernière 
soit choisie de telle sorte qu’elle tende rapidement vers o à 
partir de deux valeurs de A, l’une supérieure, l'autre infé- 
rieure à o, et qu'en dehors de ces deux limites on puisse la 
considérer comme nulle. Or la probabilité que l’erreur soit ; 
comprise entre A et une quantité A -t- r/A qui en diffère 
infiniment peu, sera exprimée par y (A). r/A, et, par suite, 
la probabilité que l’erreur est comprise entre D et D\ pat 



Cette intégrale, prise depuis la plus grande valeur néga- 
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tive de A jusqu’à sa plus grande valeur positive, ou plus 
généralement depuis A = — oo jusqu’à A = oo , devra né- 
cessairement être égale à i. On aura donc 

y (a) c/A = I . 




Supposons donc qu’on ait un système déterminé de valeurs 
des quantités p, g, r, s, etc. : la probabilité que l’obser- 
vation donnera pour V la valeur M, sera exprimée par 
f (M — V) , après qu’on aura substitué dans V les valeurs 
de g, r,s, etc.; de même <f (M' — \"), <f (M" — Y ") , etc., 
exprimeront les probabilités pour que les observations dou- 
aient aux fonctions V', V", etc., les valeurs M', M", elc. 
C est pourquoi , tant qu’on pourra considérer toutes les ob- 
servations comme des événements indépendants les uns 
des autres , le produit 

* (M — V) y (M' — V') y (M" — V").. . . — 12 , 

exprimera la probabilité que toutes' ces valeurs résulteront 
en même temps des observations. 

2 . 

De même qu’en se donnant des valeurs quelconques des 
inconnues, il en résulte, avant toute observation , une pro- 
babilité déterminée pour un système de valeurs des fouc- 
tion V, V', Y", etc., de même , après que l’observation aura 
donné pour ces fonctions des valeurs déterminées, il en 
résultera pour chaque système de valeurs des inconnues 
qui en découleront , une probabilité déterminée : car il est 
clair qu’on devra considérer comme les plus probables les 
systèmes qui donnent à l’événement observé la plus grande 
probabilité. L’appréciation de cette probabilité peut s’ob- 
tenir par le théorème suivant : 

Si, en adoptant une certaine hypothèse H, la probabilité ■ 
d’ün événement déterminé E est h, mais gu en adoptant 
une autre hypothèse H', exclusive de la première et ayani 

8 . 
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à priori in même probabilité, la probabilité du même évé- 
nement soit h'.: je dis que lorsque /’ événement E aura eu 
lieu , la probabilité que H soit la vraie hypothèse sera à la 
probabilité que H' soit la vraie hypothèse comme h est à h 1 . 

Pour le démontrer et afin de distinguer toutes les circon- 
stances d’où peut dépendre, soit que l'hypothèse H ou H', 
ou toute autre ait lieu, l’arrivée d’un événement E ou d’un 
autre événement, formons un système des cas différents qui 
peuvent se présenter et que nous regarderons comme éga- 
lement probables à priori (c’est-à-dire tant qu’il y a doute 
si c’est l’événement E ou un autre qui aura lieu). Ces cas 
peuvent être ainsi distribués : 


NOMBRE DES CAS. 

HYPOTHESE 

propre à ces cas. 

ÉVÉNEMENT 

qui (loi! en résulter. 

m 

H 

E 

ri 

fl 

Différent de E 

m' 

H' 

E 

//' 

H' 

Différent de E 

* m" 

Différente de 11 et de H' 

E 

n" 

Différente de H et de H' 

Différent de E 


On aura d'après cela 


ni n 


h' — 


Or, avant l’arrivée de l'événement, la probabilité de l’hypo- 
thèse H était 

ni -t- n 

ni -+- n -f- ni' -y- n' + ni" n " 


Après l'événement qui exclut n n' -f* n" cas, parmi ceux 
qui sont possibles, cette probabilité sera 

ni ‘ . 

m -i- m' -4 -ni' , ' . , 
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De même les probabilités de l'hypothèse II' avant et après 
l’événement sont respectivement 

m' -f- r>' m' 

m 4 - n 4- m' 4- «' - 4 - m" 4- n" m 4- /«' 4- ni " 1 

mais, comme on a supposé que les hypothèses H et H' avaient 
avant l’événement la même probabilité, ou aura 

m 4- n — ni' 4- 


d’où résulte immédiatement la vérité du théorème. 

Si maintenant on suppose qu'on n’a, pour déterminer les 
inconnues, que les observations 


V = M, V' = M', V" = M", . . , 

et que tous les systèmes de valeurs des inconnues étaient éga- 
lement probables avant ces observations, il est visible que 
la probabilité d'un certain système, après ces observations, 
sera proportionnelle à SI. C’est-à-dire que ISldpdqdr. . . 
exprimera la probabilité que les valeurs des inconnues 
soient respectivement comprises dans les limites infini- 
ment voisines p et p 4- dp, q et q 4- dq , r et r -j- dr, , . . , 
X représentant une- quantité indépendante de p, q , r , 
r, etc.; et l’on aura évidemment 

, p X) p CO p » 



De là résulte naturellement que le système le plus pro- 
bable des valeurs de p, q , r, etc., correspondra au maxi- 
mum de SI, et se tirera des v équations 


(l il 

si lou pose 
V— M = (’, V'— M': 


du (la 




i± 1*1 

f (A) . d a 
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ces équations prendront la forme suivante : 
dv, dd dv" 

_ T(l ,) + _; ? + t' K’’- 


^ ’t \ 

T q i W' 


dd 

d'I 




dv" 

dT, 




■ o , 


De là résulte qu’on pourra obtenir par l’élimination une 
solution pleinement déterminée du problème, dès que la 
nature de la fonction ç' sera connue. Mais comme cette 
lonction ne peut être définie à priori , abordons la question 
à un autre point de vue et cherchons une fonction acceptée 
tacitement comme base, en vertu d’un principe simple et 
généralement admis. Or on a coutume de regarder comme 
un axiome l’hypothèse que si une quantité a été obtenue par 
plusieurs observations immédiates, faites avec le même soin 
dans des circonstances semblables, la moyenne arithmétique 
des valeurs observées sera la valeur la plus probable de celte 
quantité , sinon en toute rigueur , du moins avec une grande * 
approximation, de telle sorte que le plus sûr soit toujours 
de s’y arrêter. Si donc l’on pose 


et 


V = V'= \"...=p, . 
M M' + M" . . 


on devra avoir en général 

<p' (M — p ) -+- f' ( M' — ? p) -H (M" — p) -+■ . . , — o 

pour toute valeur entière et positive de p. Faisant en- 
sui te 

M — pN, 


on aura généralement 


?'[({* — ') N J = (! — {'•)?'(— N), 


Digifeed by Google 



( "9 ) 

d oit l'on tire facilement que «loi t être en général une 

A 

constante A. On aura donc 

lng<p (A) = ^ AA- 4- const. = y A A ! -I- log y . , 

d’où 


Or on voit facilement que la constante A doit être négative, 
pour que fi puisse devenir maximum : posons donc 



et comme, d'après un élégant théorème de Laplace, on a 




notre fonction deviendra 


? (A) — 




i. 


La fonction que nous venons de trouver ne peut pas ex- 
primer, en toute rigueur, la probabilité des erreurs , puis- 
que les erreurs possibles étant toujours renfermées entre 
certaines limites, la probabilité d’erreurs plus grandes de- 
vrait être toujours nulle, tandis que notre fonction a toujours 
une valeur linie. Cependant ce defaut, que présenterait 
également' toute autre fonction analytique, n’a aucune im- 
portance dans les applications, parce que la valeur de notre 
fonction décroît si rapidement, pour peu que AA ait une 
valeur considérable, qu’on peut, en toute sûreté, la re- 
garder alors comme équivalente à o. D’ailleurs, la nature de 
la question ne permettra jamais d’assigner les limites des 
er reurs avec une rigueur absolue-. 

Au reste., la constante h peut être regardée comme scr- 



( iao ) 

vant de mesure à la précision des observations. Si en effet la 
probabilité de l'erreur A dans un système d’observations 
est exprimée par 

h — h ■ a 5 

~ ~r e 

\n 

et dans un autre système d’observations plus ou moins 
exactes que les premières par 

h' —h” A’ 

r ’ - 

la probabilité que dans une observation du premier sys- 
tème l’erreur soit comprise entre les limites ■. — d et -f- à, 
sera exprimée par * 

r A.-‘ ,a x 

J— a - \ K 

et de même la probabilité que l’erreur d’une observation du 
second système soit comprise entre les limites — d' et -f-d', 
sera exprimée par 

f 3 ' 

J -S’ V’f 

or ces intégrales sont manifestement égales lorsqu’on a 

i,s = w y. 

Si, par exemple, on a ' . 

fi' = 2 fl, 

une erreur double dans le premier système sera commise 
aussi facilement qu’une erreur simple dans le second , de 
sorte que les dernières observations, pour nous servir 
d’une expression consacrée par l’usage, jouissent d'un degré 
de précision deux fois plus grand. 

Voici maintenant quelques conséquences de cette loi. H 
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ésl clair qu'il faut , pour que le produit 

devienne maximum , (pie la somme 
<>’ v 7 -+- v " 7 


devienne minimum. Donc le système de valeurs des incon- 
nues p, fj , r, s , etc., le plus probable correspond au cas où 
les carrés des différences entre les valeurs observées et les 
valeurs calculées des fonctions V, V*, etc., donnent la 
somme la plus petite possible, pourvu que toutes les obser- 
vations soient également présumées précises. 

. Ce principe, qui est de la plus grande utilité dans toutes 
les applications des mathématiques à la philosophie natu- 
relle, doit être regardé comme un axiome, au même titre 
que le principe qui nous fait adopter la moyenne arithmé- 
tique des valeurs observées d’une même quantité comme la 
valeur la plus probable de cette quantité. 

Le principe s’étend sans peine au cas d’observations 
d’une précision inégale . Car si les précisions des observa- 
tions par lesquelles on a trouvé 

V = M, V' = M', V" = M" , . . . , 

sont représentées respectivement par A, h', h ", etc., c’est- 
à-dire si l’on suppose que des erreurs réciproquement pro- 
portionnelles à ces quantités puissent être commises avec la 
même facilité, il est clair que cela revient au même que si, 
pardesobservalions d’une égale précision (représentée pari), 
les valeurs des fonctions AV, A'\ ', A"V",etr., avaient été 
trouvées égales à AM, h'M.\ w h" M". etc. ^ c’est pourquoi le 
système le plus probable de valeurs des 'quantités p, r , 

s, etc., sera celui où la somme 

h 7 v 7 -t- h ' 7 v ' 7 -t- A"* v " 7 -t- ■ . . , 

* * • \ -* 
c est-à-dirc ou la somme des carrés des différences entre 
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les valeurs observées et calculées , multipliés respectivement 
parles carrés des nombres qui expriment le degré deprê- 
xision , devient un minimum. Par là , il n’est pas même né- 
cessaire que les fonctions V, V', V", etc., se rapportent à 
des quantités homogènes, mais elles pourront représenter 
des quantités hétérogènes (par exemple des secondes d’arcs 
et de temps) ; pourvu que l’on puisse estimer le rapport 
des erreurs qui, dans chacune de ces grandeurs , peuvent 
être commises avec la même facilité. 

6 . 

Le principe exposé dans l’article précédent se recom- 
mande aussi par cela qu’il réduit le calcul numérique des 
inconnues à un algorithme très-expéditif, quand les fonc- 
tions V, V', V", etc., sont linéaires. Supposons 

V — M = r = — m -+- a p b } + c r + d » + 

V' — M' = «>' — — m' a' p -f- b' q -+r d r -j- il' s -(-... , 

V" — M” — v" = — m" -f- a" p -f- b" q -f- c" r 4- d"s 4- . . . , 


Posons 

av 4- a! v' 4- 4- • . . — l 5 , 

b v 4- b' d 4- b" v" 4- . . . = Q , 
c r 4- d d 4- c i> 4- • * . — R , 
c/i' 4- d'd 4- c/V'4-. . S, 


alors les v équations de l’art. 3, qui déterminent les valeurs 
des inconnues, seront 

P = o, Q = o, R = o, S — o , ... , 

si nous supposons les observations également bonnes, cas 
auquel nous pouvons ramener tous les autres , cdmme nous 
l’avons montré dans l’article précédent. On a ainsi autant 
► d’équations linéaires que d’inconnues : on les résoudra par 
la méthode ordinaire. 
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\ oyons maintenant si cette élimination est toujours pos- 
sible ou si elle peut donner une valeur indéterminée ou im- 
possible. Il résulte de la théorie de l'élimination que le se- 
cond ou le troisième cas aura lieu si , en laissant de côté une 
des équations 

P = o, Q = o, R = u , . . . , 

on peut déduire des équations conservées une équation 
identique ou contradictoire à celle que l’on a omise , ou, ce 
qui revient au môme, si l’on peut assigner une fonction li- 
néaire 

«P-f PQ+7R-+- 

qui soit identiquement nulle ou qui ne contienne aucune 
des inconnues. Supposons donc que l’on ait 

« P 4- PQ 4- 7 R 4- . . . = x , 

on a l’équation identique 

( t> 4- m ) v H- ( / 4- m' ) v' 4- ( v" -j- m" ) v" -+- . . . . 

= p P -4- q Q 4- rR -t- s S -)-... . 

Si l'on suppose qu’en faisant 


p=.-xx, q=$x, r = qx,.. , 

les fonctions v, e', v", etc., deviennent respectivement 
— ni -f \x, — m' + X'x, — rn” 4- X".r, . . , 

on aura l’équation identique 

( X* 4- X 3 4- X 3 4- • . ) x 7 — ( X /// 4 - X 3 ft\' 4- X" ni 4 - • • ■ ) x — r.x^ 
et, par suite, 


X 3 4- X' 3 4- X"* 4- . . . 


o, y 4 - \m 4 - X 3 rn f 4- X /3 ni" 4- ... — O , 


d’où résulte 
et , par suite , 



* = 0; 
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c'est-à-dire que les fonctions V, V', V ", etc., devraient ne 
pas changer si p , q, r, s, etc., reçoivent des accroisse- 
ments quelconques proportionnels aux nombres», [3 , y,’ 

o, etc. Un pareil cas, dans lequel la détermination des in- 
connues ne serait pas possible, même si l’on donnait les 
vraies valeurs des fonctions , n’appartient pas à notre sujet ,, 
comme nous en avons averli plus haut. 

Au reste, on peut réduire facilement tous les cas à celui 
où les fonctions V, V', etc., sont linéaires. Désignons 
par ît , x, p, a, etc., des valeurs approchées des inconnues 

p, q, r, s, etc. (que nous obtiendrons en faisant usage de v 
équations prises parmi les p. équations 

V=M, V'=M', V" = M" 

et posons 

p = n+p', q = y -1- q', r = pP-r', S—o - (-/, 

il est clair que ces nouvelles inconnues seront si petites, que 
leurs carrés et leurs produits seront négligeables, et que les 
équations deviendront linéaires par suite des substitutions 
indiquées. Que si, à la fin du calcul, on trouve contre 
toute attente que les valeurs de p', q\ r 1 , s', etc., qu’on 
en tire soient trop considérables, et qu’il paraisse peu sûr 
de négliger leurs carrés et leurs produits, on remédiera 
à cet inconvénient en répétant la même opération (mais 
en prenant pour jt, y , p, !7,*etc., les valeurs corrigées de 
Pi < 7 , 'VU etc.). 

7 . 

Tant qu’on n’a qu’une seule inconnue p, pour la déter- 
mination de laquelle on a trouvé que les fonctions 

n p -t- n, a' p 4- ri', a" p -+- n", . . . , 
prenaient respectivement les valeurs 
. . M, M\ M", . . , 

et cela par des observations également exactes , l« valeur la 
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plus probable de p est 


a m + n' ni’ •+• a" m" -f- . . . 
a' -+- a " 1 -f- a”’ -h . . . * 

en posant 

m — M — », m'rir — m" — M" — n" , ... 

Pour apprécier le degré de précision qu’on doit attri- 
buer à cette valeur, supposons que la probabilité d’une 
erreur A, commise dans les observations, soit exprimée 
par 

h - h' i 1 

~T e ■ ; 

V* 

il en résultera que la probabilité que la vraie valeur de p 
soit A -f -p', sera proportionnelle à la fonction 

g — h' [{ap — ni)’ -¥-{a' p — m')’ (a" p — m")* 

( 

dans laquelle on aura fait 

p = A+y. 

L’exposant de cette fonction peut être réduit à la forme 

— h' (a 1 -+- a '* + a"' h )!p l — apA-t-B), 

'dans laquelle B désigne une quantité indépendante de p : la 
fonction sera par suite proportionnelle à 

e — hHa' 4- a-' -4-o’S-e '..)/>'* ■ 

On voit que le degré de précision qu’il faut attribuer à la 
valeur de A est le même que si cette valeur avait été 
trouvée par une observation immédiate dont la précision 
serait à la précision des observations primitives comme 


h y 'ai -+• n' ! -H n" 1 est à h, 


ou comme 


y'»’ -t- ri' n”' H- . . est a i . 
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8 . 


Avautde rechercher, dans le cas de plusieurs inconnues , 
le degré de précision qu’on doit attribuer à chacune d'elles, 
il importe d’étudier plus attentivement la fonction 


■ (i ' -4- V - 


que nous désignerons par W.. 

I. Posons 

1 rfW . , „ . 

— =/> = A 4- a/> + 8 q +7r + 5.(+ . . . 

2 dp 


et 


AV — — = AV', 


il en résulte évidemment 


P; 


comme on a 


d AV' d AV ip' dp' 


dp dp a dp 

on voit que la fonction W'sera indépendante de p. Le coeffi- 
cient 

a = a 1 + a' 1 - 1- fl" 1 -I- . . . , 
sera toujours évidemment une quantité positive. 

II. De meme, posons, * * 

i dW 


7. dq 


et 


on aura 


et 


: q' = 1' + + 7’r + I'j+.. 

o ' 2 

• AV'-?-=AV'', 

P 

, _ i jw // Æ— o_ P 

• ^ 1 dq a. dq a 


d AV" 
dq 


— O. 


Donc la fonction W" est indépendante à la fois de p et - 
de q. Ces circonstances n’auraient plus lieu si l’on pouvait 
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P' = O. 

Mais il est clair que W' se déduit de 


en remplaçant dans t', e', v'\ etc., la quantité p par sa va- 
leur tirée de l'équation 

/>' = o : 

donc J3' sera la somme des coefficients de q 5 dans (<*, v'*, 
e"’, etc., après cette substitution. Mais ces coefficients sont 
tous des carrés et ne peuvent s’évanouir tous à la fois, si ce 
n’est dans le cas, que nous excluons de nos recherches, où 
les inconnues seraient indéterminées; donc /S' doit être po- 
sitive. 

III.- Si l’on pose, enfin. 


1 d AV" 

2 dr 


= r'= r + y "r + 3"s- h. 


et 


on aura 


AV" 


— — W" 
« — ’* < 


* * 

et W'" sera indépendant de p, de q et de r. On prouvera 
comme plus haut que le coefficient y" doit être positif. O 11 
voit, en effet, facilement que y" es't la somme des coeffi- 
cients de r 2 dans e ! , e'*, e" ! , etc., après que les quantités p 
et q ont été éliminées de e, e", etc., .à l’aide des équa-* 
lions , 

P'- = 0 , q' = o. 

IV. De là même manière en posant 


1 d AV" , . 

- ; = J = À -+- l) J -t- . 

2 ds 


AV* — AV'" — — 
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on aura 


, - S , 3' , 3 " , 

' =S-- P ~ f <, - ? T>, 

W ,T sera indépendant de p, r et J, et à'" une quantité 
positive. . 

V. S’il y a un plus grand nombre d’inconnues, on conti- 
nuera delà même manière et l’on aura enfin 


w = - P ’ -, q ! -+- — r'»- +pt' , +...+ const. , 

expression où a, (3', y", d'", etc., désignent des quantités 
positives. 

\ I. On a déjà vu que la probabilité d'un système de va- 
leurs de p , q, r, s , etc., était proportionnelle à la fonction 
é~ h " : par conséquent, la valeur de p restant indétermi- 
née, la probabilité d’un certain système de valeurs de q, 
r, 5, etc., sera proportionnelle à l’intégrale 


£ 


— A* VV , 

e rfp, 


qui est égale, d’après le théorème de Laplacè, à 

h cl 2 n~ c ' 

et cette probabilité sera proportionnelle à la fonction 
• c- A ’ w '. 

De même, si l’on considère de plus q comme indéterminé, 
la probabilité d’un système de valeurs de r, s, etc., sera 
proportionnelle à 


£ 


IrW 


r/f/, 


c’est-à-dire à 

/» K' J n-'-c ’ > 
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et, par suite, proportionnelle à c~~ h '^ v '. De même si r 
est aussi regardée comme indéterminée, la probabilité 
d’un système de valeurs déterminées de s, etc., sera pro- 
portionnelle à é~ h ‘ w , et ainsi de suite. Supposons que 
le nombre des inconnues se réduise à quatre; les conclu- 
sions seraient les mêmes dans le cas général. La valeur la 
plus probable de s sera 

_ — 

— 

et la probabilité qu’elle différera de a de la véritable va- 
leur sera proportionnelle à 

— hh n ' »* . 

C i 

d’où nous concluons que \Jd'" mesure la précision relative 
à cette détermination , en prenant pour unité la précision 
des observations primitives. 


. 9 . 

Par la méthode du paragraphe précédent un certain 
degré de précision a été assigné à la seule inconnue qui, 
dans le travail de l’élimination, a été gardée la dernière. 
Pour éviter cet inconvénient, nous allons calculer à'" d’une 
autre manière. 

Des équations 


P =P', 
Q= 




9 
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on tire, en les résolvant par rapport à p', q', r\ s', 
p' = P, 

<]' = Q + ^ P > 

/ = R -H *itë/ Q -f“ «A> # P > 

i' = S+ £"R -+- ifl" Q ■+■ X"P, 


de sorte que X, X 1 , X", sont des quantités 

déterminées. On aura donc (en restreignant à quatre le 
nombre des inconnues) 


t 


r 

r 


Jl+jQ + yl + ^S, 


d’où résulte la conséquence suivante : Les valeurs des in- 
connues p, <7, r, s, etc., que l’on doit tirer des équations 

P — o , Q— -o, R = o, S = o, 

sont évidemment exprimées par des fonctions linéaires de 
P , Q , R , S , etc-. , savoir : 

p m: L -H A P -H B Q -4- C R-t-D S H — . - * , 
y = L' +A' P+B' Q+CR + D' S + ,.., # 

r = L' + A"P + B"Q-t-C"R + D"S + ..., 
s — L'" -4- A'" P+ B"Q + CR+ D"'S 


Cela posé, les valeurs les plus probables de ces incon- 
nues sont respectivement L, L', L", etc. Les degrés de 
précision qui doivent être attribués à ces déterminations 
sont respectivement 

î 1 1 i_ 

Ta’ v'b 7 ’ Vc 5 ’ ' VU* 

en prenant pour unité la précision des observations pri- 
mitives -, car ce que nous avons dit plus liaut de l’incon- 
nue s |^pour laquelle D'" répond s’applique aux au- 
tres inconnues par une simple permutation. 
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10 . 

Pour éclaircir par un exemple les- recherches qui pré- 
cèdent, supposons que par des observations pour les- 
quelles une égale précision doit être présumée, ou ait 
trouve 

p — q 4- 2 r = 3, 

3/j h- 2 q — 5r = 5, 

4 /■> 4- q -i- 4 r = 2 1 , 


mais que, par une observation à laquelle une précision 

égale à - doit être attribuée, on ait trouvé 
° 2 ’ _ < 

— IJ) ■+■ 6q -|- 6 r = 28 . 

A cette dernière nous substituerons la suivante , 


— p + - 3 q + 'ir= 14 , 

que uous supposerons provenir d’une observation aussi 
précise que les premières. De là on tire 


P = 27 p +6 q — 88, 

Q — 6 p 4- i 5 y 4 - r— 70, 

R = 7 4- 54 r — 107, 

et par l’élimination 

19899/; = 49*54 4-809P— 324 Q 4- 6R, 

737 q — 2617— t2P4- 54 Q — R, 

6633 r — 12707 4- 2P — 9Q 4- i23R. 

Les valeurs les plus probables des inconnues seront donc 

p=z 2,470 q = 3,55i, r= 1,916, 

avec des degrés de précision égaux respectivement à 


V- 


96*19 . 


809 


: 4-9*3. 


i/ih- 

V 54 - 


3 , 69 , 
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11 . 

Le sujet que nous avons traité jusqu’ici donnerait lieu 
à d’élégantes recherches, auxquelles nous ne nous arrête- 
rons pas, pour ne pas trop nous écarter de notre objet 
principal. Par la même raison, nous réservons pour une 
autre occasion l’exposition des artifices qui permettent de 
réduire le calcul à un algorithme plus expéditif. Qu’on 
nous permette seulement d’ajouter une seule observation. 

Lorsque le nombre des fonctions ou des équations pro- 
posées est considérable, le calcul est surtout rendu pénible 
par cette circonstance, que les coefficients par lesquels on 
doit multiplier les équations primitives pour obtenir P, Q, 
R , S , etc. , sont presque toujours des fractions décimales 
compliquées. Si l’on ne croit pas important dans ce cas de 
calculer ces produits avec le plus grand soin à l’aide des 
Tables de logarithmes, il suffira le plus souvent de leur 
substituer des nombres plus simples qui en diffèrent peu. 
Il ne peut en résulter d’erreurs notables qu’autant que la 
précision des inconnues devient moindre que la précision 
des observations primitives. 

12 . 

Au reste , le principe d’après lequel la somme des carrés 
des différences entre les quantités observées et les quan- 
tités calculées doit être un minimum, peut encore s’éta- 
blir sans recourir au calcul des probabilités, comme il suit. 

Lorsque le nombre des inconnues est égal au nombre 
des observations, on peut déterminer les premières de 
manière qu’elles satisfassent aux secondes. Mais lorsque 
le premier nombre est le plus petit des deux, on ne peut 
obtenir un accord absolu lorsque les observations ne sont 
pas douées d’une précision absolue. Il faut donc dans ce 
cas chercher à établir l’accord le plus satisfaisant, c’est-à- 
dire à faire en sorte que les différences soient atténuées 
le plus possible. Mais cette idée a par elle-même quelque 
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chose de vague. Eu effet, quoiqu’un système de valeurs des 
inconnues doive être sans aucun doute préféré à un autre 
système où toutes ces différences seraient respectivement 
plus grandes, le choix entre deux systèmes dans l'un des- 
quels l’accord serait plus satisfaisant pour quelques-unes 
des observations, mais moins satisfaisant pour d’autres, est 
en quelque sorte arbitraire, et l'on peut évidemment pro- 
poser plusieurs principes par lesquels la première condi- 
tion soit remplie. En désignant par A, A', A", etc., les 
différences entre le calcul et les observations , on satisfera à 
cette condition, non-seulement si A* 4- A' 3 4- A' /3 -4 . . 
devient un minimum (ce qui est notre principe), mais en- 
core si A‘ 4 - A'* 4 - A"* 4 -. • -, ou A 6 4- A' 6 4 - A"« -f- . . . , 
ou généralement une somme de puissances paires, devient 
un minimum. Mais de tous ces principes le nôtre est le plus 
simple, tous les autres nous entraînant dans des calculs 
extrêmement compliqués. Au reste, ce principe, dont nous 
avons fait usage dès l’année 179s, a été donné dernière- 
ment par Legendre dans ses Nouvelles méthodes pour la 
détermination des orbites des comètes, Paris, r 806 ; on 
trouvera dans cet ouvrage plusieurs conséquences que le 
désir d’abréger nous a fait omettre. 

Si l’exposant de la puissance paire dont nous venons de 
parler était infini, nous serions ramené au système dans 
lequel les plus grandes erreurs sont moindres que dans tout 
autre système. 

Laplacc se sert, pour la résolution d’équations linéaires 
en nombre plus grand que les inconnues, d’un autre prin- 
cipe, proposé d’abord par Boscovich, savoir, que la somme 
des valeurs absolues des différences devienne minimum. On 
peut facilement démontrer que le système des valeurs des 
inconnues trouvé par ce seul principe doit nécessaire- 
ment (*) satisfaire à autant d’équations , prises parmi les 


(") Excepté quelques cas spéciaux où il y a indétermination- 
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proposées, qu’ il y a d’inconnues , de sorte <juo les autres 
équations ne sont employées que pour décider le choix que 
l’on doit faire. 

Si , par exemple, l’équation V = M est du nombre de 
celles qui ne sont pas satisfaites, le système des valeurs 
trouvées par le principe en question ne serait pas altéré si 
au lieu deM on avait observé une autre valeur N telle, que, 
n étant la valeur calculée, les différences M — n et N — n 
fussent de même signe. Au reste, Laplace tempère en 
quelque sorte ce principe en y ajoutant cette nouvelle con- 
dition, que la somme des différences, prise avec leurs 
signes, soit nulle. Il en résulte que le nombre des équations 
satisfaites est moindre d’une unité que le nombre des in- 
connues : mais l’observation que nous venons de faire sub- 
siste encore lorsqu’il n’y a que deux inconnues. 


NOTE IL 

APPLICATION DE LA MÉTHODE DES MOINDRES CARRÉS A LA 
CORRECTION DES ÉLÉMENTS DE LA PLANÈTE PALLAS. 

1 . 

M. Gauss a donné, dans le tome I des Mémoires de 
Gottingue , l’application de sa méthode à la correction des 
éléments de la planète Pallas. L’illustre géomètre ayant 
développé sur cet exemple l’algorithme indiqué plus briè- 
vement dans son grand onvrage Theoria Motus Corporum 
cœlestium (voir la Note précédente) , nous avons cru devoir 
traduire ici celte portion de son Mémoire. La première 
partie exigeant la connaissance approfondie de la théorie 
du mouvement des planètes, nous nous dispenserons de la 
reproduire, et nous prendrons pour point de départ les 
douze équations auxquelles les corrections des six éléments 
de l’orbite doivent satisfaire. 
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En désignant ces corrections par 

AL, d 7 , rfre, dp, dQ, di, 

les équations obtenues par M. Gauss sont les suivantes : 

o = — 1 83 ",g 3 + o , 79363 d L -+- 1 43 ,66 d 7 + o , 3 94 9 3 d re 

-t-o,g5g2odp — o, i 8856 d Q -t- 0,17387 d/'; 
o = — 6", 81 — 0,02668 d L -1-46,71 d 7 -+• o,o ?.658 d jt 

— o, 2 o 858 d<p + o, i 5 g 46 d Q -t- 1 ,25782 d / ; 
o =; — o",o6 -t- o, 5888 o AL -+- 358 , iad 7 -t- o, 26208 du 

— o ,85234 dp -t-o, 14912 d Q -t- o, 17775 d<; 

o = — 3 ",og -f- o , o 1 3 1 8 d £ H- 28 , 3 g d 7 — o , o 1 3 1 8 d re 

— o ,07861 dp -f- 0,91704 ilQ + o, 54365 d / ; 
o = — o",oa -+- 1 , 73436 AL 4- 1846, 17 d 7 — o, 546 o 3 d jt 

— 2,05662 dp — o, i 8833 d Q — 0,1 7445 di ; 

o = — 8", 98 — o, i26o6dZ — 227 ,42 d 7 -t- o , 12606 d re 

— o, 38 g 3 gdf -t- o, 17176 d — 1 , 3544 1 d'* 
o = — 2 ", 3 i -1- o , gg 584 di-f- 1 579,03 d 7 -t- o,o 6456 dre 

-t- 1 ,99545 dp — o, 06040 d Q — o, 3375 o di ; 
o = -t- 2",47 — o, 08089 AL — 67,22d7 -+- 0,0808g dre 

— o,ogg7odp — o, 4635 gd Q -t- 1 ,22803 di; 
o =-f- o",oi -t- o, 653 i 1 di-t- i 32 g,ogd 7 + o, 38 gg 4 dre 

— o ,08439 dp — o,o 43 o 5 d -I- 0,34268 di ; 
o = -+- 38 ", 12 — o ,00218 AL -+- 38,47 d 7 -l- 0,00218 dre 

— o, 187 10 dp -l- o, 473 oi d Q — 1 , 1437 i.df ; 
0 = — 3 17", 73+ 0,69957 di-f- 1719,33 d 7 -t- o, 1291 3 dre 

— 1,38787 dp-t-o, 171 3 odQ — o, o 836 o d i ; 
o =-t- 1 17", 97 — o,oi 3 i 5 dZ, — . 43 , 84 ^ 7 -t-o,oi 3 i 5 dre 

-t- o ,02929 d p-t- 1 ,o 2 > 38 d Q — 0,27187 di. 

D’après la nature des observations qui ont fourni la 
dixième de ces équations, elle inspire trop peu de confiance 
pour qu’on juge utile delà faire intervenir, et c’est d’après 
les onze autres seulement que l’on déterminera les six in- 
connues. 
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Les explications suivantes sont littéralement traduites du 
Mémoire de M. Gauss. 


2 . 

Dans l’impossibilité où nous sommes de satisfaire exac- 
tement aux onze équations proposées , c’est-à-dire d’an- 
nuler tous les seconds membres , nous chercherons à rendre 
la somme de leurs carrés aussi petite que possible. 

On aperçoit facilement que si l’on considère les fonc- 
tions linéaires 

n -+- ap -(- b q cr + Js 4- . . . = tv, 
n' -+■ a' p 4- b' q 4- </ r 4 - d' s 4-. . . = tv' , 
n" 4- a" p 4- b" q 4- c"r-\- d" x 4- . 


les équations qu’il faut résoudre pour rendre 

n = iv 7 4- «''*4- w" 1 . . . 

un minimum, sont 

aw 4- a' 4- a" w" 4- . . . = o , 
bw 4- b' w' 4- b"w" 4-. ..= 0, 
civ 4- e' a/ 4- c" w" 4- . . . = o , 


ou, en posant, pour abréger, 

an 4 - a' n' 4- a" n" 4- . . . = [an ) , 
a ’ 4 -a ' 1 4 -a " 1 4-...= («a), 

ab 4 - a' b' 4 - a" b" 4 - . . . = ( ab ) , 


i’ + i” 4 -b"' 4 - = (bb), 

bc + b’c'-+- è'V'4-. : . . . = (èc) , 


p, q, r, s, etc., devront se déterminer par les équations sui- 
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(an) -+■ (aa)p-\- (ab)q -t- ( ac ) r -t- . . . = o, 
(bn) -h (ab)p + (bb)q -+- ( bc) r + . . . = o , 
( c / j ) -t- (ac)p + (bc) q + (ce) r -H. . . = o, 


L’élimination, très-pénible lorsque le nombre des incon- 
nues est considérable, peut se simplifier notablement de la 
manière suivante. Outre les coefficients (an), (an), etc., 

[dont le nombre est ^ (t’-f- 3 i) , si le nombre des incon- 
nues est t ] , supposons que l’on ait calculé la somme 
«’ + /i' J + n " 2 + . . .= (nn)y 


on voit facilement que l’on a 

Cl = ( nrt ) + 2 (an) p + 2(bn) q + 2 (en) r + .. . 

-)- (aa)p 2 H- 2 (ab)pq 2(ac)pr-f~. . . 

-+• ( bb) q 2 2(bc)qr -+- 2 ( bit ) qs H- . . . 

-f- (cc) r’ -+- 2(crf) rs +. . . ; 

et, en désignant 

(an )- (- (aa)p + (ab) q-h. ■ ■ 

par A, tous les termes de A -r qui contiennent le facteur p , 
r (aa) 1 . 

se trouvent dans l’expression fl, et, par suite, 

A’ 


il — 


(aa) 


est une fonction indépendante de p. C’est pourquoi, en 
posant 

/ \ ( an Y ! \ ' 

(^j= ( '" , ’ ,) ’ 

( an ) (M _ , 

( bn ) \ — = (t>n,i), 


(en) 


(aa) 
(an)_(cn) 
(aa) 


■ (cn,i), 
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on aura 
A* 


11 — (âa) = ( nn > l )+ + z( c 'i,i) r + 2[‘irt,l) s . . . 

-+- (b b , l) y’ -+- 3 ( bc, l'j qr -t- 2 ( bd, 1) qs 
-+- (ce, 1) r’-f- a (cd, 1) rs 


nous désignerons cette fonction par Cl 1 . 

De même, en posant 

(bn,i) + (bb,i)q + (bc,l)r ..= B, 
la différence 

B J 


n'- 


(bb,i) 


sera indépendante de q\ nous la représenterons par fl". 
En posant, de même, 

. (bn,\)(bc,\) 

(bb , 7 ) ~~ 

(cr,,) (M,I) — (cc ’ 2) ’ 


et 

(c«, 2) -t- (cc, 2) r -f- (cd, 2 )s+... = C, 
la différence 

„ C* 

il" 


sera une fonction indépendante de r. 
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En continuant ainsi , nous formerons une suite d’expres- 
sions fi, fl', fl", etc., dont la dernière sera indépendante 
des diverses inconnues, et représentée par (nn, p), si p 
désigne le nombre de ces inconnues ; nous aurons alors 



B’ 


c* 

( ce , 2) 


D J 


• . (n«,p). 


On prouvera facilement que fi étant une somme de carrés 


«>’ -t- IV 1 ■+■ w'" 1 -+- . . . , 

et ne pouvant devenir négative, les diviseurs (au), ( bb , i), 
(ce, a), etc., sont tous positifs. (Nous supprimons, pour 
abréger, le détail de la démonstration.) D’après cela, la va- 
leur minimum de fi correspond évidemment aux valeurs 
des inconnues, pour lesquelles 


A=o, B = o , C = o,..., 

et, en commençant à résoudre le système par la dernière 
équation, qui ne contient qu’une inconnue, on trouvera 
les valeurs de p, <7, r, s, etc., sans avoir aucune élimination 
à effectuer. La méthode donne, en même temps, la valeur 
minimum de fi , qui est (nn , p). 


3 . 

Appliquons ces principes à notre exemple , dans lequel p , 
<7, r, s, etc., sont remplacés par d Z, d 7 , dtt, d®, dQ , dt. 
J’ai trouvé, par des. calculs exécutés avec soin : 


(««) 

= 

148848 , 

(ac) 

— — o>° 9 3 44 > 

(an) 

= — 

3 7 io 9 > 

(ad) 

= — 2, 28516, 

(bn) 

= — 

58 o 1 04 , 

(ac) 

= — 0 , 34664 , 

(en) 

= — 

n 3 , 45 , 

(«/) 

= — 0,18194, 

(dn) 

= ■+■ 

268 , 53 , 

(bb) 

= -+- 10834225 

(en) 

= ■+■ 

94 > 26 > 

(bc) 

= — 49,06, 

(/») 

= — 

3 i ,81 , 

(bd) 

= — 3229,77, 

(an) 

= + 

5,91569, 

(be) 

= -.98,64, 

(ab) 

= + 

7203,91 

W) 

= -.43,05, 
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(cc) = + 0,71917 
(crf) = + 1 , i3382, 
(ce) = + 0,06400, 
(<f) =-}-o, 2634 i, 

(dd) = + 1 2,oo34o , 

( de ) = — 0,37137, 

(df) = — 0,11762, 
(ee) =+2,28215, 
(e/) = — o,36i36, 
(ff) = + 5,62456. 

D’où l’on déduit : 

(nn,i) =+125569, 

(bn,i) =— i38534, 
(c«,i) = — 1 19,31 , 

( dn,i ) = — I25,l8, 
(en,i) = + 72,52, 
(fn,i) = — 43,22, 
(bb,i) = + 2.458225, 
( bc , i ) = + 62,13, 
(bd,i) = — 5io,58, 
(be,i) =+ 2i3,84, 
(b/,i) =+73,45, 
(ce, 1) = + 0,71769, 

( cd , 1) =+1,09773, 
(ce,i) = . — o,o585?., 
(c/, 1) = + o,26o54, 
(dd, 1) =+11,12064, 
(de, 1) =— o,5o528 , 
(rf/,i) = — 0,18790, 

( ce , 1 ) =+2, 2.6l85, 
(c/,t) =—0,37,202, 
(/f, 1) = 4- 5,61905. 

De la même manière : 

( nn , 2) = + 117763, 
(c/1,2) = — 1 i 5 , 8 i , 
(rf/i, 2) = — i 53,95, 
(e/i, 2) = + 84,57, 


( 140 ) 

(/«> 2) = — 39,o3, 
(cc, 2) = + 0,71612, 
(ctf, 2) = + 1,1 io 63 , 
(cc, 2) = — 0,06392, 
(c/, 2) = + o, 25868 , 
(«W, 2) = + 1 1 ,01466, 
(de, 2) = — 0,46088, 
(rf/, 2) = — 0,17265, 
(ee, 2) = + 2,24325, 
(c/, 2) = — 0,37841 , 
(//, 2) = + 5 , 6 i 686 . 

D’où : 

•(/i/i, 3) = + 99034, 
( d />, 3) = + 25,66, 
( en , 3) = + 74,23, 
(/"> 3) = + 2,75, 

(iW, 3) = + 9,29213, 
(rfc,3) = — 0,36175, 
(df , 3) = — o,5 7 384, 
(ce, 3) = + 2,23754, 
(c/,3) = — 0,35632, 
( ff , 3) = + 5,52.342. 

De meme : 

(/!/!, 4) = + 98963, 
(e/i,4) =4- 75,23, 
(/«,4)= + 4,33, 

(cc, 4) = 4- 2,22346, 
(e/,4) = — 0,37766, 
(//, 4) = + 5,48798. 

D’où : 

( nn , 5) = + 96418, 
(//1, 5)= + 17,11, 
(//, 5) =+5,42383. 

D’où enfin : 

( nn , 6) = + 96364. 
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Nous avons donc les six équations suivantes : 

o = -t- 17",! 1 -H 5,42383 <1/, 
o = -I- 75", ?.3 -+- a, 25.346 d Q — o , 37766 d/ , 
o = -i- a5",66 +9,2g2i3dip — o,36i75dQ — 0,57384^/, 
o — 1 1 5",8 1 -4-o,7i6ia<ljr -t-i,no63dy — o,o639adQ 
-+• o , ?,5868 d / , 

0 = — 1 3854 " -t- 2458225 d 7 -+-62,i3djr ■ — o, 5 io 58 dy 

- 4 - 21 3 , 84 d ^ -+- 73,45 d» , 

0 = — 37 1 ",09 + 5,91 569 AL -l- 7203,91 d 7 — o,oo 344 dir 
— 2. .ao 5 i 6 dif — o, 34664 d Q — o, 1 8 1 94 d / ; 

d’où l’on déduit : 

di =— 3", i5 ; 

dR=- 34", 37; 
d<p = — 4"» 2 9» 
d jt = 4- 166", 44 ; 

d7 =■+• o", 054335; 

d L = — 3", 06. 

Telles sont les corrections qu’il faut apporter aux élé- 
ments trouvés d’abord pour la planète. 

NOTE m. 

MÉMOIRE SUR LA DÉTERMINATION DE LA PRÉCISION 
DES OBSERVATIONS, 

PAR U. GAUSS. 

(Extrait du Zeitschrift fur Astronomie und Verwandte Wissenschqften , 
tome I , page 1 85. ) 

i. 

Pour établir les principes de la méthode des moindres 
carrés, nous avons admis que la probabilité d’une erreur 


Digitized by Google 



( *4a ) 

il’observation A soit exprimée par la formule 



où 7t représente la demi-circonférence, e la base des loga- 
ri thmes hyperboliques, et A une constante que l’on peut con- 
sidérer [ Thcoria Motus Corporum cœlestiurn , art. 178 (*)] 
comme la mesure de l’exactitude des observations. Il n’est 
pas nécessaire de connaître la valeur de A pour déterminer, 
à l’aide de la méthode des moindres carrés, les valeurs les 
plus probables des quantités dont les observations dépen- 
dent , le rapport de l’exactitude des résultats à l’exaclitude 
des observations est également indépendant de A. 

Toutefois , comme la connaissance de la quantité A est 
très-intéressante et instructive, je vais montrer comment 
les observations peuvent servir à la déterminer. 


2 . 


Commençons par quelques remarques qui éclairciront 
la question, et représentons par 0 (t) l’intégrale définie 



Quelques valeurs particulières de cette fonction donneront 
une idée de sa marche : 


t — 0,4789363 = p 
t — o,5p5 1161 = px 1 ,247 790 
t = 0,73?. 8691 = p X 1 ,536 618 
t = 0,906 1939 = p x 1 ,900 o 32 
t — 1 =px 2,096716 

/ = 1 , i63 0872 = p X 2,348664 
t = i ,821 3864 — p X 3 , 8 1 8 g3o 
t = 2,327 6754 = p X 4 476 

t = 2,761 o654 = p X 5,768 204 

t = 00 =» 


9 f = o ,5 

et = 0,6 

0t = o,7 

©f = 0,8 

et = 0,8427008 

©/= 0,9 

©t = 0,99 

®'= °»999 
©' = o .9999 
et— 1 


(*) Voyez page il§ de ce volume. 

•• 


J. R. 
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La probabilité que Terreur d’une observation soit com- 
prise clans les limites •+• A et — A; ou, abstraction faite 
des signes, qu’elle ne surpasse pas A, sera égale à 

f f -*" 1 *. 

a 


elle sera le double de l’intégrale en question , lorsque celle- 
ci est prise depuis x = o jusqu’à x = A, et, par consé- 
quent, elle sera égale à 0 (h A). 

Ainsi la probabilité que l’erreur ne soit pas moindre que 

t est ésalc à - , ou é fraie à la probabilité du cas contraire. 

h ° 2 D 1 

J’appellerai donc cette quantité j T erreur probable , et je 
la désignerai par r. 

Au contraire , la probabilité que l’erreur excède 
/’X 2,438664 n’est que de ; la probabilité que Terreur 
surpasse rx 3 , 8 18390 n’est que de et ainsi de suite. 


3. 


Supposons maintenant que les erreurs de m observations 
réellement faites soient «, |3 , 7 , etc., et cherchons quelles 
conséquences on en peut tirer relativement aux valeurs 
de h et r. 

Fin faisant deux hypothèses sur la valeur exacte de h , et 
la supposant égale à H ou égale à H', les probabilités qu’elles 
soient entachées des erreurs «, ( 3 , 7 , etc., seront, pour les 
deux cas, dans le rapport de 


H er HVi X H X H e~ a "*' X .. . 


— H’,?* 


-H’-/ 


H'*- H ' S “* X H'e-H”' 3 ’ X H'c- H 'V X ■ • • , 


c’est-à-dire comme 

H m — n ’ ( *’-*- i i ‘+y l -+- ■•■) esta R m e~ H' 1 ( a'-t- ji'+ ■/+ . . .). 
11 est évident que les probabilités que H ou H' soient les 
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véritables valeurs de h, sont dans le môme rapport [ Theo- 
ria Motus Corporum cœlestium, art. 17 6 (*)]; par consé- 
quent la probabilité d’une valeui^quelconque de h est pro- 
portionnelle à 

h m e~ f 

et la valeur de h la plus probable est celle pour laquelle 
cette fonction devient un maximum. Mais on trouve par 
les règles connues que h est alors égal à 

\/ 2(a’+ p i -t~7 , + . . .) ’ 
donc la valeur de r la plus probable sera alors 



+ 7 7 + . • •) 

m 


OU 

0,6744897 X sJ 4 - P’ 4 - 7 ’ + -- •)■ 
Ce résultat est général, que m soit grand ou petit. 


4 . 


il est facile de comprendre que les valeurs trouvées 
pour h et pour r sont d’autant moins certaines que le nom- 
bre m est plus petit. 

Développons maintenant le degré d’exactitude que l’on 
doit attribuer aux valeurs de h et de r lorsque m est un 
nombre considérable. 

Désignons par H la valeur de h la plus probable que nous 
avons trouvée , 



et remarquons que la probabilité que H soit la véritable 
valeur de /», est à la probabilité que (H -|- X) soit cette véri- 


(*) Voyez pîigo. 1 16 de ce volume. 


J. B. 
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table valeur dans le rapport de 

m 

H . <? a à (H -f* /)". c 


îH’ ’ 


■ou comme 


i*m / i_ l i i' _i_ )' \ 

t . e H* \ 3‘0^4'H* 5 H <h ) 

Le second terme ne sera sensible par rapport au premier 
que si — est une petite fraction , et, dans ce cas > nous pour- 
rons remplacer le rapport indiqué par 

ï* m , 

.• . • TF- 

Ce qui veut dire : La probabilité que la valeur véritable 
de h soit comprise entre (II 4- ).) et (H 4 - 1 dl) est ap- 
proximativement égale à 

m 


K e 




où R est une constante telle, que l’intégrale 




d \ 


prise entre les limites admissibles de X, devienne égale 
à l’unité. 

Comme dans le cas actuel , à cause de la grande valeur 

A* ni 

de m , e 11 * devient excessivement petit lorsque — cesse 

d’ètre une petite fraction, il sera permis de prendre l’in- 
tégrale depuis — oo jusqu’à + oo , et l’on obtient 


i . fïîi * 

k “hV7 


Par conséquent, la probabilité que la véritable valeur de h 

G. lO 
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soit comprise entre (H — X) et (H -f- À), sera égale à 


U*'"]' 


et sera égale à ^ lorsque ^ \fm = p. 

11 y a donc un contre un à parier que la véritable valeur 
de h soit entre 

ou que la véritable valeur de r soit entre , 

- R R 

et — » 

*,_JL I + _P_ 

VÆ s/m 

où R désigne la valeur la plus probable de /■ trouvée dans 
l’article précédent. Ces limites peuvent s’appeler les limites 
' probables des véritables valeurs de h et de r. 11 est évi- 
dent que nous pourrons admettre ici pour limites probables 

R (’~i=) " R ( ,+ ;fe> : 


5 . 

Dans la discussion précédente, nous avons considéré a , 
(3, y, etc., comme des quantités définies et données, afin d’é- 
valuer la probabilité que la véritable valeur de h ou de r 
soit comprise entre certaines limites. 

On peut envisager la question sous un autre point de 
vue, en admettant que les erreurs des observations soient 
soumises à une loi /lé terminée de probabilité; on peut alors 
évaluer la probabilité pour que la somme des carrés de m 
erreurs d’observations tombe entre certaines limites. I.a- 
place a déjà résolu ce problème dans le cas où m est un 
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nombre très-grand, ainsi que le problème de déterminer 
la probabilité que la somme de m erreurs d’observations 
tombe entre certaines limites. 

Il est facile de généraliser cette recherche; je me bor-, 
lierai à indiquer ici le résultat. 

Désignons par ç (x) la probabilité d’une erreur d’obser- 
vation a:, de manière que 


/: 


- x . dx = I . 


Désignons enc.ore par K„ la valeur de l’intégrale 

r 

! tpx .x" dx. 

J — eo 

Soit ensuite 


S„= =c"-t- P" + 7 " 4 -. . 

où a, (3,y, etc., représentent m erreurs quelconques d’ob- 
servation; les termes de cette somme seront tous pris posi- 
tivement, même si n est impair. 

mK„ sera alors la valeur la plus probable de S„, et la 
probabilité que la véritable valeur de S„ tombe entre les 
limites (mK„ — et (mK„-|-l) sera égale à 

X 

par conséquent, les limites probables de S„ seront 


et 


m K„— p — K’) 

• • 

mK„-t- p y/ 2 m(K,„ — K’). 


Ce résultat s’applique , d’une manière générale, à toute 
loi des erreurs d’observation. En l’appliquant au cas parti- 
culier où 


« — h' x* 

<fX=---C , 

vV 
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nous trouverons 

— ■- - m, p = ? 

A" \Jt 

où le signe caractéristique II est pris dans la signification 
des Disqnisitiortcs generales circa serieni infini tam {*) 
( Comm . not>. Soc. Gotting., tome III; M. 5, art. 28.) 
Ainsi 

i 


K =1, 


K,= — , 

2 A’ 

K - '-1 

K * 4 a*’ 


K c = 


i.3.5 

8 A* ’ 


K,= 
K,= 
K s = 
K : = 


A \Jv 

I 

A 3 fiîr 

I . 2 
A 3 y/îr 
1 . 2.3 
A ’ 


par conséquent la valeur la plus probable de S„ sera 

mil ÿ (n — i) 

A" v ^ 

et les limites les plus probables de la véritable valeur de S„ 


seront 


mn ~ (, . A [~ n (”-T) • v'* Tjj 

A* | P Vm-[[n|(» —•!)]’ J) 


et 


m II - 


A" 


i + p 


. 

/» r n (”— t) • ^ _ ~i | 

V ; [(-rF Jj 


( * ) D’après la notatinn adoptée par tin grand nombre de géomètres, on 
aurait 

Il (n) =. r ( « i\ 

J. B. 
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Doue, si nous posons comme ci-dessus 


i=" 


où r représente l’erreur probable d'observation , la valeur 
la plus probable de 

« / S. \/it 

V-^) . 

sera évidemment r; et les limites probables de la valeur de 
cette quantité seront 

( P /a [~n( 'i — i) . \Jl ~| ) 

I “V LW)’ " Jl 

et 

/ip(— i. 

i J) 

Il y a donc aussi un contre un à parier que r soit compris 
entre les limites 
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ce qui s’accorde parfaitement avec celles que nous avons 
trouvées art. 4. 

En général on aura, si n est pair, les limites 

y »i. i .3.5.,. . (n — i) 


et 


f* V^TT 


s„ 


.3.5.7. ••(«-.) 
y I . . p . /» n» + i).(»+3)...|M-i) n j 

( «y/«L i.3.5. ..(n — 1 ) J ) ’ 

cl si n est impair, les limites suivantes : 

n J S„ \[v 

y 2 - 3 - - (^r 

• x } , _ £ x /lT , J.- 5 j . 7^( 2n - I J n _ J j 

\ n\ m | [2.4.6 (" — t)] ! J) 


et 


« j S „ . SJi 

\ / m. 1.2.3 . . 




Xi H-- 


p . / ifî .3.5.7... (« — «)*■ 


«y [ï.4-6. . . («— i)J' 

6 . 


ôl- 


J’ajoute ici encore les valeurs numériques pour les cas.les 
plus simples : 
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Limites probables de r. 

4 

'• °~>4 

H- 

,v. .,5,.s.,,xy^;(, T , '^)i 


V. o, 465 553a X ■ ( ■ qz - 0 ’ G3 ^ ) ; 


« /S t / o , r 

V ™ ■ v ,=f: 


VI. 0,429 497 a X(7^.[i+ ^7^7.7^ \ . 

y/n / 


U11 voit encore, par cetle comparaison, que la deuxième 
manière de déterminer r est la plus avantageuse; car cent 
erreurs d’observation, traitées d’après cette formule, don- 
nent un résultat aussi sur que 

1 >4 erreurs d’après la formule 1 ; 

109 » » » » III; 

1 33 » » » * » IV ; 

1 78 # » » u » V ; 

201 » » H H VI. 

Cependant la formule 1 présente l’avauiage de se prêter le 
mieux au calcul numérique; et comme son degré d’exacti- 
tude est peu inférieur à celui de la formule II , on peut tou- 
jours s’en servir, à moins que l’on ne connaisse déjà la 
somme des carrés des erreurs ou qu’on désire la connaître. 

1, 

Le procédé suivant est encore plus commode, mais beau- 
coup moins exact. 
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Ordonnons par ordre de grandeur les valeurs absolues 
des 7/1 erreurs d’observation 1 , et désignons par M le terme 
du milieu si leur nombre est impair, ou la moyenne arith- 
métique entre les deux termes du milieu si leur nombre 
' est pair. 

On peut démontrer (ce que nous ne ferons pas ici) que, 
pour un grand nombre d’observations, la valeur la plus 
probable de M est /', et que les limites probables de M 
sont • 

et 

ou que les limites probables de la valeur de r sont 
et 

' - a ('~'‘vÆ> . ' . 

ou, eu nombres, 

• «('T- 

Ce procédé n’est, par conséquent, pas beaucoup moins 
exact qui: l’application de la formule VI, et il faudrait 
mettre en compte 249 erreurs d’observation pour obtenir 
un résultat aussi certain qu’en appliquant la formule II à 
cent erreurs d’observation. 

8 . 

L’application de ces méthodes aux erreurs commises dans 
48 observations des ascensions droites de l’étoile Polaire, 


7 5au 974 \ 
\Jm ) 
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par bcsscl ('*'), a donné 

S, = 60" , 46 , 

Sj = 1 1 o" , 600 , 

S, = 250”, 34 11 18. 

De là on a déduit les valeurs les plus probables de r : 


D’après la formule I. 

i",o65. 

Erreur probable =dbo",o68. 

II. 

i", 024 , 

» =±o",o , jo. 

» III. 

i",ooi , 

u =±o", 072 , 

et d’après' l’article 7. 

i",o45, 

» = ± o", 1 1 3 : 


concordance de résultats qu’on pouvait à peine espérer. 
Bessel donne i",o6y et semble, .par conséquent, avoir cal- 
culé d’après la formule I. 

NOTE IV. 

APPLICATION DU CALCUL DES PROBABILITÉS A UN PROBLÈME 
DE GÉOMÉTRIE PRATIQUE. 

(Extrait d’une Lettre de M. Gauss à M. Schumacher; Astrononüsche 
Nachrichten , tome l , page 80. ) 

â_ 

t ' 

D’après votre désir, je vous envoie les règles relatives à 
l’emploi de la méthode des moindres carrés dans la solution 
du problème suivant : 

Déterminer lu position (C un point d'après les angles 
horizontaux observés de ce point entre d’autres points 
exactement connus. 

Cette question, très-élémentaire, ne peut embarrasser 
ceux qui ont bien saisi l’esprit de la méthode des moindres 
carrés. Je développerai néanmoins les formules auxquelles 


( * ) Bodf. , Annuaire astronomique pour 1818 , page lZf\. 
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elle conduit, en faveur des personnes qui auraient à traiter 
la question pratique sans avoir été à même d’étudier cette 
théorie. 

Soient a et b les coordonnées de l'un. des points donnés ; 
nous supposerons, que l’on compte les x positifs du nord 
au sud et les y positifs de l’ouest à l’est; soient x et y les 
coordonnées approchées du point inconnu et dx, dy les 
corrections encore inconnues qu’il faut leur attribuer. Dé- 
terminons deux quantités çp et r par les formules 

b— y a — x b — r 

tango = > r — : = — — —, 

a — x cos y sin f 

çp étant pris dans un tel quadrant , que la valeur de r sou 
positive. 

Posons de plus 

2 o 6 ?. 65 " ( b — y) 206266" (a — x) 

r ï ^ * r"* 


L’azimut du premier point est alors, pour un observateur 
placé au second (en prenant o pour azimut d’une parallèle à 
l’axe des x ), 

? 4- 4- p((r, 


les deux derniers termes étant exprimés en secondes. 

Soient ç?', /3' les quantités analogues à qp, «, (3 et re- 

latives au deuxième des points donnés, çp", a", (j" celles qui 
se rapportent au troisième, et ainsi de suite. Supposons 
<[ue, pour les mesures angulaires prises au point dont la 
position est inconnue, on ait fait usage d’un théodolite sans 
répétition, dont on dirige successivement la lunette vers 
les points connus , sans changer la place de l’instrument 
lui-mème. Si h , h' , h" sont les azimuts observés, on aurait , 
en supposant les observations rigoureusement exactes, et 
dx, dy exactement connus, 


y — i + 7 . dx + pf ly =. y' — h' -f- a! dx 4- dy 

= <p" — h" -f- %"dx ■+■ [ï’dy 
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Si donc on écrit que trois de ees différences ont la même 
valeur, on trouvera des valeurs approchées pour dx et dj ; 
si l’on n’a observé que trois points, -il n’y a rien de plus à 
faire; mais si le nombre des points considérés est plus 
grand, les erreurs seront le mieux compensées, en prenant 
la moyenne de ces diverses expressions (1), égalant à zéro 
la diilérence entre chacune d’elles et cette moyenne, et 
appliquant à ces équations la méthode des moindres carrés. 

Si toutes les mesures sont indépendantes les unes des 
autres, chacune d’elles fournit une équation entre dx et dy, 
et il faut combiner ces équations par la méthode des moin- 
dres carrés, en ayant égard, si l’on veut, à l’inégale préci- 
sion des observations. 

Soient, par exemple, i l’angle compris entre le premier et 
le deuxième point, i' l’angle compris entre le second et le 
troisième, et ainsi de suite, en comptant toujours de gauche 
à droite; on aura les équations 

f' — f — '! + ( a! — a ) dx -|- (£' — p)dy=zO, 

y" — ç' — /'+(«" — a!)dx -+- (fi" — p') dy — o. 

* 

Si les diverses mesures ont le même poids J on déduira de 
ces équations deux équations normales, en les ajoutant, 
après avoir successivement multiplié chacune d’elles par le 
coefficient de dx ou par celui de dy. 

Si, au contraire, les mesures des angles sont.d’une exac- 
titude inégale, et, par exemple, la première soit fondée 
sur p et la seconde sur p' répétitions , il faut que, dans les • 
deux cas, et avant leur addition, les équations soient encore 
multipliéespar p, p', etc. ; on trouve ensuite dx , dy , etc., 
par l'élimination entre les deux équations normales ainsi 
trouvées. 

(Les préceptes qui précèdent sont seulement destinés aux 
personnes auxquelles la méthode des moindres carrés est 
encore inconnue et pour lesquelles il sera peut-être bon de 
rappeler que, dans les multiplications, les signes de a ' — «, 
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(5' — (3, etc., doivent être rigoureusement conservés. Enfin, 
je remarque encore que l’on a seulement en vue de com- 
penser les erreurs commises sur les angles , les coordonnées 
des points donnés étant supposées exactes. ) 

Appliquons les indications qui précèdent aux- observa- 
tions que nous avons faites an commun sur le Holkensbas- 
tion, à Copenhague. Je dois prévenir que les résultats ne 
peuvent être ici d’uire exactitude rigoureuse. Les points ob- 
servés étant très-rapprochés de la station, une inexactitude 
de quelques dixièmes de pied sur leur position peut exercer 
une i nfluence beaucoup plus grande que les' erreurs habi- 
tuellement à craindre sur les angles mesurés. On ne s’éton- 
nera donc pas que la meilleure compensation des angles 
laisse subsister des différences beaucoup plus grandes que 
celles que l’on peut admettre comme possibles dans les ob- 
servations de même nature. Cette application doit être prise 
comme un exemple de la marche à suivre dans d’autres cas. 


Angles mesurés du Holkcnsbastion. 


F'riedrichsberg — Pétri 73° 35 ' 22.", 8; 

Pétri — Êrlosersthurm 104.57 . 33 ", o; 

Erlosersthurm — Friedrichsberg . 181.27. 5 ",o; 

Friedrichsberg — Fraueuthurm. . 80.37. 10", 8; 

Frauenthurm — Friedrichsthurm. iot . 1 1 . 5 o" ,8; 
Friedrichsthurm — Friedrichsberg. 178.11. i", 5 . 


Coordonnées des divers points en pieds de Paris , l'origine étant 
à l’Observatoire de Copenhague. 

Pétri 4 ^ 7 1 7 -F 1007,7; 

Frauenthurm ..... 710,0 -t- 684 , 2 ; 
Friedrichsberg . . . . 243 o ,6 -t- 8335 ,o; 

Erlosersthurm. .... '2940,0 — 3536 , o ; 
Friedrichsthurm. . . 3o59,3 — 223 i ,2. 


Les coordonnées approchées du Bastion sont 

x = -+- 2836,44 I 
x— -f- 444 ) 33 , 
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et nous trouvons ainsi les azimuts : 


Pétri i66° 3o' fa" , 5 & -t- 19, 92 rfx 4- 83, o4 rfy; 

Fruenthurm 173.33. 5o" , 54 4- 1 o , 80 rfx 4- 95 , 78 rfy ; 

Friedrichsberg 92 .56. 29", 46 4- 26,07 d J 4- 1, 34 dy; 


Erlosersthurin 27 1 . 29 . 25" , 38 — 5 1 , 79 rfx — 1 , 35 rfy ; 

Friedrichsthurm . . . 274. 45 ■4 ,W »4^ — 76,56dx — 6,38rfy. 

L’angle soiis lequel on voit la distance de Pétri à Friedrichs- 
berg est, par suite, 

73°34'3",io — 6,i5rfx4- 81.7 orf/; 
en l’égalant à l’angle observé, on a 

— 79", 70 — 6, i5 rfx 4- 8 1,70 rfy = o. - 

On' obtient, de la même manière, les équations sui- 
vantes : 

6g" ,82 — 7 1 , 7 1 rfx — 84 , 3g rfy = o , 

9", 08 4- 77 ,86 rfx 4- 2, 69 rfy =0, 
o",a8 — 15,27 rf x 4- g4 » 44 rfy = o, 

o",o4 — 87,36rfx — io2,i6rfy=o, 

— 3", 4® 4-102, 63rfx 4- 7,72rfy = o. 

En supposant les observations également précises, 011 déduit 
de là les équations normales 

29640 rfx 4- 1 4o33 rfy = 4 1 58", 
i4o33rfx 4-332igrfy= 12383", 

et, par suite, 

rfx = — o,o5, rfy =0,40; * 

les coordonnées du Bastion sont donc 
2836 , 3g , 

444,73. 

Les différences entre les valeurs observées des angles et 
celles que l’on calculerait d’après ces résultats sont trop 
grandes pour qu’on puisse les attribuer aux erreurs d’ob- 
servation ; elles tiennent , comme nous l’avons dit , à un dé- 
faut-de précision dans la détermination des points connus. 
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Les coordonnées x et y , adoptées comme premières ap- 
proximations, ont été déduites directement du quatrième 
et du cinquième angle. Quoique la méthode directe doive 
être considérée comme un sujet presque épuisé, j’indiquerai 
néanmoius , pour compléter, la méthode que j’ai l’habitude 
d’employer en pareil cas. 

Soient a et J les coordonnées du premier point connu; 
celles du second seront de la forme 

a -t- R cos E , b -f- R sin E , 
et celles du troisième, 

- - a R' cos E', b -t- R' sin E'. 

Soient 

a - f- p cos f , b -t- p sin s 

les coordonnées cherchées du point d’où l'on observe; soient 
M l’angle observé {toujours de la gauche À la droite) entre 
le premier et le deuxième point, et M' l’angle observé entre 
le premier et le troisième (en supposant qu’on en ait, s’il y 
a lieu, retranché 180 degrés); soient 

R _ R' _ , 

sin M n ’ sin M' n ’ 

E— M = N, E' — M' = N' : 

ou a les deux équations 

p — n sin (s — N ) , 
p =n' sin(e — N') : • 

qui, écrites de la manière suivante, 

P • 

n sin ( e — N ) ’ 

sin (t — N')* 

se résoudront par la méthode exposée ( Theotia Motus 
Corpofum cœlestium , page 82). 

L’une des solutions exposées en eet endroit conduit à la 
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règle suivante. Soit n' plus grand, ou au moins pas (dus 
petit que n , ce qui est évidemment permis, puisqu’on peut 
choisir le second point arbitrairement; posons 


p = tangÇ, 

tang 
tang 

on aura 

e étant connu , l’une des équations , ou même encore toutes 
deux, fourniront la valeur de p. 

Dans notre exemple, si nous considérons Frauen- 
thurm comme le premier point, Fricdriclisberg comme le 
deuxième et Friedrichsthurm comme le troisième, nous 
aurons : 

« = 710,0, b = 684,2, ' 

E = 77 0 19' 3 i" ,9?., E' = 3 o 8 ” 5 i' 45 ", 77 , 

log H = 3 , 8944 2 °5 , logR' = 3,5733549, 

M = 9g 0 22' 5 o" , 20 , M' = 101° 1 t' 5 o", 8 o, ' 

N = 337°56'42'',72, N' =207° 3 g. 54 ", 97 , 

logn = 3,9002650, logn' = 3,5817019. 

n' étant plus grand que », nous changerons ici l’ordre et 
nous poserons 


4 (N' — N) 
( 45 ° -ç) 


= tartg-s}. : 


N = 207 ”39' 54 " ,97 , N' = 337° 56 ' 4 a", 72 , 
logn =• 3,5817019, log n' — 3,goo265o; ' 


on en déduit 

> 9 ° 3 9 ' 3 ", 87, 
= 80. 45 . 3 . ",69, 

6 =353.33.5o",5 - 3, 
log p = 3 , 3303990 ; 

et les coordonnées du Holkensbastion , 

2836,44* > 444,33°. 
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NOTE V. 

SUR LA DÉTERMINATION CHRONOMÉTRIQUE DES LONGITUDES. 

( Astronomischti Nachrichlen , tome V, page 227.) 


Soient 0, 0', 6", etc., les époques (au nombre de n), 
auxquelles un chronomètre a accusé, les différences n, a ', 
a"., etc., avec les temps des lieux dont les longitudes sont 
x, x', x ff , etc., 0, 0',0", etc., étant supposées réduites au 
temps d’un seul et même lieu et u désignant l’avance jour- 
nalière du chronomètre ; on aurait, si l’instrument était 
parfaitement régulier, les équations 

a — G « — x = a! — G' it — x' — a" — G " u — x" = . ... 

Pour que ces équations suffisent à la détermination dés in- 
connues x, x', x ",... , u, il faut, d’une part, considérer l’une 
des longitudes comme donnée, et, d’autre part, il est néces- 
< saire que l’on ait observé au moins deux fois dans le même 
lieu; de telle sorte que deux au moins des inconnues x, x', 
x w , etc., soient égales entre elles. Si parmi ces quantités il 
n’y en a que deux qui soient identiques, le problème est 
entièrement déterminé; dans le cas contraire, il devient 
indéterminé , et l’on doit faire en sorte que les équations 

o = « — a' - 1 - ( G' — G ) u — x x’, 
o = -t- ( G" — G' ) « — x' x", 

o— a"— a'"- 1- ( G"- G") « — x"+ x ", 


soient satisfaites aussi exactement que possible, car les 
imperfections inévitables du chronomètre ne permettront 
jamais de satisfaire rigoureusement à toutes. Mais on ne 
doit pas accorder à ces équations des poids égaux; car les 
quantités 

a — a' + (G' — G) u — i+/, 
a' — n” ( G" — G' ) u — x' -f- x" , 
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représentent l’accumulation de tous les écarts de marche des 
chronomètres dans les intervalles 0' — 0, 0" — 0', etc., et 
s’il s’agit d'un bon chronomètre, à qui l’on puisse réelle- 
ment attribuer une marche moyenne sans variation crois- 
sante dans un seul sens, la valeur moyenne à craindre pour 
une pareille somme peut être considérée comme propor- 
tionnelle à la racine carrée du temps écoulé. 

On devra donc, dans l’application de la méthode des 
moindres carrés, regarder les équations précédentes comme 
ayant des poids inversement proportionnels aux dillcrences 
6'— 0, 0"—V, Q'"—0", etc. 

La solution n’a alors aucune difficulté, et fournira les 
valeurs les plus probables de x , x', x", etc., ainsi que le 
poids de chaque détermination. 

J’ajouterai pourtant quelques remarques. 

I. Si la première et la dernière observation ont été faites 
au même lieu , la valeur la plus plausible de u est celle qui 
résulte simplement de la comparaison de ces observations 
extrêmes. Le calcul devient alors très-simple, car, en vertu 
d’un théorème bien facile à démontrer, on peut, dans les 
équations, remplacer u par sa valeur la plus plausible, 
ou, ce qui revient au même, on peut employer cette valeur 
supposée exacte pour corriger les observations et les rame- 
ner à ce qu’elles seraient avec un chronomètre fictif dont 
1 avance serait nulle. 

II. Si l’on a attribué simplement aux diverses équations 
des poids égaux à 

i i i 

V — V — g" ’ ’ ' ‘ ’ 

l’unité de précision pour les poids obtenus sera l’exactitude 
de celle que l’on obtiendrait à l’aide du même chronomètre 
observé deux fois seulement , et à un jour d'intervalle; mais 
pour pouvoir comparer les résultats obtenus à l’aide de 
divers chronomètres inégalement précis, il faut encore in- 
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troduire clans le résultat un facteur qui dépende de la plus 
ou de la moins grande perfection de chaque chronomètre 
employé. 

Pour y parvenir, je suppose que les expressions 
n — a' - (— (.©' — 0 ) u — x -t- x', 
a' — a" + ( 0" — 0') u — x' -h x", 

deviennent respectivement A, A', A", etc., quand on substi- 
tue aux inconnues leurs valeurs les plus plausibles, soit 

V X" X " 1 _ „ 

0' — 0 + 0"— 0' ~ l " 0'"— 0" S ’ 

si v est le nombre des inconnues et que l’on pose 
y n — v — i 

le facteur spécifique relatif à chaque chronomètre est pro- 
portionnel à ou à - — ^ 1 et l’on peut regarder m 

comme l’écart de la marche moyenne qui est à craindre 
pour la marche d’une journée. 

m. i .es règles qui précèdent sont relatives à un chrono- 
mètre dont la marche n’est soumise à aucune irrégularité 
sensible et croissante avec le temps. Si cette hypothèse n’é- 
tait pas permise, on pourrait supposer, lorsque les observa- 
tions n’embrassent qu’une période qui n’est pas excessive- 
ment considérable, une variation proportionnelle au temps 
dans l’avance journalière de l’instrument, en introduisant 
ainsi une inconnue de plus. 

Les équations prendraient alors la forme suivante : 

o == a — n' -4- (0’ — 0) u -+- (0' 3 — 0’j v — x -y- x', 


” IV. Pour ce qui concerne la résolution des équations 
d’après la méthode des moindres carrés, il n’est peut-être pas 
inutile de rappeler qu’on doit commencer, dans le plus grand 
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nombre îles cas , par calculer les valeurs approchées îles in- 
connues et appliquer ensuite la méthode à la détermination 
des corrections très-petites que les valeurs doivent subir. 

Il nous a paru utile de rappeler ce conseil général , parce 
que beaucoup de calculateurs ont paru l’oublier, et ont été 
conduits à des calculs plus laborieux et peut-être moins 
exacts . 

J’ai déterminé la marche des cinq chronomètres suivants : 




1 


llKFt.lET 

BARRALD 

KESSEL 




3o56. 

9°4- 

1352. 

Greenwich, 3o Juin 
j 5 Juill 

h 1 

3.22 

3. l5 

J // 

- 8., 7,, 4 
10.44,39 

t U 

1 ■ V*7 
i.3a,i5 

/ u 

+3o . 5;> , 76 

/ // 

+48.29,20 

_ f H 
-+-5o. 29, 3 1 

.38 . 

3. i3 

11. n,G9 

1.36,96 

3o. 60,07 

4« .40,24 

50.39,69 

2 Août 

1 . 1 5 

1 1 .28,48 

• 44 44 

30.31,78 

48.58,87 

5o .62,1 4 

17 • 

10.18 

■a. 59,40 

a. 6,24 

29.35,69 

4 9 .5 7 ,83 

5i .38,66 

25 » 

7.37 

"3. 47.9# 

a. 1 5 , 84 

29. io,48 

5 o. 27, 

5a. a, 45 

10 Sept. 

7.40 

i5.a4.47 

a . 4o', 36 




Helgoland , 3 Juill 

3.40 

— 40. 8,00 

30.26,84 




23 » 

13.40 

/|3. 2,03 

3o. 3,89 

— -.0.20 , 34 

+16.47,39 

-f-l8.48, 39 

5 Août 

.-48 

/|3.l8 ; ll 

39.43,35 

1 . 10,24 

17.37,51 

en 

1 1 B 

i3- 9 

43.35,77 

39.33,4s 

1 .32,7-5 

18. i,3o 

Eh 

« 3o » 

19. 3o 

45.53,08 

39. 7,96 

3.40,67 

19.17,03 

S; BjI 

6 Sept. 

3. (i 

46. 5 i ,5G 

38.58,94 

3. 4,55 

19.43,80 

ai 6,56 

7 » 

8.4a 

46.38,73 

28.56,71 




Alloua, G Août 

5 . 55 

— 5i .38, 9.5 

— 37.55,76 

9.38,50 

+ 9.38,48 

+1 1 . 16, a5 

9 • 

12.35 

5i .57,35 

37.50,03 

9.38,81" 

94o,3o 

11.37,76 

3i » 

907 

54 . 10,33 

37.21 ,3o 

10.56,68 

1t. 5,92 

i 3 . a5 , 96 

4 Sept. 

22.13 

54.39,16 

37. l5,21 

1 1 . 1 5 , 36 

u 24,49 

12.48,10 

Brême, i3 Août 

0.2 

—47.50,65 

-33.16,49 

— 5.23,37 

+ 14.31,86 

+ 16. 5,83 


Prenons pour exemple le chronomètre de Brcguet 3o56. 
Soient o la longitude de Helgoland, — x celle de Greenwich, 
y celle d’Allona. Je ne fais pas entrer celle de Brème dans 
le calcul, car n’ayant pour celle ville qu’une seule obser- 
vation, il est impossible de la contrôler. Je compte le temps 
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à partir tic la première comparaison du chronomètre n° x 
(Greenwich, juin, 3 h 22*'). En substituant au chrono- 
mètre de Breguet un instrument fictif dont l’avance soit 


nulle, nous trouvons : 

S 

22,4 

+ 

60,20; 

25,0 

-h 

lq4q,6o — x ; 

28,0 

H- 

1950,87 XJ 

3 a, 9 

■+■ 

1 q 5 o ,29 — x ; 

35,9 


5 q,o 8 ; 

37,1 

— 

434,98 H- y ; 

4 o ,4 

— 

433 , 4 q +- y ; 

4 a , 4 


5 q ,88 ; 

48,3 


hj49,6o — x; 

56,2 


1952,74 — x; 

61 ,6 

-+- 

8'. 3 a; 

62,2 

— 

432,53 -t- y ; 


66,8 

— 434,98 + y, 

68,0 

+ 60, iq. 


Dans les équations, ci-dessus, les inconnues x et y sont 
séparées, ce qui facilite leur détermination-, nous trouvons 
pour x quatre déterminations : 



POIDS. 

i88q,4o 

-4 = o, 38 ; 
2,0 

1891 ,21 

— = o, 33 ; 
3,o I 

1889,78 

à=°’ 17; 

1891 ,4a 

3 T 4 = °’' 9 ’ 

d’où l’on déduit : 

x == 1890", 36 

1 , 07 ; 

on trouve de même ; 

11 

& 

» 

3 , 83 . 
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D’après ces valeurs, le chronomètre fictif marquera, en 
temps de Helgoland, 


a 

A 


a 

60,20 ; 

25,0 

59,24 — 0,96; 

28.0 

60, 5i 4- 1 ,27 ; 

32,9 

59,93 — o,58 ; 

35,9 

59,08 — o,85 ; 

37,i 

59, i4 4- o, 06; 

4°. 4 

60 ,63 4- 1 ,47; 

42,4 

5q,88 — 0,75 ; 

48,3 

59,24 — 0,62; 

56.2 

62,38 4- 3, i4 ; 

61 ,6 

61 .32 — t ,06; 

62,2 

61 ,5q 4- 0,27 ; 

66,8 

5g» 1 4 — 2,45; 

68,0 

60,19+ 1 ,o5; 


d’où l’on déduit : 

S = 6,00, 



et l’erreur moyenne à craindre est : 

pour x, o",75 ; 

pour y , o",4°- 

L«js résultats fournis par les cinq chronomètres donnent: 



EilliEUR MOYENNE 

POIDS. 

Breguet 

x — 1 890 , 36 

à craindre. 

0,75 

t,j8 

Kessel 

1893,39 

0,67 

2,23 

Barraud . . . . 

1892,32 

. 2 ± 4 a 

4. ‘8 


1892,39 

1892,52 

î = 1 892 , 35 

0,43 

0,35 

5 , 4 i 

4 

8,16 

Moyenne. 


ai . 74 


Digitized by Google 


( i 66 ) 

ou trouva de même, poiir y : 


Brcguet y = 4 !)4 > 1 2 o, io 6,25 

Kessel 493,89 o ,36 7,72 

Barraud .... 4g3 ,67 0,21 1 4 > 79 

1 49 3 < 9 8 °- 2 9 "189 

4- 494' '8 o, «4 >7,36 

58,01 


Le nombre placé sous le nom de poids dans la dernière 
colonne, est l'inverse du carré de l’erreur moyenne à 
craindre, en prenant ainsi , pour poids unité, celui qui cor- 
respond aux observations qui donnent une erreur moyenne 
à craindre égale à 1 ", de sorte que, pour Altona, l’erreur 

1» 

moyenne à craindre est ■ ■ — = o ,,r ,i3; mais il vaut mieux 

v58,oi 

considérer les nombres de la dernière colonne comme indi- 
quant seulement des rapports , et déduire la précision ab- 
solue de la différence entre les valeurs des derniers résultats 
trouvés pour x et y au moyen de chaque chronomètre. La 
précision trouvée de cette manière sera toujours un peu trop 
grande, puisque les déterminations de temps à Greenwich, 
à Helgoland et à Altona n’ont pas une précision absolue, 
et que, par conséquent, quel que soit le nombre des chro- 
nomètres, les erreurs qui proviennent de celte source se 
feront toujours sentir dans chaque résultat final. 

On peut obtenir, de la manière suivante, la longitude 
de Brème. 

Soit z cette longitude à l’est de Helgoland : la comparai- 
son du chronomètre de Breguet donne la position du chro- 
nomètre fictif, 

— i 65",52 2, 

et l’on déduit de la comparaison avec les résultats précé- 
dents , 

POIDS. 


H 

Z = 225,4 


O 


>4 


= 0,7 ; 
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Breguet 225,24 i,5 

Kessel 225,84 1,9 

Barraud 225,3g 3,6 

1 226,04 2,9 

4 224,86 4>3 

225,4?- 4,? 


La longitude de Brême, qui serait, d’après cela*, 268", 54 
à l’ouest d’Altona, est naturellement affectée des erreurs 
sur la détermination du temps à Brème, et cette différence 
semble trop petite de plusieurs secondes. D ! après mes trian- 
gulations, la tour de Ansgarius est de 2y3 ,/ , 5 1 en temps «à 
l’ouest de Gottingue, et l'observatoire d’Olbers à 271", 9. 


FIN, 


Sg/v 

608271 
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